Populating Chemical Space with Peptides Using a Genetic Algorithm

Design of Cyclic Peptides Targeting Protein-Protein Interactions Using AlphaFold

More than 930,000 protein-protein interactions (PPIs) have been identified in recent years, but their physicochemical properties differ from conventional drug targets, complicating the use of conventional small molecules as modalities. Cyclic peptides are a promising modality for targeting PPIs, but it is difficult to predict the structure of a target protein-cyclic peptide complex or to design a cyclic peptide sequence that binds to the target protein using computational methods. Recently, AlphaFold with a cyclic offset has enabled predicting the structure of cyclic peptides, thereby enabling de novo cyclic peptide designs. We developed a cyclic peptide complex offset to enable the structural prediction of target proteins and cyclic peptide complexes and found AlphaFold2 with a cyclic peptide complex offset can predict structures with high accuracy. We also applied the cyclic peptide complex offset to the binder hallucination protocol of AfDesign, a de novo protein design method using AlphaFold, and we could design a high predicted local-distance difference test and lower separated binding energy per unit interface area than the native MDM2/p53 structure. Furthermore, the method was applied to 12 other protein-peptide complexes and one protein-protein complex. Our approach shows that it is possible to design putative cyclic peptide sequences targeting PPI.

Machine Learning Guided Discovery of Non-Hemolytic Membrane Disruptive Anti-Cancer Peptides

Most antimicrobial peptides (AMPs) and anticancer peptides (ACPs) fold into membrane disruptive cationic amphiphilic α‐helices, many of which are however also unpredictably hemolytic and toxic. Here we exploited the ability of recurrent neural networks (RNN) to distinguish active from inactive and non‐hemolytic from hemolytic AMPs and ACPs to discover new non‐hemolytic ACPs. Our discovery pipeline involved: 1) sequence generation using either a generative RNN or a genetic algorithm, 2) RNN classification for activity and hemolysis, 3) selection for sequence novelty, helicity and amphiphilicity, and 4) synthesis and testing. Experimental evaluation of thirty‐three peptides resulted in eleven active ACPs, four of which were non‐hemolytic, with properties resembling those of the natural ACP lasioglossin III. These experiments show the first example of direct machine learning guided discovery of non‐hemolytic ACPs.

Solubility-Aware Protein Binding Peptide Design Using AlphaFold

New protein–protein interactions (PPIs) are identified, but PPIs have different physicochemical properties compared with conventional targets, making it difficult to use small molecules. Peptides offer a new modality to target PPIs, but designing appropriate peptide sequences by computation is challenging. Recently, AlphaFold and RoseTTAFold have made it possible to predict protein structures from amino acid sequences with ultra-high accuracy, enabling de novo protein design. We designed peptides likely to have PPI as the target protein using the “binder hallucination” protocol of AfDesign, a de novo protein design method using AlphaFold. However, the solubility of the peptides tended to be low. Therefore, we designed a solubility loss function using solubility indices for amino acids and developed a solubility-aware AfDesign binder hallucination protocol. The peptide solubility in sequences designed using the new protocol increased with the weight of the solubility loss function; moreover, they captured the characteristics of the solubility indices. Moreover, the new protocol sequences tended to have higher affinity than random or single residue substitution sequences when evaluated by docking binding affinity. Our approach shows that it is possible to design peptide sequences that can bind to the interface of PPI while controlling solubility.

Approaches for enhancing the analysis of chemical space for drug discovery

INTRODUCTION

Chemical space is a powerful, general, and practical conceptual framework in drug discovery and other areas in chemistry that addresses the diversity of molecules and it has various applications. Moreover, chemical space is a cornerstone of chemoinformatics as a scientific discipline. In response to the increase in the set of chemical compounds in databases, generators of chemical structures, and tools to calculate molecular descriptors, novel approaches to generate visual representations of chemical space in low dimensions are emerging and evolving. Such approaches include a wide range of commercial and free applications, software, and open-source methods.

AREAS COVERED

The current state of chemical space in drug design and discovery is reviewed. The topics discussed herein include advances for efficient navigation in chemical space, the use of this concept in assessing the diversity of different data sets, exploring structure-property/activity relationships for one or multiple endpoints, and compound library design. Recent advances in methodologies for generating visual representations of chemical space have been highlighted, thereby emphasizing open-source methods.

EXPERT OPINION

Quantitative and qualitative generation and analysis of chemical space require novel approaches for handling the increasing number of molecules and their information available in chemical databases (including emerging ultra-large libraries). In addition, it is of utmost importance to note that chemical space is a conceptual framework that goes beyond visual representation in low dimensions. However, the graphical representation of chemical space has several practical applications in drug discovery and beyond.

Pareto Optimization of Oligomer Polarizability and Dipole Moment Using a Genetic Algorithm

High-performance electronic components are highly sought after in order to produce increasingly smaller and cheaper electronic devices. Drawing inspiration from inorganic dielectric materials, in which both polarizability and polarization contribute, organic materials can also maximize both. For a large set of small molecules drawn from PubChem, a Pareto-like front appears between the polarizability and dipole moment, indicating the presence of an apparent trade-off between these two properties. We tested this balance in π-conjugated materials by searching for novel conjugated hexamers with simultaneously large polarizabilities and dipole moments with potential use for dielectric materials. Using a genetic algorithm (GA) screening technique in conjunction with an approximate density functional tight-binding method for property calculations, we were able to efficiently search chemical space for optimal hexamers. Given the scope of chemical space, using the GA technique saves considerable time and resources by speeding up molecular searches compared to a systematic search. We also explored the underlying structure-function relationships, including sequence and monomer properties, that characterize large polarizability and dipole moment regimes.

PepSeA: Peptide Sequence Alignment and Visualization Tools to Enable Lead Optimization

Therapeutic peptides offer potential advantages over small molecules in terms of selectivity, affinity, and their ability to target "undruggable" proteins that are associated with a wide range of pathologies. Despite their importance, current molecular design capabilities that inform medicinal chemistry decisions on peptide programs are limited. More specifically, there are unmet needs for structure-activity relationship (SAR) analysis and visualization of linear, cyclic, and cross-linked peptides containing non-natural motifs, which are widely used in drug discovery. To bridge this gap, we developed PepSeA (Peptide Sequence Alignment and Visualization), an open-source, freely available package of sequence-based tools (https://github.com/Merck/PepSeA). PepSeA enables multiple sequence alignment of non-natural amino acids and enhanced visualization with the hierarchical editing language for macromolecules (HELM). Via stepwise SAR analysis of a ChEMBL peptide data set, we demonstrate the utility of PepSeA to accelerate decision making in lead optimization campaigns in pharmaceutical setting. PepSeA represents an initial attempt to expand cheminformatics capabilities for therapeutic peptides and to enable rapid and more efficient design-make-test cycles.

Machine learning & deep learning in data-driven decision making of drug discovery and challenges in high-quality data acquisition in the pharmaceutical industry

Predicting novel small molecule bioactivities for the target deconvolution, hit-to-lead optimization in drug discovery research, requires molecular representation. Previous reports have demonstrated that machine learning (ML) and deep learning (DL) have substantial implications in virtual screening, peptide synthesis, drug ADMET screening and biomarker discovery. These strategies can increase the positive outcomes in the drug discovery process without false-positive rates and can be achieved in a cost-effective way with a minimum duration of time by high-quality data acquisition. This review substantially discusses the recent updates in AI tools as cheminformatics application in medicinal chemistry for the data-driven decision making of drug discovery and challenges in high-quality data acquisition in the pharmaceutical industry while improving small-molecule bioactivities and properties.

Systematic Evaluation of Liquid Chromatography (LC) Column Combinations for Application in Two-Dimensional LC Metabolomic Studies

In comparison to proteomics, the application of two-dimensional liquid chromatography (2D LC) in the field of metabolomics is still premature. One reason might be the elevated chemical complexity and the associated challenge of selecting proper separation conditions in each dimension. As orthogonality of dimensions is a major issue, the present study aimed for the identification of successful stationary phase combinations. To determine the degree of orthogonality, first, six different metrics, namely, Pearson's correlation coefficient (1 - |R|), the nearest-neighbor distances (H̅_NND), the "asterisk equations" (A_O), and surface coverage by bins (SC_G), convex hulls (SC_CH), and α-convex hulls (SC_αH), were critically assessed by 15 artificial 2D data sets, and a systematic parameter optimization of α-convex hulls was conducted. SG_G, SC_αH with α = 0.1, and H̅_NND generated valid results with sensitivity toward space utilization and data distribution and, therefore, were applied to pairs of experimental retention time sets obtained for >350 metabolites, selected to represent the chemical space of human urine. Normalized retention data were obtained for 23 chromatographic setups, comprising reversed-phase (RP), hydrophilic interaction liquid chromatography (HILIC), and mixed-mode separation systems with an ion exchange (IEX) contribution. As expected, no single LC setting provided separation of all considered analytes, but while conventional RP×HILIC combinations appeared rather complementary than orthogonal, the incorporation of IEX properties into the RP dimension substantially increased the 2D potential. Eventually, one of the most promising column combinations was implemented for an offline 2D LC time-of-flight mass spectrometry analysis of a lyophilized urine sample. Targeted screening resulted in a total of 164 detected metabolites and confirmed the outstanding coverage of the 2D retention space.

Progress on open chemoinformatic tools for expanding and exploring the chemical space

The concept of chemical space is a cornerstone in chemoinformatics, and it has broad conceptual and practical applicability in many areas of chemistry, including drug design and discovery. One of the most considerable impacts is in the study of structure-property relationships where the property can be a biological activity or any other characteristic of interest to a particular chemistry discipline. The chemical space is highly dependent on the molecular representation that is also a cornerstone concept in computational chemistry. Herein, we discuss the recent progress on chemoinformatic tools developed to expand and characterize the chemical space of compound data sets using different types of molecular representations, generate visual representations of such spaces, and explore structure-property relationships in the context of chemical spaces. We emphasize the development of methods and freely available tools focusing on drug discovery applications. We also comment on the general advantages and shortcomings of using freely available and easy-to-use tools and discuss the value of using such open resources for research, education, and scientific dissemination.

Assigning the Origin of Microbial Natural Products by Chemical Space Map and Machine Learning

Microbial natural products (NPs) are an important source of drugs, however, their structural diversity remains poorly understood. Here we used our recently reported MinHashed Atom Pair fingerprint with diameter of four bonds (MAP4), a fingerprint suitable for molecules across very different sizes, to analyze the Natural Products Atlas (NPAtlas), a database of 25,523 NPs of bacterial or fungal origin. To visualize NPAtlas by MAP4 similarity, we used the dimensionality reduction method tree map (TMAP). The resulting interactive map organizes molecules by physico-chemical properties and compound families such as peptides and glycosides. Remarkably, the map separates bacterial and fungal NPs from one another, revealing that these two compound families are intrinsically different despite their related biosynthetic pathways. We used these differences to train a machine learning model capable of distinguishing between NPs of bacterial or fungal origin.

D-Peptide Builder: A Web Service to Enumerate, Analyze, and Visualize the Chemical Space of Combinatorial Peptide Libraries

Peptide-based drug discovery is re-gaining attention in drug discovery. Similarly, combinatorial chemistry continues to be a useful technique for the rapid exploration of chemical space. A current challenge, however, is the enumeration of combinatorial peptide libraries using freely accessible tools. To facilitate the swift enumeration of combinatorial peptide libraries, we introduce herein D-Peptide Builder. In the current version, the user can build up to pentapeptides, linear or cyclic, using the natural pool of 20 amino acids. The user can use non- and/or N-methylated amino acids. The server also enables the rapid visualization of the chemical space of the newly enumerated peptides in comparison with other libraries relevant to drug discovery and preloaded in the server. D-Peptide Builder is freely accessible at http://dpeptidebuilder. quimica.unam.mx:4000/. It is also accessible through the open D-Tools platform (DIFACQUIM Tools for Chemoinformatics https://www.difacquim.com/d-tools/).

One molecular fingerprint to rule them all: drugs, biomolecules, and the metabolome

Background

Molecular fingerprints are essential cheminformatics tools for virtual screening and mapping chemical space. Among the different types of fingerprints, substructure fingerprints perform best for small molecules such as drugs, while atom-pair fingerprints are preferable for large molecules such as peptides. However, no available fingerprint achieves good performance on both classes of molecules.

Results

Here we set out to design a new fingerprint suitable for both small and large molecules by combining substructure and atom-pair concepts. Our quest resulted in a new fingerprint called MinHashed atom-pair fingerprint up to a diameter of four bonds (MAP4). In this fingerprint the circular substructures with radii of r = 1 and r = 2 bonds around each atom in an atom-pair are written as two pairs of SMILES, each pair being combined with the topological distance separating the two central atoms. These so-called atom-pair molecular shingles are hashed, and the resulting set of hashes is MinHashed to form the MAP4 fingerprint. MAP4 significantly outperforms all other fingerprints on an extended benchmark that combines the Riniker and Landrum small molecule benchmark with a peptide benchmark recovering BLAST analogs from either scrambled or point mutation analogs. MAP4 furthermore produces well-organized chemical space tree-maps (TMAPs) for databases as diverse as DrugBank, ChEMBL, SwissProt and the Human Metabolome Database (HMBD), and differentiates between all metabolites in HMBD, over 70% of which are indistinguishable from their nearest neighbor using substructure fingerprints.

Conclusion

MAP4 is a new molecular fingerprint suitable for drugs, biomolecules, and the metabolome and can be adopted as a universal fingerprint to describe and search chemical space. The source code is available at https://github.com/reymond-group/map4 and interactive MAP4 similarity search tools and TMAPs for various databases are accessible at http://map-search.gdb.tools/ and http://tm.gdb.tools/map4/.

Wasp venom peptide as a new antichagasic agent

Chagas disease is caused by Trypanosoma cruzi and affects approximately 10 million people a year worldwide. The only two treatment options, benznidazole and nifurtimox, have low efficacy and high toxicity towards human cells. Mastoporan peptide (MP) a small cationic AMP from the venom of the wasp Polybia paulista has been reported as a potent trypanocidal agent. Thus, we evaluated the antichagasic effect of another AMP from the venom of the same wasp Polybia paulista, polybia-CP (ILGTILGLLSKL-NH₂), and investigated its mechanism of action against different stages of the trypanosomal cells life cycle. Polybia-CP was tested against the epimastigote, trypomastigote and amastigote forms of the T. cruzi Y strain (benznidazole-resistant strain) and inhibited the development of these forms. We also assessed the selectivity of the AMP against mammalian cells by exposing LLC-MK2 cells to polybia-CP, the peptide presented a high selectivity index (>106). The mechanism of action of polybia-CP on trypanosomal cells was investigated by flow cytometry, scanning electron microscopy (SEM) and enzymatic assays with T. cruzi GAPDH (tcGAPDH), enzyme that catalyzes the sixth step of glycolysis. Polybia-CP induced phosphatidylserine exposure, it also increased the formation of reactive species of oxigen (ROS) and reduced the transmembrane mitochondrial potential. Polybia-CP also led to cell shrinkage, evidencing apoptotic cell death. We did not observe the inhibition of tcGAPDH or autophagy induction. Altogether, polybia-CP has shown the features of a promising template for the development of new antichagasic agents.