Distribution and diversity of 'Tectomicrobia', a deep-branching uncultivated bacterial lineage harboring rich producers of bioactive metabolites.

Genomic and functional analyses of bacterial sponge symbionts belonging to the uncultivated candidate genus 'Entotheonella' has revealed them as the prolific producers of bioactive compounds previously identified from their invertebrate hosts. These studies also suggested 'Entotheonella' as the first members of a new candidate phylum, 'Tectomicrobia'. Here we analyzed the phylogenetic structure and environmental distribution of this as-yet sparsely populated phylum-like lineage. The data show that 'Entotheonella' and other 'Tectomicrobia' are not restricted to marine habitats but widely distributed among terrestrial locations. The inferred phylogenetic trees suggest several intra-phylum lineages with diverse lifestyles. Of these, the previously described 'Entotheonella' lineage can be more accurately divided into at least three different candidate genera with the terrestrial 'Candidatus Prasianella', the largely terrestrial 'Candidatus Allonella', the 'Candidatus Thalassonella' comprising sponge-associated members, and the more widely distributed 'Candidatus Entotheonella'. Genomic characterization of 'Thalassonella' members from a range of sponge hosts did not suggest a role as providers of natural products, despite high genomic similarity to 'Entotheonella' regarding primary metabolism and implied lifestyle. In contrast, the analysis revealed a correlation between the revised 'Entotheonella' 16S rRNA gene phylogeny and a specific association with sponges and their natural products. This feature might serve as a discovery method to accelerate the identification of new chemically rich 'Entotheonella' variants, and led to the identification of the first 'Entotheonella' symbiont in a non-tetractinellid sponge, Psammocinia sp., indicating a wide host distribution of 'Entotheonella'-based chemical symbiosis.

Marine Sponge Endosymbionts: Structural and Functional Specificity of the Microbiome within Euryspongia arenaria Cells.

Sponge microbiomes are typically profiled by analyzing the community DNA of whole tissues, which does not distinguish the taxa residing within sponge cells from extracellular microbes. To uncover the endosymbiotic microbiome, we separated the sponge cells to enrich the intracellular microbes. The intracellular bacterial community of sponge Euryspongia arenaria was initially assessed by amplicon sequencing, which indicated that it hosts three unique phyla not found in the extracellular and bulk tissue microbiomes. These three phyla account for 66% of the taxonomically known genera in the intracellular microbiome. The shotgun metagenomic analysis extended the taxonomic coverage to viruses and eukaryotes, revealing the most abundant signature taxa specific to the intracellular microbiome. Functional KEGG pathway annotation demonstrated that the endosymbiotic microbiome hosted the greatest number of unique gene orthologs. The pathway profiles distinguished the intra- and extracellular microbiomes from the tissue and seawater microbiomes. Carbohydrate-active enzyme analysis further discriminated each microbiome based on their representative and dominant enzyme families. One pathway involved in digestion system and family esterase had a consistently higher level in intracellular microbiome and could statistically differentiate the intracellular microbiome from the others, suggesting that triacylglycerol lipases could be the key functional component peculiar to the endosymbionts. The identified higher abundance of lipase-related eggNOG categories further supported the lipid-hydrolyzing metabolism of endosymbiotic microbiota. Pseudomonas members, reported as lipase-producing bacteria, were only in the endosymbiotic microbiome, meanwhile Pseudomonas also showed a greater abundance intracellularly. Our study aided a comprehensive sponge microbiome that demonstrated the taxonomic and functional specificity of endosymbiotic microbiota. IMPORTANCE Sponges host abundant microbial symbionts that can produce an impressive number of novel bioactive metabolites. However, knowledge on intracellular (endosymbiotic) microbiota is scarce. We characterize the composition and function of the endosymbiotic microbiome by separation of sponge cells and enrichment of intracellular microbes. We uncover a noteworthy number of taxa exclusively in the endosymbiotic microbiome. We unlock the unique pathways and enzymes of endosymbiotic taxa. This study achieves a more comprehensive sponge microbial community profile, which demonstrates the structural and functional specificity of the endosymbiotic microbiome. Our findings not only open the possibility to reveal the low abundant and the likely missed microbiota when directly sequencing the sponge bulk tissues, but also warrant future in-depth exploration within single sponge cells.

Single-cell metabolite detection and genomics reveals uncultivated talented producer.

The production of bioactive metabolites is increasingly recognized as an important function of host-associated bacteria. An example is defensive symbiosis that might account for much of the chemical richness of marine invertebrates including sponges (Porifera), 1 of the oldest metazoans. However, most bacterial members of sponge microbiomes have not been cultivated or sequenced, and therefore, remain unrecognized. Unequivocally linking metabolic functions to a cellular source in sponge microbiomes is, therefore, a challenge. Here, we report an analysis pipeline of microfluidic encapsulation, Raman microscopy, and integrated digital genomics (MERMAID) for an efficient identification of uncultivated producers. We applied this method to the chemically rich bacteriosponge (sponge that hosts a rich bacterial community) Theonella swinhoei, previously shown to contain 'Entotheonella' symbionts that produce most of the bioactive substances isolated from the sponge. As an exception, the antifungal aurantosides had remained unassigned to a source. Raman-guided single-bacterial analysis and sequencing revealed a cryptic, distinct multiproducer, 'Candidatus Poriflexus aureus' from a new Chloroflexi lineage as the aurantoside producer. Its exceptionally large genome contains numerous biosynthetic loci and suggested an even higher chemical richness of this sponge than previously appreciated. This study highlights the importance of complementary technologies to uncover microbiome functions, reveals remarkable parallels between distantly related symbionts of the same host, and adds functional support for diverse chemically prolific lineages being present in microbial dark matter.

Opening up the Single-Cell Toolbox for Microbial Natural Products Research.

The diverse microbes that produce natural products represent an important source of novel therapeutics, drug leads, and scientific tools. However, the vast majority have not been grown in axenic culture and are members of complex communities. While meta-'omic methods such as metagenomics, -transcriptomics, and -proteomics reveal collective molecular features of this "microbial dark matter", the study of individual microbiome members can be challenging. To address these limits, a number of techniques with single-bacterial resolution have been developed in the last decade and a half. While several of these are embraced by microbial ecologists, there has been less use by researchers interested in mining microbes for natural products. In this review, we discuss the available and emerging techniques for targeted single-cell analysis with a particular focus on applications to the discovery and study of natural products.

A monodomain class II terpene cyclase assembles complex isoprenoid scaffolds.

Class II terpene cyclases, such as oxidosqualene and squalene-hopene cyclases, catalyse some of the most complex polycyclization reactions. They minimally exhibit a ß,γ-didomain architecture that has been evolutionarily repurposed in a wide range of terpene-processing enzymes and likely resulted from a fusion of unidentified monodomain proteins. Although single domain class I terpene cyclases have already been identified, the corresponding class II counterparts have not been previously reported. Here we present high-resolution X-ray structures of a monodomain class II cyclase, merosterolic acid synthase (MstE). With a minimalistic ß-domain architecture, this cyanobacterial enzyme is able to construct four rings in cytotoxic meroterpenoids with a sterol-like topology. The structures with bound substrate, product, and inhibitor provide detailed snapshots of a cyclization mechanism largely governed by residues located in a noncanonical enzyme region. Our results complement the few known class II cyclase crystal structures, while also indicating that archaic monodomain cyclases might have already catalyzed complex reaction cascades.

Characterization of an Orphan Type†III Polyketide Synthase Conserved in Uncultivated "Entotheonella" Sponge Symbionts.

Uncultivated bacterial symbionts from the candidate genus "Entotheonella" have been shown to produce diverse natural products previously attributed to their sponge hosts. In addition to these known compounds, "Entotheonella" genomes contain rich sets of biosynthetic gene clusters that lack identified natural products. Among these is a small type III polyketide synthase (PKS) cluster, one of only three clusters present in all known "Entotheonella" genomes. This conserved "Entotheonella" PKS (cep) cluster encodes the type III PKS CepA and the putative methyltransferase CepB. Herein, the characterization of CepA as an enzyme involved in phenolic lipid biosynthesis is reported. In vitro analysis showed a specificity for alkyl starter substrates and the production of tri- and tetraketide pyrones and tetraketide resorcinols. The conserved distribution of the cep cluster suggests an important role for the phenolic lipid polyketides produced in "Entotheonella" variants.

The impact of cryosolution thermal contraction on proteins and protein crystals: volumes, conformation and order.

Cryocooling of macromolecular crystals is commonly employed to limit radiation damage during X-ray diffraction data collection. However, cooling itself affects macromolecular conformation and often damages crystals via poorly understood processes. Here, the effects of cryosolution thermal contraction on macromolecular conformation and crystal order in crystals ranging from 32 to 67% solvent content are systematically investigated. It is found that the solution thermal contraction affects macromolecule configurations and volumes, unit-cell volumes, crystal packing and crystal order. The effects occur through not only thermal contraction, but also pressure caused by the mismatched contraction of cryosolvent and pores. Higher solvent-content crystals are more affected. In some cases the solvent contraction can be adjusted to reduce mosaicity and increase the strength of diffraction. Ice formation in some crystals is found to cause damage via a reduction in unit-cell volume, which is interpreted through solvent transport out of unit cells during cooling. The results point to more deductive approaches to cryoprotection optimization by adjusting the cryosolution composition to reduce thermal contraction-induced stresses in the crystal with cooling.

Directed Evolution Mimics Allosteric Activation by Stepwise Tuning of the Conformational Ensemble.

Allosteric enzymes contain a wealth of catalytic diversity that remains distinctly underutilized for biocatalysis. Tryptophan synthase is a model allosteric system and a valuable enzyme for the synthesis of noncanonical amino acids (ncAA). Previously, we evolved the ß-subunit from Pyrococcus furiosus, PfTrpB, for ncAA synthase activity in the absence of its native partner protein PfTrpA. However, the precise mechanism by which mutation activated TrpB to afford a stand-alone catalyst remained enigmatic. Here, we show that directed evolution caused a gradual change in the rate-limiting step of the catalytic cycle. Concomitantly, the steady-state distribution of the intermediates shifts to favor covalently bound Trp adducts, which have increased thermodynamic stability. The biochemical properties of these evolved, stand-alone TrpBs converge on those induced in the native system by allosteric activation. High-resolution crystal structures of the wild-type enzyme, an intermediate in the lineage, and the final variant, encompassing five distinct chemical states, show that activating mutations have only minor structural effects on their immediate environment. Instead, mutation stabilizes the large-scale motion of a subdomain to favor an otherwise transiently populated closed conformational state. This increase in stability enabled the first structural description of Trp covalently bound in a catalytically active TrpB, confirming key features of catalysis. These data combine to show that sophisticated models of allostery are not a prerequisite to recapitulating its complex effects via directed evolution, opening the way to engineering stand-alone versions of diverse allosteric enzymes.

Single-bacterial genomics validates rich and varied specialized metabolism of uncultivated Entotheonella sponge symbionts.

Marine sponges are prolific sources of unique bioactive natural products. The sponge Theonella swinhoei is represented by several distinct variants with largely nonoverlapping chemistry. For the Japanese chemotype Y harboring diverse complex polyketides and peptides, we previously provided genomic and functional evidence that a single symbiont, the filamentous, multicellular organism "Candidatus Entotheonella factor," produces almost all of these compounds. To obtain further insights into the chemistry of "Entotheonella," we investigated another phylotype, "Candidatus Entotheonella serta," present in the T. swinhoei WA sponge chemotype, a source of theonellamide- and misakinolide-type compounds. Unexpectedly, considering the lower chemical diversity, sequencing of individual bacterial filaments revealed an even larger number of biosynthetic gene regions than for Ca E. factor, with virtually no overlap. These included genes for misakinolide and theonellamide biosynthesis, the latter assigned by comparative genomic and metabolic analysis of a T. swinhoei chemotype from Israel, and by biochemical studies. The data suggest that both compound families, which were among the earliest model substances to study bacterial producers in sponges, originate from the same bacterium in T. swinhoei WA. They also add evidence that metabolic richness and variability could be a more general feature of Entotheonella symbionts.

Enzyme Nicotinamide Cofactor Specificity Reversal Guided by Automated Structural Analysis and Library Design.

The specificity of enzymes for nicotinamide adenine dinucleotide (NAD) or nicotinamide adenine dinucleotide phosphate (NADP) as redox carriers can pose a significant hurdle for metabolic engineering and synthetic biology applications, where switching the specificity might be beneficial. We have developed an easy-to-use computational tool (CSR-SALAD) for the design of mutant libraries to simplify the process of reversing the cofactor specificity of an enzyme. Here, we describe the optimal use of this tool and present methods for its application in a laboratory setting.

4S-Hydroxylation of Insulin at ProB28 Accelerates Hexamer Dissociation and Delays Fibrillation.

Daily injections of insulin provide lifesaving benefits to millions of diabetics. But currently available prandial insulins are suboptimal: The onset of action is delayed by slow dissociation of the insulin hexamer in the subcutaneous space, and insulin forms amyloid fibrils upon storage in solution. Here we show, through the use of noncanonical amino acid mutagenesis, that replacement of the proline residue at position 28 of the insulin B-chain (ProB28) by (4S)-hydroxyproline (Hzp) yields an active form of insulin that dissociates more rapidly, and fibrillates more slowly, than the wild-type protein. Crystal structures of dimeric and hexameric insulin preparations suggest that a hydrogen bond between the hydroxyl group of Hzp and a backbone amide carbonyl positioned across the dimer interface may be responsible for the altered behavior. The effects of hydroxylation are stereospecific; replacement of ProB28 by (4R)-hydroxyproline (Hyp) causes little change in the rates of fibrillation and hexamer disassociation. These results demonstrate a new approach that fuses the concepts of medicinal chemistry and protein design, and paves the way to further engineering of insulin and other therapeutic proteins.

Directed Evolution of a Bright Near-Infrared Fluorescent Rhodopsin Using a Synthetic Chromophore.

By engineering a microbial rhodopsin, Archaerhodopsin-3 (Arch), to bind a synthetic chromophore, merocyanine retinal, in place of the natural chromophore all-trans-retinal (ATR), we generated a protein with exceptionally bright and unprecedentedly red-shifted near-infrared (NIR) fluorescence. We show that chromophore substitution generates a fluorescent Arch complex with a 200-nm bathochromic excitation shift relative to ATR-bound wild-type Arch and an emission maximum at 772 nm. Directed evolution of this complex produced variants with pH-sensitive NIR fluorescence and molecular brightness 8.5-fold greater than the brightest ATR-bound Arch variant. The resulting proteins are well suited to bacterial imaging; expression and stability have not been optimized for mammalian cell imaging. By targeting both the protein and its chromophore, we overcome inherent challenges associated with engineering bright NIR fluorescence into Archaerhodopsin. This work demonstrates an efficient strategy for engineering non-natural, tailored properties into microbial opsins, properties relevant for imaging and interrogating biological systems.

A General Tool for Engineering the NAD/NADP Cofactor Preference of Oxidoreductases.

The ability to control enzymatic nicotinamide cofactor utilization is critical for engineering efficient metabolic pathways. However, the complex interactions that determine cofactor-binding preference render this engineering particularly challenging. Physics-based models have been insufficiently accurate and blind directed evolution methods too inefficient to be widely adopted. Building on a comprehensive survey of previous studies and our own prior engineering successes, we present a structure-guided, semirational strategy for reversing enzymatic nicotinamide cofactor specificity. This heuristic-based approach leverages the diversity and sensitivity of catalytically productive cofactor binding geometries to limit the problem to an experimentally tractable scale. We demonstrate the efficacy of this strategy by inverting the cofactor specificity of four structurally diverse NADP-dependent enzymes: glyoxylate reductase, cinnamyl alcohol dehydrogenase, xylose reductase, and iron-containing alcohol dehydrogenase. The analytical components of this approach have been fully automated and are available in the form of an easy-to-use web tool: Cofactor Specificity Reversal-Structural Analysis and Library Design (CSR-SALAD).

Discovery of a regioselectivity switch in nitrating P450s guided by molecular dynamics simulations and Markov models.

The dynamic motions of protein structural elements, particularly flexible loops, are intimately linked with diverse aspects of enzyme catalysis. Engineering of these loop regions can alter protein stability, substrate binding and even dramatically impact enzyme function. When these flexible regions are unresolvable structurally, computational reconstruction in combination with large-scale molecular dynamics simulations can be used to guide the engineering strategy. Here we present a collaborative approach that consists of both experiment and computation and led to the discovery of a single mutation in the F/G loop of the nitrating cytochrome P450 TxtE that simultaneously controls loop dynamics and completely shifts the enzyme's regioselectivity from the C4 to the C5 position of L-tryptophan. Furthermore, we find that this loop mutation is naturally present in a subset of homologous nitrating P450s and confirm that these uncharacterized enzymes exclusively produce 5-nitro-L-tryptophan, a previously unknown biosynthetic intermediate.

Artificial domain duplication replicates evolutionary history of ketol-acid reductoisomerases.

The duplication of protein structural domains has been proposed as a common mechanism for the generation of new protein folds. A particularly interesting case is the class II ketol-acid reductoisomerase (KARI), which putatively arose from an ancestral class I KARI by duplication of the C-terminal domain and corresponding loss of obligate dimerization. As a result, the class II enzymes acquired a deeply embedded figure-of-eight knot. To test this evolutionary hypothesis we constructed a novel class II KARI by duplicating the C-terminal domain of a hyperthermostable class I KARI. The new protein is monomeric, as confirmed by gel filtration and X-ray crystallography, and has the deeply knotted class II KARI fold. Surprisingly, its catalytic activity is nearly unchanged from the parent KARI. This provides strong evidence in support of domain duplication as the mechanism for the evolution of the class II KARI fold and demonstrates the ability of domain duplication to generate topological novelty in a function-neutral manner.

Cofactor specificity motifs and the induced fit mechanism in class I ketol-acid reductoisomerases.

Although most sequenced members of the industrially important ketol-acid reductoisomerase (KARI) family are class I enzymes, structural studies to date have focused primarily on the class II KARIs, which arose through domain duplication. In the present study, we present five new crystal structures of class I KARIs. These include the first structure of a KARI with a six-residue ß2αB (cofactor specificity determining) loop and an NADPH phosphate-binding geometry distinct from that of the seven- and 12-residue loops. We also present the first structures of naturally occurring KARIs that utilize NADH as cofactor. These results show insertions in the specificity loops that confounded previous attempts to classify them according to loop length. Lastly, we explore the conformational changes that occur in class I KARIs upon binding of cofactor and metal ions. The class I KARI structures indicate that the active sites close upon binding NAD(P)H, similar to what is observed in the class II KARIs of rice and spinach and different from the opening of the active site observed in the class II KARI of Escherichia coli. This conformational change involves a decrease in the bending of the helix that runs between the domains and a rearrangement of the nicotinamide-binding site.

Structural, functional, and spectroscopic characterization of the substrate scope of the novel nitrating cytochrome P450 TxtE.

A novel cytochrome P450 enzyme, TxtE, was recently shown to catalyze the direct aromatic nitration of L-tryptophan. This unique chemistry inspired us to ask whether TxtE could serve as a platform for engineering new nitration biocatalysts to replace current harsh synthetic methods. As a first step toward this goal, and to better understand the wild-type enzyme, we obtained high-resolution structures of TxtE in its substrate-free and substrate-bound forms. We also screened a library of substrate analogues for spectroscopic indicators of binding and for production of nitrated products. From these results, we found that the wild-type enzyme accepts moderate decoration of the indole ring, but the amino acid moiety is crucial for binding and correct positioning of the substrate and therefore less amenable to modification. A nitrogen atom is essential for catalysis, and a carbonyl must be present to recruit the αB'1 helix of the protein to seal the binding pocket.

General approach to reversing ketol-acid reductoisomerase cofactor dependence from NADPH to NADH.

To date, efforts to switch the cofactor specificity of oxidoreductases from nicotinamide adenine dinucleotide phosphate (NADPH) to nicotinamide adenine dinucleotide (NADH) have been made on a case-by-case basis with varying degrees of success. Here we present a straightforward recipe for altering the cofactor specificity of a class of NADPH-dependent oxidoreductases, the ketol-acid reductoisomerases (KARIs). Combining previous results for an engineered NADH-dependent variant of Escherichia coli KARI with available KARI crystal structures and a comprehensive KARI-sequence alignment, we identified key cofactor specificity determinants and used this information to construct five KARIs with reversed cofactor preference. Additional directed evolution generated two enzymes having NADH-dependent catalytic efficiencies that are greater than the wild-type enzymes with NADPH. High-resolution structures of a wild-type/variant pair reveal the molecular basis of the cofactor switch.