High-throughput functional annotation of natural products by integrated activity profiling.

Determining mechanism of action (MOA) is one of the biggest challenges in natural products discovery. Here, we report a comprehensive platform that uses Similarity Network Fusion (SNF) to improve MOA predictions by integrating data from the cytological profiling high-content imaging platform and the gene expression platform Functional Signature Ontology, and pairs these data with untargeted metabolomics analysis for de novo bioactive compound discovery. The predictive value of the integrative approach was assessed using a library of target-annotated small molecules as benchmarks. Using Kolmogorov-Smirnov (KS) tests to compare in-class to out-of-class similarity, we found that SNF retains the ability to identify significant in-class similarity across a diverse set of target classes, and could find target classes not detectable in either platform alone. This confirmed that integration of expression-based and image-based phenotypes can accurately report on MOA. Furthermore, we integrated untargeted metabolomics of complex natural product fractions with the SNF network to map biological signatures to specific metabolites. Three examples are presented where SNF coupled with metabolomics was used to directly functionally characterize natural products and accelerate identification of bioactive metabolites, including the discovery of the azoxy-containing biaryl compounds parkamycins A and B. Our results support SNF integration of multiple phenotypic screening approaches along with untargeted metabolomics as a powerful approach for advancing natural products drug discovery.

Membrane Permeability in a Large Macrocyclic Peptide Driven by a Saddle-Shaped Conformation.

The effort to modulate challenging protein targets has stimulated interest in ligands that are larger and more complex than typical small-molecule drugs. While combinatorial techniques such as mRNA display routinely produce high-affinity macrocyclic peptides against classically undruggable targets, poor membrane permeability has limited their use toward primarily extracellular targets. Understanding the passive membrane permeability of macrocyclic peptides would, in principle, improve our ability to design libraries whose leads can be more readily optimized against intracellular targets. Here, we investigate the permeabilities of over 200 macrocyclic 10-mers using the thioether cyclization motif commonly found in mRNA display macrocycle libraries. We identified the optimal lipophilicity range for achieving permeability in thioether-cyclized 10-mer cyclic peptide-peptoid hybrid scaffolds and showed that permeability could be maintained upon extensive permutation in the backbone. In one case, changing a single amino acid from d-Pro to d-NMe-Ala, representing the loss of a single methylene group in the side chain, resulted in a highly permeable scaffold in which the low-dielectric conformation shifted from the canonical cross-beta geometry of the parent compounds into a novel saddle-shaped fold in which all four backbone NH groups were sequestered from the solvent. This work provides an example by which pre-existing physicochemical knowledge of a scaffold can benefit the design of macrocyclic peptide mRNA display libraries, pointing toward an approach for biasing libraries toward permeability by design. Moreover, the compounds described herein are a further demonstration that geometrically diverse, highly permeable scaffolds exist well beyond conventional drug-like chemical space.

High-Content Image-Based Screening and Deep Learning for the Detection of Anti-Inflammatory Drug Leads.

We developed a high-content image-based screen that utilizes the pro-inflammatory stimulus lipopolysaccharide (LPS) and murine macrophages (RAW264.7) with the goal of enabling the identification of novel anti-inflammatory lead compounds. We screened 2,259 bioactive compounds with annotated mechanisms of action (MOA) to identify compounds that block the LPS-induced phenotype in macrophages. We utilized a set of seven fluorescence microscopy probes to generate images that were used to train and optimize a deep neural network classifier to distinguish between unstimulated and LPS-stimulated macrophages. The top hits from the deep learning classifier were validated using a linear classifier trained on individual cells and subsequently investigated in a multiplexed cytokine secretion assay. All 12 hits significantly modulated the expression of at least one cytokine upon LPS stimulation. Seven of these were allosteric inhibitors of the mitogen-activated protein kinase kinase (MEK1/2) and showed similar effects on cytokine expression. This deep learning morphological assay identified compounds that modulate the innate immune response to LPS and may aid in identifying new anti-inflammatory drug leads.

Evidence for aggregation-independent, PrPC-mediated Aß cellular internalization.

Evidence linking amyloid beta (Aß) cellular uptake and toxicity has burgeoned, and mechanisms underlying this association are subjects of active research. Two major, interconnected questions are whether Aß uptake is aggregation-dependent and whether it is sequence-specific. We recently reported that the neuronal uptake of Aß depends significantly on peptide chirality, suggesting that the process is predominantly receptor-mediated. Over the past decade, the cellular prion protein (PrPC) has emerged as an important mediator of Aß-induced toxicity and of neuronal Aß internalization. Here, we report that the soluble, nonfibrillizing Aß (1-30) peptide recapitulates full-length Aß stereoselective cellular uptake, allowing us to decouple aggregation from cellular, receptor-mediated internalization. Moreover, we found that Aß (1-30) uptake is also dependent on PrPC expression. NMR-based molecular-level characterization identified the docking site on PrPC that underlies the stereoselective binding of Aß (1-30). Our findings therefore identify a specific sequence within Aß that is responsible for the recognition of the peptide by PrPC, as well as PrPC-dependent cellular uptake. Further uptake stereodifferentiation in PrPC-free cells points toward additional receptor-mediated interactions as likely contributors for Aß cellular internalization. Taken together, our results highlight the potential of targeting cellular surface receptors to inhibit Aß cellular uptake as an alternative route for future therapeutic development for Alzheimer's disease.

Geometrically Diverse Lariat Peptide Scaffolds Reveal an Untapped Chemical Space of High Membrane Permeability.

Constrained, membrane-permeable peptides offer the possibility of engaging challenging intracellular targets. Structure-permeability relationships have been extensively studied in cyclic peptides whose backbones are cyclized from head to tail, like the membrane permeable and orally bioavailable natural product cyclosporine A. In contrast, the physicochemical properties of lariat peptides, which are cyclized from one of the termini onto a side chain, have received little attention. Many lariat peptide natural products exhibit interesting biological activities, and some, such as griselimycin and didemnin B, are membrane permeable and have intracellular targets. To investigate the structure-permeability relationships in the chemical space exemplified by these natural products, we generated a library of scaffolds using stable isotopes to encode stereochemistry and determined the passive membrane permeability of over 1000 novel lariat peptide scaffolds with molecular weights around 1000. Many lariats were surprisingly permeable, comparable to many known orally bioavailable drugs. Passive permeability was strongly dependent on N-methylation, stereochemistry, and ring topology. A variety of structure-permeability trends were observed including a relationship between alternating stereochemistry and high permeability, as well as a set of highly permeable consensus sequences. For the first time, robust structure-permeability relationships are established in synthetic lariat peptides exceeding 1000 compounds.

Developing Cyclic Peptomers as Broad-Spectrum Type III Secretion System Inhibitors in Gram-Negative Bacteria.

Antibiotic-resistant bacteria are an emerging global health threat. New antimicrobials are urgently needed. The injectisome type III secretion system (T3SS), required by dozens of Gram-negative bacteria for virulence but largely absent from nonpathogenic bacteria, is an attractive antimicrobial target. We previously identified synthetic cyclic peptomers, inspired by the natural product phepropeptin D, that inhibit protein secretion through the Yersinia Ysc and Pseudomonas aeruginosa Psc T3SSs but do not inhibit bacterial growth. Here, we describe the identification of an isomer, 4EpDN, that is 2-fold more potent (50% inhibitory concentration [IC50] of 4 µM) than its parental compound. Furthermore, 4EpDN inhibited the Yersinia Ysa and the Salmonella SPI-1 T3SSs, suggesting that this cyclic peptomer has broad efficacy against evolutionarily distant injectisome T3SSs. Indeed, 4EpDN strongly inhibited intracellular growth of Chlamydia trachomatis in HeLa cells, which requires the T3SS. 4EpDN did not inhibit the unrelated twin arginine translocation (Tat) system, nor did it impact T3SS gene transcription. Moreover, although the injectisome and flagellar T3SSs are evolutionarily and structurally related, the 4EpDN cyclic peptomer did not inhibit secretion of substrates through the Salmonella flagellar T3SS, indicating that cyclic peptomers broadly but specifically target the injectisome T3SS. 4EpDN reduced the number of T3SS needles detected on the surface of Yersinia pseudotuberculosis as detected by microscopy. Collectively, these data suggest that cyclic peptomers specifically inhibit the injectisome T3SS from a variety of Gram-negative bacteria, possibly by preventing complete T3SS assembly.

Cyclosporin A: Conformational Complexity and Chameleonicity.

The chameleonic behavior of cyclosporin A (CsA) was investigated through conformational ensembles employing multicanonical molecular dynamics simulations that could sample the cis and trans isomers of N-methylated amino acids; these assessments were conducted in explicit water, dimethyl sulfoxide, acetonitrile, methanol, chloroform, cyclohexane (CHX), and n-hexane (HEX) using AMBER ff03, AMBER10:EHT, AMBER12:EHT, and AMBER14:EHT force fields. The conformational details were discussed employing the free-energy landscapes (FELs) at T = 300 K; it was observed that the experimentally determined structures of CsA were only a part of the conformational space. Comparing the ROESY measurements in CHX-d12 and HEX-d14, the major conformations in those apolar solvents were essentially the same as that in CDCl3 except for the observation of some sidechain rotamers. The effects of the metal ions on the conformations, including the cis/trans isomerization, were also investigated. Based on the analysis of FELs, it was concluded that the AMBER ff03 force field best described the experimentally derived conformations, indicating that CsA intrinsically formed membrane-permeable conformations and that the metal ions might be the key to the cis/trans isomerization of N-methylated amino acids before binding a partner protein.

Connecting the conformational behavior of cyclic octadepsipeptides with their ionophoric property and membrane permeability.

Cyclic octadepsipeptides such as PF1022A and its synthetic derivative emodepside exhibit anthelmintic activity with the latter sold as a commercial drug treatment against gastrointestinal nematodes for animal health use. The structure-permeability relationship of these cyclic depsipeptides that could ultimately provide insights into the compound bioavailability is not yet well understood. The fully N-methylated amide backbone and apolar sidechain residues do not allow for the formation of intramolecular hydrogen bonds, normally observed in the membrane-permeable conformations of cyclic peptides. Hence, any understanding gained on these depsipeptides would serve as a prototype for future design strategies. In previous nuclear magnetic resonance (NMR) studies, two macrocyclic core conformers of emodepside were detected, one with all backbone amides in trans-configuration (hereon referred as the symmetric conformer) and the other with one amide in cis-configuration (hereon referred as the asymmetric conformer). In addition, these depsipeptides were also reported to be ionophores with a preference of potassium over sodium. In this study, we relate the conformational behavior of PF1022A, emodepside, and closely related analogs with their ionophoric characteristic probed using NMR and molecular dynamics (MD) simulations and finally evaluated their passive membrane permeability using PAMPA. We find that the equilibrium between the two core conformers shifts more towards the symmetric conformer upon addition of monovalent cations with selectivity for potassium over sodium. Both the NMR experiments and the theoretical Markov state models based on extensive MD simulations indicate a more rigid backbone for the asymmetric conformation, whereas the symmetric conformation shows greater flexibility. The experimental results further advocate for the symmetric conformation binding the cation. The PAMPA results suggest that the investigated depsipeptides are retained in the membrane, which may be advantageous for the likely target, a membrane-bound potassium channel.

Retrospective analysis of natural products provides insights for future discovery trends.

Understanding of the capacity of the natural world to produce secondary metabolites is important to a broad range of fields, including drug discovery, ecology, biosynthesis, and chemical biology, among others. Both the absolute number and the rate of discovery of natural products have increased significantly in recent years. However, there is a perception and concern that the fundamental novelty of these discoveries is decreasing relative to previously known natural products. This study presents a quantitative examination of the field from the perspective of both number of compounds and compound novelty using a dataset of all published microbial and marine-derived natural products. This analysis aimed to explore a number of key questions, such as how the rate of discovery of new natural products has changed over the past decades, how the average natural product structural novelty has changed as a function of time, whether exploring novel taxonomic space affords an advantage in terms of novel compound discovery, and whether it is possible to estimate how close we are to having described all of the chemical space covered by natural products. Our analyses demonstrate that most natural products being published today bear structural similarity to previously published compounds, and that the range of scaffolds readily accessible from nature is limited. However, the analysis also shows that the field continues to discover appreciable numbers of natural products with no structural precedent. Together, these results suggest that the development of innovative discovery methods will continue to yield compounds with unique structural and biological properties.

Drug-Like Properties in Macrocycles above MW 1000: Backbone Rigidity versus Side-Chain Lipophilicity.

Large macrocyclic peptides can achieve surprisingly high membrane permeability, although the properties that govern permeability in this chemical space are only beginning to come into focus. We generated two libraries of cyclic decapeptides with stable cross-ß conformations, and found that peptoid substitutions within the ß-turns of the macrocycle preserved the rigidity of the parent scaffold, whereas peptoid substitutions in the opposing ß-strands led to "chameleonic" species that were rigid in nonpolar media but highly flexible in water. Both rigid and chameleonic compounds showed high permeability over a wide lipophilicity range, with peak permeabilities differing significantly depending on scaffold rigidity. Our findings indicate that modulating lipophilicity can be used to engineer favorable ADME properties into both rigid and flexible macrocyclic peptides, and that scaffold rigidity can be used to tune optimal lipophilicity.

Conformation and Permeability: Cyclic Hexapeptide Diastereomers.

Conformational ensembles of eight cyclic hexapeptide diastereomers in explicit cyclohexane, chloroform, and water were analyzed by multicanonical molecular dynamics (McMD) simulations. Free-energy landscapes (FELs) for each compound and solvent were obtained from the molecular shapes and principal component analysis at T = 300 K; detailed analysis of the conformational ensembles and flexibility of the FELs revealed that permeable compounds have different structural profiles even for a single stereoisomeric change. The average solvent-accessible surface area (SASA) in cyclohexane showed excellent correlation with the cell permeability, whereas this correlation was weaker in chloroform. The average SASA in water correlated with the aqueous solubility. The average polar surface area did not correlate with cell permeability in these solvents. A possible strategy for designing permeable cyclic peptides from FELs obtained from McMD simulations is proposed.

CycLS: Accurate, whole-library sequencing of cyclic peptides using tandem mass spectrometry.

Cyclic peptides are of great interest as therapeutic compounds due to their potential for specificity and intracellular activity, but specific compounds can be difficult to identify from large libraries without resorting to molecular encoding techniques. Large libraries of cyclic peptides are often DNA-encoded or linearized before sequencing, but both of those deconvolution strategies constrain the chemistry, assays, and quantification methods which can be used. We developed an automated sequencing program, CycLS, to identify cyclic peptides contained within large synthetic libraries. CycLS facilitates quick and easy identification of all library-members via tandem mass spectrometry data without requiring any specific chemical moieties or modifications within the library. Validation of CycLS against a library of 400 cyclic hexapeptide peptoid hybrids (peptomers) of unique mass yielded a result of 95% accuracy when compared against a simulated library size of 234,256 compounds. CycLS was also evaluated by resynthesizing pure compounds from a separate 1800-member library of cyclic hexapeptides and hexapeptomers with high mass redundancy. Of 22 peptides resynthesized, 17 recapitulated the retention times and fragmentation patterns assigned to them from the whole-library bulk assay results. Implementing a database-matching approach, CycLS is fast and provides a robust method for sequencing cyclic peptides that is particularly applicable to the deconvolution of synthetic libraries.

Using 1H and 13C NMR chemical shifts to determine cyclic peptide conformations: a combined molecular dynamics and quantum mechanics approach.

Solving conformations of cyclic peptides can provide insight into structure-activity and structure-property relationships, which can help in the design of compounds with improved bioactivity and/or ADME characteristics. The most common approaches for determining the structures of cyclic peptides are based on NMR-derived distance restraints obtained from NOESY or ROESY cross-peak intensities, and 3J-based dihedral restraints using the Karplus relationship. Unfortunately, these observables are often too weak, sparse, or degenerate to provide unequivocal, high-confidence solution structures, prompting us to investigate an alternative approach that relies only on 1H and 13C chemical shifts as experimental observables. This method, which we call conformational analysis from NMR and density-functional prediction of low-energy ensembles (CANDLE), uses molecular dynamics (MD) simulations to generate conformer families and density functional theory (DFT) calculations to predict their 1H and 13C chemical shifts. Iterative conformer searches and DFT energy calculations on a cyclic peptide-peptoid hybrid yielded Boltzmann ensembles whose predicted chemical shifts matched the experimental values better than any single conformer. For these compounds, CANDLE outperformed the classic NOE- and 3J-coupling-based approach by disambiguating similar ß-turn types and also enabled the structural elucidation of the minor conformer. Through the use of chemical shifts, in conjunction with DFT and MD calculations, CANDLE can help illuminate conformational ensembles of cyclic peptides in solution.

Genome-Based Identification of a Plant-Associated Marine Bacterium as a Rich Natural Product Source.

The large number of sequenced bacterial genomes provides the opportunity to bioinformatically identify rich natural product sources among previously neglected microbial groups. Testing this discovery strategy, unusually high biosynthetic potential was suggested for the Oceanospirillales member Gynuella sunshinyii, a Gram-negative marine bacterium from the rhizosphere of the halophilic plant Carex scabrifolia. Its genome contains numerous unusual biosynthetic gene clusters for diverse types of metabolites. Genome-guided isolation yielded representatives of four different natural product classes, of which only alteramide A was known. Cytotoxic lacunalides were identified as products of a giant trans-acyltransferase polyketide synthase gene cluster, one of six present in this strain. Cytological profiling against HeLa cells suggested that lacunalide A disrupts CDK signaling in the cell cycle. In addition, chemical studies on model compounds were conducted, suggesting the structurally unusual ergoynes as products of a conjugated diyne-thiourea cyclization reaction.

Naphthablins B and C, Meroterpenoids Identified from the Marine Sediment-Derived Streptomyces sp. CP26-58 Using HeLa Cell-Based Cytological Profiling.

HeLa cell-based cytological profiling (CP) was applied to an extract library of marine sediment-derived actinomycetes to discover new cytotoxic secondary metabolites. Among the hit strains, Streptomyces sp. CP26-58 was selected for further investigation to identify its cytotoxic metabolites. CP revealed that the known ionophore tetronasin (1) was responsible for the cytotoxic effect found in the extract. Furthermore, three naphthoquinone meroterpenoids, naphthablin A (2) and two new derivatives designated as naphthablins B (3) and C (4), were isolated from other cytotoxic fractions. The structures of the new compounds were elucidated based on analysis of their HRESIMS and comprehensive NMR data. The absolute configurations of the new compounds were deduced by simulating ECD spectra and calculating potential energies for the model compounds using density function theory (DFT) calculations. Compound 1 showed a significant cytotoxic effect against HeLa cells with an IC50 value of 0.23 µM, and CP successfully clustered 1 with calcium ionophores.

Cell-permeable cyclic peptides from synthetic libraries inspired by natural products.

Drug design efforts are turning to a new generation of therapeutic targets, such as protein-protein interactions (PPIs), that had previously been considered "undruggable" by typical small molecules. There is an emerging view that accessing these targets will require molecules that are larger and more complex than typical small molecule drugs. Here, we present a methodology for the discovery of geometrically diverse, membrane permeable cyclic peptide scaffolds based on the synthesis and permeability screening of a combinatorial library, followed by deconvolution of membrane-permeable scaffolds to identify cyclic peptides with good to excellent passive cell permeabilities. We use a combination of experimental and computational approaches to investigate structure-permeability relationships in one of these scaffolds, and uncover structural and conformational factors that govern passive membrane diffusion in a related set of cyclic peptide diastereomers. Further, we investigate the dependency of permeability on side-chain identity of one of these scaffolds through single-point diversifications to show the adaptability of these scaffolds toward development of permeability-biased libraries suitable for bioactivity screens. Overall, our results demonstrate that many novel, cell permeable scaffolds exist beyond those found in extant natural products, and that such scaffolds can be rapidly identified using a combination of synthesis and deconvolution which can, in principle, be applied to any type of macrocyclic template.

Phenotype-Guided Natural Products Discovery Using Cytological Profiling.

Phenotype-guided natural products discovery is emerging as a useful new discovery tool that addresses challenges in early, unbiased natural product biological annotation. These high-content approaches yield screening results that report directly on the impact of test compounds on cellular processes in target organisms and can be used to predict the modes of action of bioactive constituents from primary screening data. In this study we explored the use of our recently implemented cytological profiling platform for the isolation of compounds with a specific, predefined mode of action, namely, induction of mitotic arrest. Screening of a microbially derived extract library revealed six extracts whose cytological profiles clustered closely with those of known antimitotic agents from the pure compound training set. Subsequent examination of one of these extracts revealed the presence of two separate bioactive constituents, each of which possessed a unique cytological profile. The first, diketopiperazine XR334 (3), recapitulated the observed antimitotic phenotype of the original extract, demonstrating that cytological profiling can be used for the targeted isolation of compounds with specific modes of action. The second, nocapyrone L (6), possessed a cytological profile that clustered with known calcium channel modulators, in line with previous published activities for this compound class, indicating that cytological profiling is a flexible and powerful platform for the de novo characterization of compound modes of action.

Biosynthetic products from a nearshore-derived gram-negative bacterium enable reassessment of the kailuin depsipeptides.

Sampling of California nearshore sediments resulted in the isolation of a Gram-negative bacterium, Photobacterium halotolerans, capable of producing unusual biosynthetic products. Liquid culture in artificial seawater-based media provided cyclic depsipeptides including four known compounds, kailuins B-E (2-5), and two new analogues, kailuins G and H (7 and 8). The structures of the new and known compounds were confirmed through extensive spectroscopic and Marfey's analyses. During the course of these studies, a correction was made to the previously reported double-bond geometry of kailuin D (4). Additionally, through the application of a combination of derivatization with Mosher's reagent and extensive (13)C NMR shift analysis, the previously unassigned chiral center at position C-3 of the ß-acyloxy group of all compounds was determined. To evaluate bioactivity and structure-activity relationships, the kailuin core (13) and kailuin lactam (14) were prepared by chiral synthesis using an Fmoc solid-phase peptide strategy followed by solution-phase cyclization. All isolated compounds and synthetic cores were assayed for solid tumor cell cytotoxicity and showed only minimal activity, contrary to other published reports. Additional phenotypic screenings were done on 4 and 5, with little evidence of activity.

A cell-based screen reveals that the albendazole metabolite, albendazole sulfone, targets Wolbachia.

Wolbachia endosymbionts carried by filarial nematodes give rise to the neglected diseases African river blindness and lymphatic filariasis afflicting millions worldwide. Here we identify new Wolbachia-disrupting compounds by conducting high-throughput cell-based chemical screens using a Wolbachia-infected, fluorescently labeled Drosophila cell line. This screen yielded several Wolbachia-disrupting compounds including three that resembled Albendazole, a widely used anthelmintic drug that targets nematode microtubules. Follow-up studies demonstrate that a common Albendazole metabolite, Albendazole sulfone, reduces intracellular Wolbachia titer both in Drosophila melanogaster and Brugia malayi, the nematode responsible for lymphatic filariasis. Significantly, Albendazole sulfone does not disrupt Drosophila microtubule organization, suggesting that this compound reduces titer through direct targeting of Wolbachia. Accordingly, both DNA staining and FtsZ immunofluorescence demonstrates that Albendazole sulfone treatment induces Wolbachia elongation, a phenotype indicative of binary fission defects. This suggests that the efficacy of Albendazole in treating filarial nematode-based diseases is attributable to dual targeting of nematode microtubules and their Wolbachia endosymbionts.