Exploring the chemical space of peptides for drug discovery: a focus on linear and cyclic penta-peptides

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.

The Dimeric Peptide (KKYRYHLKPF)2K Shows Broad-Spectrum Antiviral Activity by Inhibiting Different Steps of Chikungunya and Zika Virus Infection

Chikungunya virus (CHIKV) and Zika virus (ZIKV) are important disease-causing agents worldwide. Currently, there are no antiviral drugs or vaccines approved to treat these viruses. However, peptides have shown great potential for new drug development. A recent study described (p-BthTX-I)₂K [(KKYRYHLKPF)₂K], a peptide derived from the Bothropstoxin-I toxin in the venom of the Bothrops jararacussu snake, showed antiviral activity against SARS-CoV-2. In this study, we assessed the activity of this peptide against CHIKV and ZIKV and its antiviral action in the different stages of the viral replication cycle in vitro. We observed that (p-BthTX-I)₂K impaired CHIKV infection by interfering with the early steps of the viral replication cycle, reducing CHIKV entry into BHK-21 cells specifically by reducing both the attachment and internalization steps. (p-BthTX-I)₂K also inhibited the ZIKV replicative cycle in Vero cells. The peptide protected the cells against ZIKV infection and decreased the levels of the viral RNA and the NS3 protein of this virus at viral post-entry steps. In conclusion, this study highlights the potential of the (p-BthTX-I)₂K peptide to be a novel broad-spectrum antiviral candidate that targets different steps of the replication cycle of both CHIKV and ZIKV.

In Silico Prediction of Anti-Infective and Cell-Penetrating Peptides from Thalassophryne nattereri Natterin Toxins

The therapeutic potential of venom-derived peptides, such as bioactive peptides (BAPs), is determined by specificity, stability, and pharmacokinetics properties. BAPs, including anti-infective or antimicrobial peptides (AMPs) and cell-penetrating peptides (CPPs), share several physicochemical characteristics and are potential alternatives to antibiotic-based therapies and drug delivery systems, respectively. This study used in silico methods to predict AMPs and CPPs derived from natterins from the venomous fish Thalassophryne nattereri. Fifty-seven BAPs (19 AMPs, 8 CPPs, and 30 AMPs/CPPs) were identified using the web servers CAMP, AMPA, AmpGram, C2Pred, and CellPPD. The physicochemical properties were analyzed using ProtParam, PepCalc, and DispHred tools. The membrane-binding potential and cellular location of each peptide were analyzed using the Boman index by APD3, and TMHMM web servers. All CPPs and two AMPs showed high membrane-binding potential. Fifty-four peptides were located in the plasma membrane. Peptide immunogenicity, toxicity, allergenicity, and ADMET parameters were evaluated using several web servers. Sixteen antiviral peptides and 37 anticancer peptides were predicted using the web servers Meta-iAVP and ACPred. Secondary structures and helical wheel projections were predicted using the PEP-FOLD3 and Heliquest web servers. Fifteen peptides are potential lead compounds and were selected to be further synthesized and tested experimentally in vitro to validate the in silico screening. The use of computer-aided design for predicting peptide structure and activity is fast and cost-effective and facilitates the design of potent therapeutic peptides. The results demonstrate that toxins form a natural biotechnological platform in drug discovery, and the presence of CPP and AMP sequences in toxin families opens new possibilities in toxin biochemistry research.

Biological Membrane-Penetrating Peptides: Computational Prediction and Applications

Peptides comprise a versatile class of biomolecules that present a unique chemical space with diverse physicochemical and structural properties. Some classes of peptides are able to naturally cross the biological membranes, such as cell membrane and blood-brain barrier (BBB). Cell-penetrating peptides (CPPs) and blood-brain barrier-penetrating peptides (B3PPs) have been explored by the biotechnological and pharmaceutical industries to develop new therapeutic molecules and carrier systems. The computational prediction of peptides’ penetration into biological membranes has been emerged as an interesting strategy due to their high throughput and low-cost screening of large chemical libraries. Structure- and sequence-based information of peptides, as well as atomistic biophysical models, have been explored in computer-assisted discovery strategies to classify and identify new structures with pharmacokinetic properties related to the translocation through biomembranes. Computational strategies to predict the permeability into biomembranes include cheminformatic filters, molecular dynamics simulations, artificial intelligence algorithms, and statistical models, and the choice of the most adequate method depends on the purposes of the computational investigation. Here, we exhibit and discuss some principles and applications of these computational methods widely used to predict the permeability of peptides into biomembranes, exhibiting some of their pharmaceutical and biotechnological applications.

Latin American databases of natural products: biodiversity and drug discovery against SARS-CoV-2

In this study, we evaluated 3444 Latin American natural products using cheminformatic tools. We also characterized 196 compounds for the first time from the flora of El Salvador that were compared with the databases of secondary metabolites from Brazil, Mexico, and Panama, and 42 969 compounds (natural, semi-synthetic, synthetic) from different regions of the world. The overall analysis was performed using drug-likeness properties, molecular fingerprints of different designs, two parameters similarity, molecular scaffolds, and molecular complexity metrics. It was found that, in general, Salvadoran natural products have a large diversity based on fingerprints. Simultaneously, those belonging to Mexico and Panama present the greatest diversity of scaffolds compared to the other databases. This study provided evidence of the high structural complexity that Latin America's natural products have as a benchmark. The COVID-19 pandemic has had a negative effect on a global level. Thus, in the search for substances that may influence the coronavirus life cycle, the secondary metabolites from El Salvador and Panama were evaluated by docking against the endoribonuclease NSP-15, an enzyme involved in the SARS CoV-2 viral replication. We propose in this study three natural products as potential inhibitors of NSP-15.

Predicting cell-penetrating peptides using machine learning algorithms and navigating in their chemical space

Cell-penetrating peptides (CPPs) are naturally able to cross the lipid bilayer membrane that protects cells. These peptides share common structural and physicochemical properties and show different pharmaceutical applications, among which drug delivery is the most important. Due to their ability to cross the membranes by pulling high-molecular-weight polar molecules, they are termed Trojan horses. In this study, we proposed a machine learning (ML)-based framework named BChemRF-CPPred (beyond chemical rules-based framework for CPP prediction) that uses an artificial neural network, a support vector machine, and a Gaussian process classifier to differentiate CPPs from non-CPPs, using structure- and sequence-based descriptors extracted from PDB and FASTA formats. The performance of our algorithm was evaluated by tenfold cross-validation and compared with those of previously reported prediction tools using an independent dataset. The BChemRF-CPPred satisfactorily identified CPP-like structures using natural and synthetic modified peptide libraries and also obtained better performance than those of previously reported ML-based algorithms, reaching the independent test accuracy of 90.66% (AUC = 0.9365) for PDB, and an accuracy of 86.5% (AUC = 0.9216) for FASTA input. Moreover, our analyses of the CPP chemical space demonstrated that these peptides break some molecular rules related to the prediction of permeability of therapeutic molecules in cell membranes. This is the first comprehensive analysis to predict synthetic and natural CPP structures and to evaluate their chemical space using an ML-based framework. Our algorithm is freely available for academic use at http://comptools.linc.ufpa.br/BChemRF-CPPred .

First total synthesis of versicotide A, B and C

The syntheses of versicotides A-C, natural products containing anthranilic acid and NMe-Ala, were achieved by solid phase peptide synthesis on 2-chlorotrityl resin followed by solution phase macrocyclization. Using an oxyma additive, the difficult coupling reactions to the deactivated aromatic amine of o-aminobenzoic acid, were performed in high yield, avoiding anthranilic rearrangements or side reactions.

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/).

PubChem and ChEMBL beyond Lipinski

Seven million of the currently 94 million entries in the PubChem database break at least one of the four Lipinski constraints for oral bioavailability, 183,185 of which are also found in the ChEMBL database. These non-Lipinski PubChem (NLP) and ChEMBL (NLC) subsets are interesting because they contain new modalities that can display biological properties not accessible to small molecule drugs. Unfortunately, the current search tools in PubChem and ChEMBL are designed for small molecules and are not well suited to explore these subsets, which therefore remain poorly appreciated. Herein we report MXFP (macromolecule extended atom-pair fingerprint), a 217-D fingerprint tailored to analyze large molecules in terms of molecular shape and pharmacophores. We implement MXFP in two web-based applications, the first one to visualize NLP and NLC interactively using Faerun (http://faerun.gdb.tools/), the second one to perform MXFP nearest neighbor searches in NLP and NLC (http://similaritysearch.gdb.tools/). We show that these tools provide a meaningful insight into the diversity of large molecules in NLP and NLC. The interactive tools presented here are publicly available at http://gdb.unibe.ch and can be used freely to explore and better understand the diversity of non-Lipinski molecules in PubChem and ChEMBL.