Classification of the adenylation and acyl-transferase activity of NRPS and PKS systems using ensembles of substrate specific hidden Markov models

Structural, biochemical and bioinformatic analyses of nonribosomal peptide synthetase adenylation domains

Covering: 1997 to July 2023The adenylation reaction has been a subject of scientific intrigue since it was first recognized as essential to many biological processes, including the homeostasis and pathogenicity of some bacteria and the activation of amino acids for protein synthesis in mammals. Several foundational studies on adenylation (A) domains have facilitated an improved understanding of their molecular structures and biochemical properties, in particular work on nonribosomal peptide synthetases (NRPSs). In NRPS pathways, A domains activate their respective acyl substrates for incorporation into a growing peptidyl chain, and many nonribosomal peptides are bioactive. From a natural product drug discovery perspective, improving existing bioinformatics platforms to predict unique NRPS products more accurately from genomic data is desirable. Here, we summarize characterization efforts of A domains primarily from NRPS pathways from July 1997 up to July 2023, covering protein structure elucidation, in vitro assay development, and in silico tools for improved predictions.

Discovery and Biosynthesis of Anthrochelin, a Growth-Promoting Metallophore of the Human Pathogen Luteibacter anthropi

Various human pathogens have emerged from environmental strains by adapting to higher growth temperatures and the ability to produce virulence factors. A remarkable example of a pathoadapted bacterium is found in the genus Luteibacter, which typically comprises harmless soil microbes, yet Luteibacter anthropi was isolated from the blood of a diseased child. Up until now, nothing has been known about the specialized metabolism of this pathogen. By comparative genome analyses we found that L. anthropi has a markedly higher biosynthetic potential than other bacteria of this genus and uniquely bears an NRPS gene locus tentatively coding for the biosynthesis of a metallophore. By metabolic profiling, stable isotope labeling, and NMR investigation of a gallium complex, we identified a new family of salicylate-derived nonribosomal peptides named anthrochelins A-D. Surprisingly, anthrochelins feature a C-terminal homocysteine tag, which might be introduced during peptide termination. Mutational analyses provided insight into the anthrochelin assembly and revealed the unexpected involvement of a cytochrome P450 monooxygenase in oxazole formation. Notably, this heterocycle plays a key role in the binding of metals, especially copper(II). Bioassays showed that anthrochelin significantly promotes the growth of L. anthropi in the presence of low and high copper concentrations, which occur during infections.

Discovery and Biosynthesis of the Cytotoxic Polyene Terpenomycin in Human Pathogenic Nocardia

Nocardia are opportunistic human pathogens that can cause a range of debilitating and difficult to treat infections of the lungs, brain, skin, and soft tissues. Despite their close relationship to the well-known secondary metabolite-producing genus, Streptomyces, comparatively few natural products are known from the Nocardia, and even less is known about their involvement in the pathogenesis. Here, we combine chemistry, genomics, and molecular microbiology to reveal the production of terpenomycin, a new cytotoxic and antifungal polyene from a human pathogenic Nocardia terpenica isolate. We unveil the polyketide synthase (PKS) responsible for terpenomycin biosynthesis and show that it combines several unusual features, including "split", skipped, and iteratively used modules, and the use of the unusual extender unit methoxymalonate as a starter unit. To link genes to molecules, we constructed a transposon mutant library in N. terpenica, identifying a terpenomycin-null mutant with an inactivated terpenomycin PKS. Our findings show that the neglected actinomycetes have an unappreciated capacity for the production of bioactive molecules with unique biosynthetic pathways waiting to be uncovered and highlights these organisms as producers of diverse natural products.

Knowledge-guided data mining on the standardized architecture of NRPS: Subtypes, novel motifs, and sequence entanglements

Non-ribosomal peptide synthetase (NRPS) is a diverse family of biosynthetic enzymes for the assembly of bioactive peptides. Despite advances in microbial sequencing, the lack of a consistent standard for annotating NRPS domains and modules has made data-driven discoveries challenging. To address this, we introduced a standardized architecture for NRPS, by using known conserved motifs to partition typical domains. This motif-and-intermotif standardization allowed for systematic evaluations of sequence properties from a large number of NRPS pathways, resulting in the most comprehensive cross-kingdom C domain subtype classifications to date, as well as the discovery and experimental validation of novel conserved motifs with functional significance. Furthermore, our coevolution analysis revealed important barriers associated with re-engineering NRPSs and uncovered the entanglement between phylogeny and substrate specificity in NRPS sequences. Our findings provide a comprehensive and statistically insightful analysis of NRPS sequences, opening avenues for future data-driven discoveries.

Bioinformatic Analysis Reveals both Oversampled and Underexplored Biosynthetic Diversity in Nonribosomal Peptides

The traditional natural product discovery approach has accessed only a fraction of the chemical diversity in nature. The use of bioinformatic tools to interpret the instructions encoded in microbial biosynthetic genes has the potential to circumvent the existing methodological bottlenecks and greatly expand the scope of discovery. Structural prediction algorithms for nonribosomal peptides (NRPs), the largest family of microbial natural products, lie at the heart of this new approach. To understand the scope and limitation of the existing prediction algorithms, we evaluated their performances on NRP synthetase biosynthetic gene clusters. Our systematic analysis shows that the NRP biosynthetic landscape is uneven. Phenylglycine and its derivatives as a group of NRP building blocks (BBs), for example, have been oversampled, reflecting an extensive historical interest in the glycopeptide antibiotics family. In contrast, the benzoyl BB, including 2,3-dihydroxybenzoate (DHB), has been the most underexplored, hinting at the possibility of a reservoir of as yet unknown DHB containing NRPs with functional roles other than a siderophore. Our results also suggest that there is still vast unexplored biosynthetic diversity in nature, and the analysis presented herein shall help guide and strategize future natural product discovery campaigns. We also discuss possible ways bioinformaticians and biochemists could work together to improve the existing prediction algorithms.

Unique Initiation and Termination Mechanisms Involved in the Biosynthesis of a Hybrid Polyketide-Nonribosomal Peptide Lyngbyapeptin B Produced by the Marine Cyanobacterium Moorena bouillonii

Lyngbyapeptin B is a hybrid polyketide-nonribosomal peptide isolated from particular marine cyanobacteria. In this report, we carried out genome sequence analysis of a producer cyanobacterium Moorena bouillonii to understand the biosynthetic mechanisms that generate the unique structural features of lyngbyapeptin B, including the (E)-3-methoxy-2-butenoyl starter unit and the C-terminal thiazole moiety. We identified a putative lyngbyapeptin B biosynthetic (lynB) gene cluster comprising nine open reading frames that include two polyketide synthases (PKSs: LynB1 and LynB2), four nonribosomal peptide synthetases (NRPSs: LynB3, LynB4, LynB5, and LynB6), a putative nonheme diiron oxygenase (LynB7), a type II thioesterase (LynB8), and a hypothetical protein (LynB9). In vitro enzymatic analysis of LynB2 with methyltransferase (MT) and acyl carrier protein (ACP) domains revealed that the LynB2 MT domain (LynB2-MT) catalyzes O-methylation of the acetoacetyl-LynB2 ACP domain (LynB2-ACP) to yield (E)-3-methoxy-2-butenoyl-LynB2-ACP. In addition, in vitro enzymatic analysis of LynB7 revealed that LynB7 catalyzes the oxidative decarboxylation of (4R)-2-methyl-2-thiazoline-4-carboxylic acid to yield 2-methylthiazole in the presence of Fe²⁺ and molecular oxygen. This result indicates that LynB7 is responsible for the last post-NRPS modification to give the C-terminal thiazole moiety in lyngbyapeptin B biosynthesis. Overall, we identified and characterized a new marine cyanobacterial hybrid PKS-NRPS biosynthetic gene cluster for lyngbyapeptin B production, revealing two unique enzymatic logics.

Endofungal bacteria boost anthelminthic host protection with the biosurfactant symbiosin

Effective protection of soil fungi from predators is crucial for their survival in the niche. Thus, fungi have developed efficient defence strategies. We discovered that soil beneficial Mortierella fungi employ a potent cytotoxin (necroxime) against fungivorous nematodes. Interestingly, this anthelminthic agent is produced by bacterial endosymbionts (Candidatus Mycoavidus necroximicus) residing within the fungus. Analysis of the symbiont's genome indicated a rich biosynthetic potential, yet nothing has been known about additional metabolites and their potential synergistic functions. Here we report that two distinct Mortierella endosymbionts produce a novel cyclic lipodepsipeptide (symbiosin), that is clearly of bacterial origin, but has striking similarities to various fungal specialized metabolites. The structure and absolute configuration of symbiosin were fully elucidated. By comparative genomics of symbiosin-positive strains and in silico analyses of the deduced non-ribosomal synthetases, we assigned the (sym) biosynthetic gene cluster and proposed an assembly line model. Bioassays revealed that symbiosin is not only an antibiotic, in particular against mycobacteria, but also exhibits marked synergistic effects with necroxime in anti-nematode tests. By functional analyses and substitution experiments we found that symbiosin is a potent biosurfactant and that this particular property confers a boost in the anthelmintic action, similar to formulations of therapeutics in human medicine. Our findings illustrate that "combination therapies" against parasites already exist in ecological contexts, which may inspire the development of biocontrol agents and therapeutics.

Recent advances in the structural analysis of adenylation domains in natural product biosynthesis

Adenylation (A) domains catalyze the biosynthetic incorporation of acyl building blocks into nonribosomal peptides and related natural products by selectively transferring acyl substrates onto cognate carrier proteins (CP). The use of noncanonical acyl units, such as nonproteinogenic amino acids and keto acids, by A domains expands the structural diversity of natural products. Furthermore, interrupted A domains, which have embedded auxiliary domains, are able to modify the incorporated acyl units. Structural information on A domains is important for rational protein engineering to generate unnatural compounds. In this review, we summarize recent advances in the structural analysis of A domains. First, we discuss the mechanisms by which A domains recognize noncanonical acyl units. We then focus on the interactions of A domains with CP domains and embedded auxiliary domains.

Evolutionary Genome Mining for the Discovery and Engineering of Natural Product Biosynthesis

Genome mining has become an invaluable tool in natural products research to quickly identify and characterize the biosynthetic pathways that assemble secondary or specialized metabolites. Recently, evolutionary principles have been incorporated into genome mining strategies in an effort to better assess and prioritize novelty and understand their chemical diversification for engineering purposes. Here, we provide an introduction to the principles underlying evolutionary genome mining, including bioinformatic strategies and natural product biosynthetic databases. We introduce workflows for traditional genome mining, focusing on the popular pipeline antiSMASH, and methods to predict enzyme substrate specificity from genomic information. We then provide an in-depth discussion of evolutionary genome mining workflows, including EvoMining, CORASON, ARTS, and others, as adopted by our group for the discovery and prioritization of natural products biosynthetic gene clusters and their products.

Genetic Relationships in the Toxin-Producing Fungal Endophyte, Alternaria oxytropis Using Polyketide Synthase and Non-Ribosomal Peptide Synthase Genes

The legume Oxytropis sericea hosts a fungal endophyte, Alternaria oxytropis, which produces secondary metabolites (SM), including the toxin swainsonine. Polyketide synthase (PKS) and non-ribosomal peptide synthase (NRPS) enzymes are associated with biosynthesis of fungal SM. To better understand the origins of the SM, an unannotated genome of A. oxytropis was assessed for protein sequences similar to known PKS and NRPS enzymes of fungi. Contigs exhibiting identity with known genes were analyzed at nucleotide and protein levels using available databases. Software were used to identify PKS and NRPS domains and predict identity and function. Confirmation of sequence for selected gene sequences was accomplished using PCR. Thirteen PKS, 5 NRPS, and 4 PKS-NRPS hybrids were identified and characterized with functions including swainsonine and melanin biosynthesis. Phylogenetic relationships among closest amino acid matches with Alternaria spp. were identified for seven highly conserved PKS and NRPS, including melanin synthesis. Three PKS and NRPS were most closely related to other fungi within the Pleosporaceae family, while five PKS and PKS-NRPS were closely related to fungi in the Pleosporales order. However, seven PKS and PKS-NRPS showed no identity with fungi in the Pleosporales or the class Dothideomycetes, suggesting a different evolutionary origin for those genes.

Computationally-guided exchange of substrate selectivity motifs in a modular polyketide synthase acyltransferase

Polyketides, one of the largest classes of natural products, are often clinically relevant. The ability to engineer polyketide biosynthesis to produce analogs is critically important. Acyltransferases (ATs) of modular polyketide synthases (PKSs) catalyze the installation of malonyl-CoA extenders into polyketide scaffolds. ATs have been targeted extensively to site-selectively introduce various extenders into polyketides. Yet, a complete inventory of AT residues responsible for substrate selection has not been established, limiting the scope of AT engineering. Here, molecular dynamics simulations are used to prioritize ~50 mutations within the active site of EryAT6 from erythromycin biosynthesis, leading to identification of two previously unexplored structural motifs. Exchanging both motifs with those from ATs with alternative extender specificities provides chimeric PKS modules with expanded and inverted substrate specificity. Our enhanced understanding of AT substrate selectivity and application of this motif-swapping strategy are expected to advance our ability to engineer PKSs towards designer polyketides.

Diversity of Polyketide Synthases and Nonribosomal Peptide Synthetases Revealed Through Metagenomic Analysis of a Deep Oligotrophic Cave

Caves are considered to be extreme and challenging environments. It is believed that the ability of microorganisms to produce secondary metabolites enhances their survivability and adaptiveness in the energy-starved cave environment. Unfortunately, information on the genetic potential for the production of secondary metabolites, such as polyketides and nonribosomal peptides, is limited. In the present study, we aimed to identify and characterize genes responsible for the production of secondary metabolites in the microbial community of one of the deepest caves in the world, Krubera-Voronja Cave (43.4184 N 40.3083 E, Western Caucasus). The analysed sample materials included sediments, drinkable water from underground camps, soil and clay from the cave walls, speleothems and coloured spots from the cave walls. The type II polyketide synthases (PKSs) ketosynthases α and β and the adenylation domains of nonribosomal peptide synthetases (NRPSs) were investigated using a metagenomic approach. Taxonomic diversity analysis showed that most PKS sequences could be attributed to Actinobacteria followed by unclassified bacteria and Acidobacteria, while the NRPS sequences were more taxonomically diverse and could be assigned to Proteobacteria, Actinobacteria, Cyanobacteria, Firmicutes, Chloroflexi, etc. Only three putative metabolites could be predicted: an angucycline group polyketide, a massetolide A-like cyclic lipopeptide and a surfactin-like lipopeptide. The absolute majority of PKS and NRPS sequences showed low similarity with the sequences of the reference biosynthetic pathways, suggesting that these sequences could be involved in the production of novel secondary metabolites.

SeMPI 2.0-A Web Server for PKS and NRPS Predictions Combined with Metabolite Screening in Natural Product Databases

Microorganisms produce secondary metabolites with a remarkable range of bioactive properties. The constantly increasing amount of published genomic data provides the opportunity for efficient identification of biosynthetic gene clusters by genome mining. On the other hand, for many natural products with resolved structures, the encoding biosynthetic gene clusters have not been identified yet. Of those secondary metabolites, the scaffolds of nonribosomal peptides and polyketides (type I modular) can be predicted due to their building block-like assembly. SeMPI v2 provides a comprehensive prediction pipeline, which includes the screening of the scaffold in publicly available natural compound databases. The screening algorithm was designed to detect homologous structures even for partial, incomplete clusters. The pipeline allows linking of gene clusters to known natural products and therefore also provides a metric to estimate the novelty of the cluster if a matching scaffold cannot be found. Whereas currently available tools attempt to provide comprehensive information about a wide range of gene clusters, SeMPI v2 aims to focus on precise predictions. Therefore, the cluster detection algorithm, including building block generation and domain substrate prediction, was thoroughly refined and benchmarked, to provide high-quality scaffold predictions. In a benchmark based on 559 gene clusters, SeMPI v2 achieved comparable or better results than antiSMASH v5. Additionally, the SeMPI v2 web server provides features that can help to further investigate a submitted gene cluster, such as the incorporation of a genome browser, and the possibility to modify a predicted scaffold in a workbench before the database screening.

Genome mining as a biotechnological tool for the discovery of novel marine natural products

Compared to terrestrial environments, the oceans harbor a variety of environments, creating higher biodiversity, which gives marine natural products a high occurrence of significant biology and novel chemistry. However, traditional bioassay-guided isolation and purification strategies are severely limiting the discovery of additional novel natural products from the ocean. With an increasing number of marine microorganisms being sequenced, genome mining is gradually becoming a powerful tool to retrieve novel marine natural products. In this review, we have summarized genome mining approaches used to analyze key enzymes of biosynthetic pathways and predict the chemical structure of new gene clusters by introducing successful stories that used genome mining strategy to identify new marine-derived compounds. Furthermore, we also put forward challenges for genome mining techniques and their proposed solutions. The detailed analysis of the genome mining strategy will help researchers to understand this novel technique and its application. With the development of a genome sequence, genome mining strategies will be applied more widely, which will drive rapid development in the field of marine natural product development.

Enhancement of metabolic flux toward ε-poly-l-lysine biosynthesis by targeted inactivation of concomitant polyene macrolide biosynthesis in Streptomyces albulus

ε-Poly-l-lysine (ε-PL) produced as a secondary metabolite of Streptomyces albulus has long been used as a natural food preservative in a number of countries, including Japan, the United States, South Korea, and China. To date, numerous studies employing classical biotechnological approaches have been carried out to improve its productivity. Here we report a modern and rational genetic approach to enhancing metabolic flux toward ε-PL biosynthesis. Based on in silico genome analyses, we revealed that S. albulus NBRC14147 produces five antifungal polyene antibiotics-tetramycin A and B, tetrin A and B, and a trace amount of nystatin A1-concomitantly with antimicrobial ε-PL. Targeted inactivation of the biosynthetic gene cluster for tetramycins and tetrins in a nystatin A1 production-deficient mutant completely abolished the production of polyene macrolides, which in turn led to an approximately 20% improvement in ε-PL production that closely correlated with the polyene defects. The biosynthetic flux for ε-PL was thus successfully enhanced by inactivation of the concomitant secondary metabolite biosynthetic pathways. Since this elimination of concomitantly produced metabolites also allows for simpler purification after fermentation production of ε-PL, the rational strain engineering strategy we show here will improve its industrial production.

The Genome of the Plant-Associated Lactic Acid Bacterium Lactococcus lactis KF147 Harbors a Hybrid NRPS-PKS System Conserved in Strains of the Dental Cariogenic Streptococcus mutans

Lactococcus lactis subsp. lactis KF147 as a non-dairy strain from lactic acid bacteria (LAB) can inhabit plant tissues. It can grow on complex carbohydrates derived from plant cell walls. Its genome size is one of the largest among the sequenced lactococcal strains, possessing many genes that do not have homologues in the published genome sequences of dairy-associated L. lactis strains. In silico analysis has identified a gene cluster encoding a hybrid NRPS-PKS system (composed of non-ribosomal peptide synthetases and polyketide synthases) in the L. lactis KF147 genome, as first example of a LAB possessing such hybrid mega-enzymes. Hybrid systems produce hybrid NRP-PK secondary metabolites (natural products) in a wide variety of bacteria, fungi, and plants. In the hybrid NRPS-PKS system of L. lactis KF147, a total of 21 NRPS and 8 PKS domains were identified that are arranged into 6 NRPS modules, 3 PKS modules, and two single functional domains (trans-acyl-transferase "transAT" and thioesterase). We found homologous hybrid systems having similar gene, module, and domain organization in six other L. lactis strains and 25 strains of the dental cariogenic Streptococcus mutans. This study mainly aimed to predict the structure and function of the hybrid NRP-PK product of L. lactis KF147 using comparative genomics techniques, and included a detailed analysis of the regulatory system. Various bioinformatical approaches were used to predict the substrate specificity of the six A domains and the iterative transAT domain. Functional conservation of the A domains within different-niche-associated strains supported the prediction of the primary core structure of the putative hybrid natural product to be Leu-DLeu-Asp-DAsn-Gly-MC-MC-MC-DAsp (MC = Malonyl-CoA). Oxidative stress resistance and biofilm formation are the most probable functions of this hybrid system. The need for such a system in two different niches is argued, as an adaptation of L. lactis and S. mutans to adhere to plant tissues and human teeth, respectively, in an oxidative environment.

Metagenome Driven Discovery of Nonribosomal Peptides

Declining rates of novel natural product discovery and exponential rates of rediscovery heralded the end of the 1940s to 1960s "golden era" of antibiotic discovery. Fifty years later, the implementation of molecular screening methodologies revealed that standard culture-based screening approaches had failed to capture the vast majority of environmental bacteria and that even for the cultivable isolates only a small fraction of the biosynthetic potential had been tapped. A diversity of metagenomic screening and synthetic biology approaches have been developed to address these issues. The nonribosomal peptides have received particular focus, owing to their high levels of bioactivity and the predictability of the biosynthetic logic of the genetically encoded assembly lines that produce them. By uniting advances in next-generation sequencing and bioinformatic analysis with a diversity of traditional disciplines, several pioneering teams have proven that this previously inaccessible resource is no longer out of reach.

Directed Evolution Reveals the Functional Sequence Space of an Adenylation Domain Specificity Code

Nonribosomal peptides are important natural products biosynthesized by nonribosomal peptide synthetases (NRPSs). Adenylation (A) domains of NRPSs are highly specific for the substrate they recognize. This recognition is determined by 10 residues in the substrate-binding pocket, termed the specificity code. This finding led to the proposal that nonribosomal peptides could be altered by specificity code swapping. Unfortunately, this approach has proven, with few exceptions, to be unproductive; changing the specificity code typically results in broadened specificity or poor function. To enhance our understanding of A domain substrate selectivity, we carried out a detailed analysis of the specificity code from the A domain of EntF, an NRPS involved in enterobactin biosynthesis in Escherichia coli. Using directed evolution and a genetic selection, we determined which sites in the code have strict residue requirements and which are tolerant of variation. We showed that the EntF A domain, and other l-Ser-specific A domains, have a functional sequence space for l-Ser recognition, rather than a single code. This functional space is more expansive than the aggregate of all characterized l-Ser-specific A domains: we identified 152 new l-Ser specificity codes. Together, our data provide essential insights into how to overcome the barriers that prevent rational changes to A domain specificity.

Nonribosomal peptides in fungal cell factories: from genome mining to optimized heterologous production

Fungi are notoriously prolific producers of secondary metabolites including nonribosomal peptides (NRPs). The structural complexity of NRPs grants them interesting activities such as antibiotic, anti-cancer, and anti-inflammatory properties. The discovery of these compounds with attractive activities can be achieved by using two approaches: either by screening samples originating from various environments for their biological activities, or by identifying the related clusters in genomic sequences thanks to bioinformatics tools. This genome mining approach has grown tremendously due to recent advances in genome sequencing, which have provided an incredible amount of genomic data from hundreds of microbial species. Regarding fungal organisms, the genomic data have revealed the presence of an unexpected number of putative NRP-related gene clusters. This highlights fungi as a goldmine for the discovery of putative novel bioactive compounds. Recent development of NRP dedicated bioinformatics tools have increased the capacity to identify these gene clusters and to deduce NRPs structures, speeding-up the screening process for novel metabolites discovery. Unfortunately, the newly identified compound is frequently not or poorly produced by native producers due to a lack of expression of the related genes cluster. A frequently employed strategy to increase production rates consists in transferring the related biosynthetic pathway in heterologous hosts. This review aims to provide a comprehensive overview about the topic of NRPs discovery, from gene cluster identification by genome mining to the heterologous production in fungal hosts. The main computational tools and methods for genome mining are herein presented with an emphasis on the particularities of the fungal systems. The different steps of the reconstitution of NRP biosynthetic pathway in heterologous fungal cell factories will be discussed, as well as the key factors to consider for maximizing productivity. Several examples will be developed to illustrate the potential of heterologous production to both discover uncharacterized novel compounds predicted in silico by genome mining, and to enhance the productivity of interesting bio-active natural products.

The many faces and important roles of protein-protein interactions during non-ribosomal peptide synthesis

Covering: up to July 2018 Non-ribosomal peptide synthetase (NRPS) machineries are complex, multi-domain proteins that are responsible for the biosynthesis of many important, peptide-derived compounds. By decoupling peptide synthesis from the ribosome, NRPS assembly lines are able to access a significant pool of amino acid monomers for peptide synthesis. This is combined with a modular protein architecture that allows for great variation in stereochemistry, peptide length, cyclisation state and further modifications. The architecture of NRPS assembly lines relies upon a repetitive set of catalytic domains, which are organised into modules responsible for amino acid incorporation. Central to NRPS-mediated biosynthesis is the carrier protein (CP) domain, to which all intermediates following initial monomer activation are bound during peptide synthesis up until the final handover to the thioesterase domain that cleaves the mature peptide from the NRPS. This mechanism makes understanding the protein-protein interactions that occur between different NRPS domains during peptide biosynthesis of crucial importance to understanding overall NRPS function. This endeavour is also highly challenging due to the inherent flexibility and dynamics of NRPS systems. In this review, we present the current state of understanding of the protein-protein interactions that govern NRPS-mediated biosynthesis, with a focus on insights gained from structural studies relating to CP domain interactions within these impressive peptide assembly lines.

Functional and structural characterization of IdnL7, an adenylation enzyme involved in incednine biosynthesis

Adenylation enzymes play an important role in the selective incorporation of the cognate carboxylate substrates in natural product biosynthesis. Here, the biochemical and structural characterization of the adenylation enzyme IdnL7, which is involved in the biosynthesis of the macrolactam polyketide antibiotic incednine, is reported. Biochemical analysis showed that IdnL7 selects and activates several small amino acids. The structure of IdnL7 in complex with an L-alanyl-adenylate intermediate mimic, 5'-O-[N-(L-alanyl)sulfamoyl]adenosine, was determined at 2.1 Å resolution. The structure of IdnL7 explains the broad substrate specificity of IdnL7 towards small L-amino acids.

Structural basis of the nonribosomal codes for nonproteinogenic amino acid selective adenylation enzymes in the biosynthesis of natural products

Nonproteinogenic amino acids are the unique building blocks of nonribosomal peptides (NRPs) and hybrid nonribosomal peptide–polyketides (NRP–PKs) and contribute to their diversity of chemical structures and biological activities. In the biosynthesis of NRPs and NRP–PKs, adenylation enzymes select and activate an amino acid substrate as an aminoacyl adenylate, which reacts with the thiol of the holo form of the carrier protein to afford an aminoacyl thioester as the electrophile for the condensation reaction. Therefore, the substrate specificity of adenylation enzymes is a key determinant of the structure of NRPs and NRP–PKs. Here, we focus on nonproteinogenic amino acid selective adenylation enzymes, because understanding their unique selection mechanisms will lead to accurate functional predictions and protein engineering toward the rational biosynthesis of designed molecules containing amino acids. Based on recent progress in the structural analysis of adenylation enzymes, we discuss the nonribosomal codes of nonproteinogenic amino acid selective adenylation enzymes.

Diversified secondary metabolite biosynthesis gene repertoire revealed in symbiotic dinoflagellates

Symbiodiniaceae dinoflagellates possess smaller nuclear genomes than other dinoflagellates and produce structurally specialized, biologically active, secondary metabolites. Till date, little is known about the evolution of secondary metabolism in dinoflagellates as comparative genomic approaches have been hampered by their large genome sizes. Here, we overcome this challenge by combining genomic and metabolomics approaches to investigate how chemical diversity arises in three decoded Symbiodiniaceae genomes (clades A3, B1 and C). Our analyses identify extensive diversification of polyketide synthase and non-ribosomal peptide synthetase genes from two newly decoded genomes of Symbiodinium tridacnidorum (A3) and Cladocopium sp. (C). Phylogenetic analyses indicate that almost all the gene families are derived from lineage-specific gene duplications in all three clades, suggesting divergence for environmental adaptation. Few metabolic pathways are conserved among the three clades and we detect metabolic similarity only in the recently diverged clades, B1 and C. We establish that secondary metabolism protein architecture guides substrate specificity and that gene duplication and domain shuffling have resulted in diversification of secondary metabolism genes.

Emerging evolutionary paradigms in antibiotic discovery

Antibiotics revolutionized medicine and remain its cornerstone. Despite their global importance and the continuous threat of resistant pathogens, few antibiotics have been discovered in recent years. Natural products, especially the secondary metabolites of Actinobacteria, have been the traditional discovery source of antibiotics. In nature, the chemistry of antibiotic natural products is shaped by the unique evolution and ecology of their producing organisms, yet these influences remain largely unknown. Here, we highlight the ecology of antibiotics employed by microbes in defensive symbioses and review the evolutionary processes underlying the chemical diversity and activity of microbe-derived antibiotics, including the dynamics of vertical and lateral transmission of biosynthetic pathways and the evolution of efficacy, targeting specificity, and toxicity. We argue that a deeper understanding of the ecology and evolution of microbial interactions and the metabolites that mediate them will allow for an alternative, rational approach to discover new antibiotics.

Acyltransferases as Tools for Polyketide Synthase Engineering

Polyketides belong to the most valuable natural products, including diverse bioactive compounds, such as antibiotics, anticancer drugs, antifungal agents, immunosuppressants and others. Their structures are assembled by polyketide synthases (PKSs). Modular PKSs are composed of modules, which involve sets of domains catalysing the stepwise polyketide biosynthesis. The acyltransferase (AT) domains and their “partners”, the acyl carrier proteins (ACPs), thereby play an essential role. The AT loads the building blocks onto the “substrate acceptor”, the ACP. Thus, the AT dictates which building blocks are incorporated into the polyketide structure. The precursor- and occasionally the ACP-specificity of the ATs differ across the polyketide pathways and therefore, the ATs contribute to the structural diversity within this group of complex natural products. Those features make the AT enzymes one of the most promising tools for manipulation of polyketide assembly lines and generation of new polyketide compounds. However, the AT-based PKS engineering is still not straightforward and thus, rational design of functional PKSs requires detailed understanding of the complex machineries. This review summarizes the attempts of PKS engineering by exploiting the AT attributes for the modification of polyketide structures. The article includes 253 references and covers the most relevant literature published until May 2018.

Biosynthesis of the 15-Membered Ring Depsipeptide Neoantimycin

Antimycins are a family of natural products possessing outstanding biological activities and unique structures, which have intrigued chemists for over a half century. Of particular interest are the ring-expanded antimycins that show promising anticancer potential and whose biosynthesis remains uncharacterized. Specifically, neoantimycin and its analogs have been shown to be effective regulators of the oncogenic proteins GRP78/BiP and K-Ras. The neoantimycin structural skeleton is built on a 15-membered tetralactone ring containing one methyl, one hydroxy, one benzyl, and three alkyl moieties, as well as an amide linkage to a conserved 3-formamidosalicylic acid moiety. Although the biosynthetic gene cluster for neoantimycins was recently identified, the enzymatic logic that governs the synthesis of neoantimycins has not yet been revealed. In this work, the neoantimycin gene cluster is identified, and an updated sequence and annotation is provided delineating a nonribosomal peptide synthetase/polyketide synthase (NRPS/PKS) hybrid scaffold. Using cosmid expression and CRISPR/Cas-based genome editing, several heterologous expression strains for neoantimycin production are constructed in two separate Streptomyces species. A combination of in vivo and in vitro analysis is further used to completely characterize the biosynthesis of neoantimycins including the megasynthases and trans-acting domains. This work establishes a set of highly tractable hosts for producing and engineering neoantimycins and their C11 oxidized analogs, paving the way for neoantimycin-based drug discovery and development.

Identification of cyclosporin C from Amphichorda felina using a Cryptococcus neoformans differential temperature sensitivity assay

We used a temperature differential assay with the opportunistic fungal pathogen Cryptococcus neoformans as a simple screening platform to detect small molecules with antifungal activity in natural product extracts. By screening of a collection extracts from two different strains of the coprophilous fungus, Amphichorda felina, we detected strong, temperature-dependent antifungal activity using a two-plate agar zone of inhibition assay at 25 and 37 °C. Bioassay-guided fractionation of the crude extract followed by liquid chromatography-mass spectrometry (LC-MS) and nuclear magnetic resonance spectroscopy (NMR) identified cyclosporin C (CsC) as the main component of the crude extract responsible for growth inhibition of C. neoformans at 37 °C. The presence of CsC was confirmed by comparison with a commercial standard. We sequenced the genome of A. felina to identify and annotate the CsC biosynthetic gene cluster. The only previously characterized gene cluster for the biosynthesis of similar compounds is that of the related immunosuppressant drug cyclosporine A (CsA). The CsA and CsC gene clusters share a high degree of synteny and sequence similarity. Amino acid changes in the adenylation domain of the CsC nonribosomal peptide synthase's sixth module may be responsible for the substitution of L-threonine compared to L-α-aminobutyric acid in the CsA peptide core. This screening strategy promises to yield additional antifungal natural products with a focused spectrum of antimicrobial activity.

SANDPUMA: ensemble predictions of nonribosomal peptide chemistry reveal biosynthetic diversity across Actinobacteria

SUMMARY

Nonribosomally synthesized peptides (NRPs) are natural products with widespread applications in medicine and biotechnology. Many algorithms have been developed to predict the substrate specificities of nonribosomal peptide synthetase adenylation (A) domains from DNA sequences, which enables prioritization and dereplication, and integration with other data types in discovery efforts. However, insufficient training data and a lack of clarity regarding prediction quality have impeded optimal use. Here, we introduce prediCAT, a new phylogenetics-inspired algorithm, which quantitatively estimates the degree of predictability of each A-domain. We then systematically benchmarked all algorithms on a newly gathered, independent test set of 434 A-domain sequences, showing that active-site-motif-based algorithms outperform whole-domain-based methods. Subsequently, we developed SANDPUMA, a powerful ensemble algorithm, based on newly trained versions of all high-performing algorithms, which significantly outperforms individual methods. Finally, we deployed SANDPUMA in a systematic investigation of 7635 Actinobacteria genomes, suggesting that NRP chemical diversity is much higher than previously estimated. SANDPUMA has been integrated into the widely used antiSMASH biosynthetic gene cluster analysis pipeline and is also available as an open-source, standalone tool.

AVAILABILITY AND IMPLEMENTATION

SANDPUMA is freely available at https://bitbucket.org/chevrm/sandpuma and as a docker image at https://hub.docker.com/r/chevrm/sandpuma/ under the GNU Public License 3 (GPL3).

CONTACT

chevrette@wisc.edu or marnix.medema@wur.nl.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Chrysogine Biosynthesis Is Mediated by a Two-Module Nonribosomal Peptide Synthetase

Production of chrysogine has been reported from several fungal genera including Penicillium, Aspergillus, and Fusarium. Anthranilic acid and pyruvic acid, which are expected precursors of chrysogine, enhance production of this compound. A possible route for the biosynthesis using these substrates is via a nonribosomal peptide synthetase (NRPS). Through comparative analysis of the NRPSs from genome-sequenced producers of chrysogine we identified a candidate NRPS cluster comprising five additional genes named chry2-6. Deletion of the two-module NRPS (NRPS14 = chry1) abolished chrysogine production in Fusarium graminearum, indicating that the gene cluster is responsible for chrysogine biosynthesis. Overexpression of NRPS14 enhanced chrysogine production, suggesting that the NRPS is the bottleneck in the biosynthetic pathway.

antiSMASH 4.0-improvements in chemistry prediction and gene cluster boundary identification

Many antibiotics, chemotherapeutics, crop protection agents and food preservatives originate from molecules produced by bacteria, fungi or plants. In recent years, genome mining methodologies have been widely adopted to identify and characterize the biosynthetic gene clusters encoding the production of such compounds. Since 2011, the 'antibiotics and secondary metabolite analysis shell-antiSMASH' has assisted researchers in efficiently performing this, both as a web server and a standalone tool. Here, we present the thoroughly updated antiSMASH version 4, which adds several novel features, including prediction of gene cluster boundaries using the ClusterFinder method or the newly integrated CASSIS algorithm, improved substrate specificity prediction for non-ribosomal peptide synthetase adenylation domains based on the new SANDPUMA algorithm, improved predictions for terpene and ribosomally synthesized and post-translationally modified peptides cluster products, reporting of sequence similarity to proteins encoded in experimentally characterized gene clusters on a per-protein basis and a domain-level alignment tool for comparative analysis of trans-AT polyketide synthase assembly line architectures. Additionally, several usability features have been updated and improved. Together, these improvements make antiSMASH up-to-date with the latest developments in natural product research and will further facilitate computational genome mining for the discovery of novel bioactive molecules.

Diversity Profile and Dynamics of Peptaibols Produced by Green Mould Trichoderma Species in Interactions with Their Hosts Agaricus bisporus and Pleurotus ostreatus

Certain Trichoderma species are causing serious losses in mushroom production worldwide. Trichoderma aggressivum and Trichoderma pleuroti are among the major causal agents of the green mould diseases affecting Agaricus bisporus and Pleurotus ostreatus, respectively. The genus Trichoderma is well-known for the production of bioactive secondary metabolites, including peptaibols, which are short, linear peptides containing unusual amino acid residues and being synthesised via non-ribosomal peptide synthetases (NRPSs). The aim of this study was to get more insight into the peptaibol production of T. aggressivum and T. pleuroti. HPLC/MS-based methods revealed the production of peptaibols closely related to hypomurocins B by T. aggressivum, while tripleurins representing a new group of 18-residue peptaibols were identified in T. pleuroti. Putative NRPS genes enabling the biosynthesis of the detected peptaibols could be found in the genomes of both Trichoderma species. In vitro experiments revealed that peptaibols are potential growth inhibitors of mushroom mycelia, and that the host mushrooms may have an influence on the peptaibol profiles of green mould agents.

Non-ribosomal peptide synthetases: Identifying the cryptic gene clusters and decoding the natural product

Non-ribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs) present in bacteria and fungi are the major multi-modular enzyme complexes which synthesize secondary metabolites like the pharmacologically important antibiotics and siderophores. Each of the multiple modules of an NRPS activates a different amino or aryl acid, followed by their condensation to synthesize a linear or cyclic natural product. The studies on NRPS domains, the knowledge of their gene cluster architecture and tailoring enzymes have helped in the in silico genetic screening of the ever-expanding sequenced microbial genomic data for the identification of novel NRPS/PKS clusters and thus deciphering novel non-ribosomal peptides (NRPs). Adenylation domain is an integral part of the NRPSs and is the substrate selecting unit for the final assembled NRP. In some cases, it also requires a small protein, the MbtH homolog, for its optimum activity. The presence of putative adenylation domain and MbtH homologs in a sequenced genome can help identify the novel secondary metabolite producers. The role of the adenylation domain in the NRPS gene clusters and its characterization as a tool for the discovery of novel cryptic NRPS gene clusters are discussed.

Polyketide and nonribosomal peptide retro-biosynthesis and global gene cluster matching

Polyketides (PKs) and nonribosomal peptides (NRPs) are profoundly important natural products, forming the foundations of many therapeutic regimes. Decades of research have revealed over 11,000 PK and NRP structures, and genome sequencing is uncovering new PK and NRP gene clusters at an unprecedented rate. However, only ∼10% of PK and NRPs are currently associated with gene clusters, and it is unclear how many of these orphan gene clusters encode previously isolated molecules. Therefore, to efficiently guide the discovery of new molecules, we must first systematically de-orphan emergent gene clusters from genomes. Here we provide to our knowledge the first comprehensive retro-biosynthetic program, generalized retro-biosynthetic assembly prediction engine (GRAPE), for PK and NRP families and introduce a computational pipeline, global alignment for natural products cheminformatics (GARLIC), to uncover how observed biosynthetic gene clusters relate to known molecules, leading to the identification of gene clusters that encode new molecules.

Standardization for natural product synthetic biology

Standardization is one of the foundational features of modern-day engineering, and the use of standardized parts and processes is a key element that distinguishes bona fide synthetic biology from traditional genetic engineering. Here, we discuss the role of standardization in natural product synthetic biology, focusing on standardization of data on biosynthetic pathways and gene clusters, as well as the role of standardization in the process of biosynthetic gene cluster engineering.

Characterization of Three Tailoring Enzymes in Dutomycin Biosynthesis and Generation of a Potent Antibacterial Analogue

The anthracycline natural product dutomycin and its precursor POK-MD1 were isolated from Streptomyces minoensis NRRL B-5482. The dutomycin biosynthetic gene cluster was identified by genome sequencing and disruption of the ketosynthase gene. Two polyketide synthase (PKS) systems are present in the gene cluster, including a type II PKS and a rare highly reducing iterative type I PKS. The type I PKS DutG repeatedly uses its active sites to create a nine-carbon triketide chain that is subsequently transferred to the α-l-axenose moiety of POK-MD1 at 4″-OH to yield dutomycin. Using a heterologous recombination approach, we disrupted a putative methyltransferase gene (dutMT1) and two glycosyltransferase genes (dutGT1 and dutGT2). Analysis of the metabolites of these mutants revealed the functions of these genes and yielded three dutomycin analogues SW140, SW91, and SW75. The major product SW91 in Streptomyces minoensis NRRL B-5482-ΔDutMT1 was identified as 12-desmethyl-dutomycin, suggesting that DutMT1 is the dedicated 12-methyltransferase. This was confirmed by the in vitro enzymatic assay. DutGT1 and DutGT2 were found to be responsible for the introduction of β-d-amicetose and α-l-axenose, respectively. Dutomycin and SW91 showed strong antibacterial activity against Staphylococcus aureus and methicillin-resistant S. aureus, whereas POK-MD1 and SW75 had no obvious inhibition, which revealed the essential role of the C-4″ triketide chain in antibacterial activity. The minimal inhibitory concentration of SW91 against the two strains was 0.125 μg mL(-1), lower than that of dutomycin (0.25 μg mL(-1)), indicating that the antibacterial activity of dutomycin can be improved through biosynthetic structural modification.

Spotlight on Antimicrobial Metabolites from the Marine Bacteria Pseudoalteromonas: Chemodiversity and Ecological Significance

This review is dedicated to the antimicrobial metabolite-producing Pseudoalteromonas strains. The genus Pseudoalteromonas hosts 41 species, among which 16 are antimicrobial metabolite producers. To date, a total of 69 antimicrobial compounds belonging to 18 different families have been documented. They are classified into alkaloids, polyketides, and peptides. Finally as Pseudoalteromonas strains are frequently associated with macroorganisms, we can discuss the ecological significance of antimicrobial Pseudoalteromonas as part of the resident microbiota.

Discovery Strategies of Bioactive Compounds Synthesized by Nonribosomal Peptide Synthetases and Type-I Polyketide Synthases Derived from Marine Microbiomes

Considering that 70% of our planet's surface is covered by oceans, it is likely that undiscovered biodiversity is still enormous. A large portion of marine biodiversity consists of microbiomes. They are very attractive targets of bioprospecting because they are able to produce a vast repertoire of secondary metabolites in order to adapt in diverse environments. In many cases secondary metabolites of pharmaceutical and biotechnological interest such as nonribosomal peptides (NRPs) and polyketides (PKs) are synthesized by multimodular enzymes named nonribosomal peptide synthetases (NRPSes) and type-I polyketide synthases (PKSes-I), respectively. Novel findings regarding the mechanisms underlying NRPS and PKS evolution demonstrate how microorganisms could leverage their metabolic potential. Moreover, these findings could facilitate synthetic biology approaches leading to novel bioactive compounds. Ongoing advances in bioinformatics and next-generation sequencing (NGS) technologies are driving the discovery of NRPs and PKs derived from marine microbiomes mainly through two strategies: genome-mining and metagenomics. Microbial genomes are now sequenced at an unprecedented rate and this vast quantity of biological information can be analyzed through genome mining in order to identify gene clusters encoding NRPSes and PKSes of interest. On the other hand, metagenomics is a fast-growing research field which directly studies microbial genomes and their products present in marine environments using culture-independent approaches. The aim of this review is to examine recent developments regarding discovery strategies of bioactive compounds synthesized by NRPS and type-I PKS derived from marine microbiomes and to highlight the vast diversity of NRPSes and PKSes present in marine environments by giving examples of recently discovered bioactive compounds.

Identification of the non-ribosomal peptide synthetase responsible for biosynthesis of the potential anti-cancer drug sansalvamide in Fusarium solani

Sansalvamide is a cyclic pentadepsipeptide produced by Fusarium solani and has shown promising results as potential anti-cancer drug. The biosynthetic pathway has until now remained unidentified, but here we used an Agrobacterium tumefaciens-mediated transformation (ATMT) approach to generate knockout mutants of two candidate non-ribosomal peptide synthetases (NRPS29 and NRPS30). Comparative studies of secondary metabolites in the two deletion mutants and wild type confirmed the absence of sansalvamide in the NRPS30 deletion mutant, implicating this synthetase in the biosynthetic pathway for sansalvamide. Sansalvamide is structurally related to the cyclic hexadepsipeptide destruxin, which both contain an α-hydroxyisocaproic acid (HICA) unit. A gene cluster responsible for destruxin production has previously been identified in Metarhizium robertsii together with a hypothetical biosynthetic pathway. Using comparative bioinformatic analyses of the catalytic domains in the destruxin and sansalvamide NRPSs, we were able to propose a model for sansalvamide biosynthesis. Orthologues of the gene clusters were also identified in species from several other genera including Acremonium chrysogenum and Trichoderma virens, which suggests that the ability to produce compounds related to destruxin and sansalvamide is widespread.

The secondary metabolite bioinformatics portal: Computational tools to facilitate synthetic biology of secondary metabolite production

Natural products are among the most important sources of lead molecules for drug discovery. With the development of affordable whole-genome sequencing technologies and other ‘omics tools, the field of natural products research is currently undergoing a shift in paradigms. While, for decades, mainly analytical and chemical methods gave access to this group of compounds, nowadays genomics-based methods offer complementary approaches to find, identify and characterize such molecules. This paradigm shift also resulted in a high demand for computational tools to assist researchers in their daily work. In this context, this review gives a summary of tools and databases that currently are available to mine, identify and characterize natural product biosynthesis pathways and their producers based on ‘omics data. A web portal called Secondary Metabolite Bioinformatics Portal (SMBP at http://www.secondarymetabolites.org) is introduced to provide a one-stop catalog and links to these bioinformatics resources. In addition, an outlook is presented how the existing tools and those to be developed will influence synthetic biology approaches in the natural products field.

Characterization of Nonribosomal Peptide Synthetases with NRPSsp

Bioinformatic sequence analysis allows the functional characterization of newly sequenced proteins. Nonribosomal peptide synthetases (NRPSs) are multi-modular enzymes involved in the biosynthesis of natural products. The current omics era has enabled the exponential growth of the sequenced NRPS, and it is important to characterize the final product of these synthetases. Here, how to achieve the prediction of substrates which bind to adenylation domains in NRPS with NRPSsp (www.nrpssp.com) bioinformatic tool is described.

Identification of the Polyketide Biosynthetic Machinery for the Indolizidine Alkaloid Cyclizidine

The cyclizidine biosynthetic gene cluster was identified from Streptomyces NCIB 11649, which revealed the polyketide biosynthetic machinery for cyclizidine alkaloid biosynthesis. Both in vivo mutagenesis study and in vitro biochemical analysis provided insight into the timing and mechanism of the biosynthetic enzymes that produce cyclizidine-type indolizidine alkaloids.

Genomes to natural products PRediction Informatics for Secondary Metabolomes (PRISM)

Microbial natural products are an invaluable source of evolved bioactive small molecules and pharmaceutical agents. Next-generation and metagenomic sequencing indicates untapped genomic potential, yet high rediscovery rates of known metabolites increasingly frustrate conventional natural product screening programs. New methods to connect biosynthetic gene clusters to novel chemical scaffolds are therefore critical to enable the targeted discovery of genetically encoded natural products. Here, we present PRISM, a computational resource for the identification of biosynthetic gene clusters, prediction of genetically encoded nonribosomal peptides and type I and II polyketides, and bio- and cheminformatic dereplication of known natural products. PRISM implements novel algorithms which render it uniquely capable of predicting type II polyketides, deoxygenated sugars, and starter units, making it a comprehensive genome-guided chemical structure prediction engine. A library of 57 tailoring reactions is leveraged for combinatorial scaffold library generation when multiple potential substrates are consistent with biosynthetic logic. We compare the accuracy of PRISM to existing genomic analysis platforms. PRISM is an open-source, user-friendly web application available at http://magarveylab.ca/prism/.

An automated Genomes-to-Natural Products platform (GNP) for the discovery of modular natural products

Bacterial natural products are a diverse and valuable group of small molecules, and genome sequencing indicates that the vast majority remain undiscovered. The prediction of natural product structures from biosynthetic assembly lines can facilitate their discovery, but highly automated, accurate, and integrated systems are required to mine the broad spectrum of sequenced bacterial genomes. Here we present a genome-guided natural products discovery tool to automatically predict, combinatorialize and identify polyketides and nonribosomal peptides from biosynthetic assembly lines using LC–MS/MS data of crude extracts in a high-throughput manner. We detail the directed identification and isolation of six genetically predicted polyketides and nonribosomal peptides using our Genome-to-Natural Products platform. This highly automated, user-friendly programme provides a means of realizing the potential of genetically encoded natural products.

Microbial natural products represent a large reservoir of potential pharmaceutical agents. Here, Johnston et al. describe a computer-automated programme for connecting genome sequences with identified and isolated natural products.

Prediction of the substrate for nonribosomal peptide synthetase (NRPS) adenylation domains by virtual screening

Nonribosomal peptide synthetases (NRPSs) synthesize a diverse array of bioactive small peptides, many of which are used in medicine. There is considerable interest in predicting NRPS substrate specificity in order to facilitate investigation of the many "cryptic" NRPS genes that have not been linked to any known product. However, the current sequence similarity-based methods are unable to produce reliable predictions when there is a lack of prior specificity data, which is a particular problem for fungal NRPSs. We conducted virtual screening on the specificity-determining domain of NRPSs, the adenylation domain, and found that virtual screening using experimentally determined structures results in good enrichment of the cognate substrate. Our results indicate that the conformation of the adenylation domain and in particular the conformation of a key conserved aromatic residue is important in determining the success of the virtual screening. When homology models of NRPS adenylation domains of known specificity, rather than experimentally determined structures, were built and used for virtual screening, good enrichment of the cognate substrate was also achieved in many cases. However, the accuracy of the models was key to the reliability of the predictions and there was a large variation in the results when different models of the same domain were used. This virtual screening approach is promising and is able to produce enrichment of the cognate substrates in many cases, but improvements in building and assessing homology models are required before the approach can be reliably applied to these models.