Fast and accurate protein structure search with Foldseek.

As structure prediction methods are generating millions of publicly available protein structures, searching these databases is becoming a bottleneck. Foldseek aligns the structure of a query protein against a database by describing tertiary amino acid interactions within proteins as sequences over a structural alphabet. Foldseek decreases computation times by four to five orders of magnitude with 86%, 88% and 133% of the sensitivities of Dali, TM-align and CE, respectively.

Clustering predicted structures at the scale of the known protein universe.

Proteins are key to all cellular processes and their structure is important in understanding their function and evolution. Sequence-based predictions of protein structures have increased in accuracy1, and over 214 million predicted structures are available in the AlphaFold database2. However, studying protein structures at this scale requires highly efficient methods. Here, we developed a structural-alignment-based clustering algorithm-Foldseek cluster-that can cluster hundreds of millions of structures. Using this method, we have clustered all of the structures in the AlphaFold database, identifying 2.30 million non-singleton structural clusters, of which 31% lack annotations representing probable previously undescribed structures. Clusters without annotation tend to have few representatives covering only 4% of all proteins in the AlphaFold database. Evolutionary analysis suggests that most clusters are ancient in origin but 4% seem to be species specific, representing lower-quality predictions or examples of de novo gene birth. We also show how structural comparisons can be used to predict domain families and their relationships, identifying examples of remote structural similarity. On the basis of these analyses, we identify several examples of human immune-related proteins with putative remote homology in prokaryotic species, illustrating the value of this resource for studying protein function and evolution across the tree of life.

Stereodivergent Hydroxylation of Berkeleylactones by Penicillium turbatum.

Penicillium turbatum has previously been reported to produce A26771B, a 16-membered macrocyclic polyketide with activity against Gram-positive bacteria, mycoplasma, and fungi, as well as the structurally related compounds berkeleylactone E and berkeleylactones I-O. In this work, large-scale cultivation of P. turbatum NRRL 5630 on rice yielded seven new berkeleylactone analogues, berkeleylactone E methyl ester, 14-epi-berkeleylactone F, berkeleylactones P-R, 12-epi-berkeleylactone Q, and 13-epi-berkeleylactone R, and six previously reported analogues, A26771B and berkeleylactones E-G and J-K. The structures of the berkeleylactones were elucidated by detailed analysis of spectroscopic data, molecular modeling, and comparison with literature values. Interestingly, six of the berkeleylactone analogues were isolated as pairs of hydroxy epimers, highlighting how Nature can exploit stereodivergence in biosynthetic pathways to increase chemical diversity. The genome of P. turbatum was sequenced, and a putative gene cluster (bekl) responsible for the biosynthesis of the berkeleylactones was identified. The new berkeleylactone analogues exhibited no significant biological activity against a panel of bacteria, fungi, the parasite Giardia duodenalis, or NS-1 murine myeloma cells, suggesting a hitherto undiscovered biological role.

UFCG: database of universal fungal core genes and pipeline for genome-wide phylogenetic analysis of fungi.

In phylogenomics the evolutionary relationship of organisms is studied by their genomic information. A common approach to phylogenomics is to extract related genes from each organism, build a multiple sequence alignment and then reconstruct evolution relations through a phylogenetic tree. Often a set of highly conserved genes occurring in single-copy, called core genes, are used for this analysis, as they allow efficient automation within a taxonomic clade. Here we introduce the Universal Fungal Core Genes (UFCG) database and pipeline for genome-wide phylogenetic analysis of fungi. The UFCG database consists of 61 curated fungal marker genes, including a novel set of 41 computationally derived core genes and 20 canonical genes derived from literature, as well as marker gene sequences extracted from publicly available fungal genomes. Furthermore, we provide an easy-to-use, fully automated and open-source pipeline for marker gene extraction, training and phylogenetic tree reconstruction. The UFCG pipeline can identify marker genes from genomic, proteomic and transcriptomic data, while producing phylogenies consistent with those previously reported, and is publicly available together with the UFCG database at https://ufcg.steineggerlab.com.

Combining bio-telemetry and underwater imagery to elucidate the reproductive behaviour of a large, long-lived Australian freshwater teleost.

Murray cod Maccullochella peelii (Mitchell) have a key ecological role in ensuring the health of Australia's largest inland waterway, but many aspects surrounding its reproductive strategies in the wild are unknown. From 2015 to 2019 within the Northern Murray-Darling Basin, Australia, we used a combination of bio-telemetry and underwater imagery to quantify the behaviour of Murray cod across their breeding cycle in a natural riverine environment. In most years, breeding behaviour including nest site selection was observed from early-August and spawning from late-August through to late-October, which is considerably earlier than previously reported. There was a positive correlation between the onset of breeding behaviour and week-of-year, and spawning was correlated with moon-phase. Whilst some nesting sites were amongst woody debris and in hollow logs, the majority were located in shallow water on hard substrate underneath undercuts along the riverbank edge. Nests were frequently established in isolated and disconnected pools with little or no measurable flow, suggesting that river hydraulics is not a key component driving spawning of Murray cod across at least some areas of its range. Larvae were observed actively swimming and controlling their position within and near nests and used a scatter tactic when dispersing. We also established that disturbing nesting Murray cod had a negative impact on egg and larval survival. We suggest a review of current regulations to safeguard the long-term conservation of the species across all sections of its range.

Synthaser: a CD-Search enabled Python toolkit for analysing domain architecture of fungal secondary metabolite megasynth(et)ases.

BACKGROUND: Fungi are prolific producers of secondary metabolites (SMs), which are bioactive small molecules with important applications in medicine, agriculture and other industries. The backbones of a large proportion of fungal SMs are generated through the action of large, multi-domain megasynth(et)ases such as polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs). The structure of these backbones is determined by the domain architecture of the corresponding megasynth(et)ase, and thus accurate annotation and classification of these architectures is an important step in linking SMs to their biosynthetic origins in the genome. RESULTS: Here we report synthaser, a Python package leveraging the NCBI's conserved domain search tool for remote prediction and classification of fungal megasynth(et)ase domain architectures. Synthaser is capable of batch sequence analysis, and produces rich textual output and interactive visualisations which allow for quick assessment of the megasynth(et)ase diversity of a fungal genome. Synthaser uses a hierarchical rule-based classification system, which can be extensively customised by the user through a web application ( http://gamcil.github.io/synthaser ). We show that synthaser provides more accurate domain architecture predictions than comparable tools which rely on curated profile hidden Markov model (pHMM)-based approaches; the utilisation of the NCBI conserved domain database also allows for significantly greater flexibility compared to pHMM approaches. In addition, we demonstrate how synthaser can be applied to large scale genome mining pipelines through the construction of an Aspergillus PKS similarity network. CONCLUSIONS: Synthaser is an easy to use tool that represents a significant upgrade to previous domain architecture analysis tools. It is freely available under a MIT license from PyPI ( https://pypi.org/project/synthaser ) and GitHub ( https://github.com/gamcil/synthaser ).

clinker & clustermap.js: automatic generation of gene cluster comparison figures.

SUMMARY: Genes involved in biological pathways are often collocalised in gene clusters, the comparison of which can give valuable insights into their function and evolutionary history. However, comparison and visualization of gene cluster similarity is a tedious process, particularly when many clusters are being compared. Here, we present clinker, a Python based tool and clustermap.js, a companion JavaScript visualization library, which used together can automatically generate accurate, interactive, publication-quality gene cluster comparison figures directly from sequence files. AVAILABILITY AND IMPLEMENTATION: Source code and documentation for clinker and clustermap.js is available on GitHub (github.com/gamcil/clinker and github.com/gamcil/clustermap.js, respectively) under the MIT license. clinker can be installed directly from the Python Package Index via pip. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

cblaster: a remote search tool for rapid identification and visualization of homologous gene clusters.

Motivation: Genes involved in coordinated biological pathways, including metabolism, drug resistance and virulence, are often collocalized as gene clusters. Identifying homologous gene clusters aids in the study of their function and evolution, however, existing tools are limited to searching local sequence databases. Tools for remotely searching public databases are necessary to keep pace with the rapid growth of online genomic data. Results: Here, we present cblaster, a Python-based tool to rapidly detect collocated genes in local and remote databases. cblaster is easy to use, offering both a command line and a user-friendly graphical user interface. It generates outputs that enable intuitive visualizations of large datasets and can be readily incorporated into larger bioinformatic pipelines. cblaster is a significant update to the comparative genomics toolbox. Availability and implementation: cblaster source code and documentation is freely available from GitHub under the MIT license (github.com/gamcil/cblaster). Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Hancockiamides: phenylpropanoid piperazines from Aspergillus hancockii are biosynthesised by a versatile dual single-module NRPS pathway.

The hancockiamides are an unusual new family of N-cinnamoylated piperazines from the Australian soil fungus Aspergillus hancockii. Genomic analyses of A. hancockii identified a biosynthetic gene cluster (hkm) of 12 genes, including two single-module nonribosomal peptide synthetase (NRPS) genes. Heterologous expression of the hkm cluster in A. nidulans confirmed its role in the biosynthesis of the hancockiamides. We further demonstrated that a novel cytochrome P450, Hkm5, catalyses the methylenedioxy bridge formation, and that the PAL gene hkm12 is dispensable, but contributes to increased production of the cinnamoylated hancockiamides. In vitro enzymatic assays and substrate feeding studies demonstrated that NRPS Hkm11 activates and transfers trans-cinnamate to the piperazine scaffold and has flexibility to accept bioisosteric thienyl and furyl analogues. This is the first reported cinnamate-activating fungal NRPS. Expression of a truncated cluster lacking the acetyltransferase gene led to seven additional congeners, including an unexpected family of 2,5-dibenzylpiperazines. These pleiotropic effects highlight the plasticity of the pathway and the power of this approach for accessing novel natural product scaffolds.

Victorin, the host-selective cyclic peptide toxin from the oat pathogen Cochliobolus victoriae, is ribosomally encoded.

The necrotrophic fungal pathogen Cochliobolus victoriae produces victorin, a host-selective toxin (HST) essential for pathogenicity to certain oat cultivars with resistance against crown rust. Victorin is a mixture of highly modified heterodetic cyclic hexapeptides, previously assumed to be synthesized by a nonribosomal peptide synthetase. Herein, we demonstrate that victorin is a member of the ribosomally synthesized and posttranslationally modified peptide (RiPP) family of natural products. Analysis of a newly generated long-read assembly of the C. victoriae genome revealed three copies of precursor peptide genes (vicA1-3) with variable numbers of "GLKLAF" core peptide repeats corresponding to the victorin peptide backbone. vicA1-3 are located in repeat-rich gene-sparse regions of the genome and are loosely clustered with putative victorin biosynthetic genes, which are supported by the discovery of compact gene clusters harboring corresponding homologs in two distantly related plant-associated Sordariomycete fungi. Deletion of at least one copy of vicA resulted in strongly diminished victorin production. Deletion of a gene encoding a DUF3328 protein (VicYb) abolished the production altogether, supporting its predicted role in oxidative cyclization of the core peptide. In addition, we uncovered a copper amine oxidase (CAO) encoded by vicK, in which its deletion led to the accumulation of new glycine-containing victorin derivatives. The role of VicK in oxidative deamination of the N-terminal glycyl moiety of the hexapeptides to the active glyoxylate forms was confirmed in vitro. This study finally unraveled the genetic and molecular bases for biosynthesis of one of the first discovered HSTs and expanded our understanding of underexplored fungal RiPPs.

Comprehensive chemotaxonomic and genomic profiling of a biosynthetically talented Australian fungus, Aspergillus burnettii sp. nov.

Aspergillus burnettii is a new species belonging to the A. alliaceus clade in Aspergillus subgenus Circumdati section Flavi isolated from peanut-growing properties in southern Queensland, Australia. A. burnettii is a fast-growing, floccose fungus with distinctive brown conidia and is a talented producer of biomass-degrading enzymes and secondary metabolites. Chemical profiling of A. burnettii revealed the metabolites ochratoxin A, kotanins, isokotanins, asperlicin E, anominine and paspalinine, which are common to subgenus Circumdati, together with burnettiene A, burnettramic acids, burnettides, and high levels of 14α-hydroxypaspalinine and hirsutide. The genome of A. burnettii was sequenced and an annotated draft genome is presented. A. burnettii is rich in secondary metabolite biosynthetic gene clusters, containing 51 polyketide synthases, 28 non-ribosomal peptide synthetases and 19 genes related to terpene biosynthesis. Functional annotation of digestive enzymes of A. burnettii and A. alliaceus revealed overlapping carbon utilisation profiles, consistent with a close phylogenetic relationship.

CRISPR-Mediated Activation of Biosynthetic Gene Clusters for Bioactive Molecule Discovery in Filamentous Fungi.

Accessing the full biosynthetic potential encoded in the genomes of fungi is limited by the low expression of most biosynthetic gene clusters (BGCs) under common laboratory culture conditions. CRISPR-mediated transcriptional activation (CRISPRa) of fungal BGCs could accelerate genomics-driven bioactive secondary metabolite discovery. In this work, we established the first CRISPRa system for filamentous fungi. First, we constructed a CRISPR/dLbCas12a-VPR-based system and demonstrated the activation of a fluorescent reporter in Aspergillus nidulans. Then, we targeted the native nonribosomal peptide synthetase-like (NRPS-like) gene micA in both chromosomal and episomal contexts, achieving increased production of the compound microperfuranone. Finally, multigene CRISPRa led to the discovery of the mic cluster product as dehydromicroperfuranone. Additionally, we demonstrated the utility of the variant dLbCas12aD156R-VPR for CRISPRa at room temperature culture conditions. Different aspects that influence the efficiency of CRISPRa in fungi were investigated, providing a framework for the further development of fungal artificial transcription factors based on CRISPR/Cas.

Biosynthesis of a New Benzazepine Alkaloid Nanangelenin A from Aspergillus nanangensis Involves an Unusual l-Kynurenine-Incorporating NRPS Catalyzing Regioselective Lactamization.

1-Benzazepine is a pharmaceutically important scaffold but is rare among natural products. Nanangelenin A (1), containing an unprecedented 3,4-dihydro-1-benzazepine-2,5-dione-N-prenyl-N-acetoxy-anthranilamide scaffold, was isolated from a novel species of Australian fungus, Aspergillus nanangensis. Genomic and retrobiosynthetic analyses identified a putative nonribosomal peptide synthetase (NRPS) gene cluster (nan). The detailed biosynthetic pathway to 1 was established by heterologous pathway reconstitution in A. nidulans, which led to biosynthesis of intermediates nanagelenin B-F (2-5 and 7). We demonstrated that the NRPS NanA incorporates anthranilic acid (Ant) and l-kynurenine (l-Kyn), which is supplied by a dedicated indoleamine-2,3-dioxygenase NanC encoded in the gene cluster. Using heterologous in vivo assays and mutagenesis, we demonstrated that the C-terminal condensation (CT) and thiolation (T3) domains of NanA are responsible for the regioselective cyclization of the tethered Ant-l-Kyn dipeptide to form the unusual benzazepine scaffold in 1. We also showed that NanA-CT catalyzes the regioselective cyclization of a surrogate synthetic substrate, Ant-l-Kyn-N-acetylcysteamine, to give the benzazepine scaffold, while spontaneous cyclization of the dipeptide yielded the alternative kinetically favored benzodiazepine scaffold. The discovery of 1 and the characterization of NanA have expanded the chemical and functional diversities of fungal NRPSs.

Nanangenines: drimane sesquiterpenoids as the dominant metabolite cohort of a novel Australian fungus, Aspergillus nanangensis.

Chemical investigation of an undescribed Australian fungus, Aspergillus nanangensis, led to the identification of the nanangenines - a family of seven new and three previously reported drimane sesquiterpenoids. The structures of the nanangenines were elucidated by detailed spectroscopic analysis supported by single crystal X-ray diffraction studies. The compounds were assayed for in vitro activity against bacteria, fungi, mammalian cells and plants. Bioinformatics analysis, including comparative analysis with other acyl drimenol-producing Aspergilli, led to the identification of a putative nanangenine biosynthetic gene cluster that corresponds to the proposed biosynthetic pathway for nanangenines.

Discovery and Heterologous Biosynthesis of the Burnettramic Acids: Rare PKS-NRPS-Derived Bolaamphiphilic Pyrrolizidinediones from an Australian Fungus, Aspergillus burnettii.

The burnettramic acids are a new class of antibiotics from an Australian fungus Aspergillus burnettii. The rare bolaamphiphilic scaffold consists of ß-d-mannose linked to a pyrrolizidinedione unit via a 26-carbon chain. The most abundant metabolite displayed potent in vitro antifungal activity. Comparative genomics identified the hybrid PKS-NRPS bua gene cluster, which was verified by heterologous pathway reconstitution in Aspergillus nidulans.

Panning for gold in mould: can we increase the odds for fungal genome mining?

Fungi are a rich source of bioactive small molecules. However, the large number of biosynthetic gene clusters (BGCs) encoding these molecules in their genomes suggests their biosynthetic potential is far greater than we previously appreciated. The mining of fungal genomes therefore holds great promise for the discovery of new chemical entities for pharmaceutical and agricultural applications. As more and more fungal genomes become available, the accompanying number of BGCs is quickly becoming unmanageable. Along with improving molecular genetic tools to accelerate the translation of BGCs to small molecules, we must devise strategies to prioritise BGCs most likely to encode the biosynthesis of novel small molecules and molecules with new or improved bioactivities or functions. In this perspective, we discuss existing and emerging strategies for prioritisation of BGCs to increase the odds of fruitful genome mining in fungi.

The cutest little baby face: a hormonal link to sensitivity to cuteness in infant faces.

We used computer image manipulation to develop a test of perception of subtle gradations in cuteness between infant faces. We found that young women (19-26 years old) were more sensitive to differences in infant cuteness than were men (19-26 and 53-60 years old). Women aged 45 to 51 years performed at the level of the young women, whereas cuteness sensitivity in women aged 53 to 60 years was not different from that of men (19-26 and 53-60 years old). Because average age at menopause is 51 years in Britain, these findings suggest the possible involvement of reproductive hormones in cuteness sensitivity. Therefore, we compared cuteness discrimination in pre- and postmenopausal women matched for age and in women taking and not taking oral contraceptives (progestogen and estrogen). Premenopausal women and young women taking oral contraceptives (which raise hormone levels artificially) were more sensitive to variations of cuteness than their respective comparison groups. We suggest that cuteness sensitivity is modulated by female reproductive hormones.

Transient-state reduction and steady-state kinetic studies of menaquinol oxidase from Bacillus subtilis, cytochrome aa3-600 nm. Spectroscopic characterization of the steady-state species.

Cytochrome aa3-600 or menaquinol oxidase, from Bacillus subtilis, is a member of the heme-copper oxidase family. Cytochrome aa3-600 contains cytochrome a, cytochrome a3, and CuB, and each is coordinated via histidine residues to subunit I. Subunit II of cytochrome aa3-600 lacks CuA, which is a common feature of the cytochrome c oxidase family members. Anaerobic reduction of cytochrome aa3-600 by the substrate analogue 2,3-dimethyl-1,4-naphthoquinone (DMN) resolves two distinct kinetic phases by stopped-flow, single-wavelength spectrometry. Global analysis of time-resolved, multiwavelength spectra shows that during these distinct phases cytochromes a and a3 are both reduced. Cyanide binding to cytochrome a3 enhances the fast phase rate, which in the presence of cyanide can be assigned to cytochrome a reduction, whereas cytochrome a3-cyanide reduction is slow. The steady-state activity of cytochrome aa3-600 exhibits saturation kinetics as a function of DMN concentration with a Km of 300 microM and a maximal turnover of 63.5 s(-1). Global kinetic analysis of steady-state spectra reveals a species that is characteristic of a partially reduced oxygen adduct of cytochrome a3-CuB, whereas cytochrome a remains oxidized. Electron paramagnetic resonance (EPR) spectroscopy of the oxidase in the steady state shows the expected signal from ferricytochrome a, and a new EPR signal at g = 2.01. A model of the catalytic cycle for cytochrome aa3-600 proposes initial electron delivery from DMN to cytochrome a, followed by rapid heme to heme electron transfer, and suggests possible origins of the radical signal in the steady-state form of the enzyme.

Investigation of CO binding and release from Mo-nitrogenase during catalytic turnover.

During enzymatic turnover in the presence of CO, Mo-nitrogenase has been shown to generate two different EPR signals termed lo-CO (PCO = 0.08 atm) and hi-CO (PCO = 0.5 atm). When the formation of hi-CO is monitored under the conditions of very low electron flux, a 2 min lag is observed prior to the initial detection of the signal followed by a near-linear rate of formation during which the S = 3/2 cofactor signal exhibits similar decay kinetics. Increasing the electron flux produces a significant increase in the rate of both the formation of hi-CO and the decay of the S = 3/2 cofactor. These results are interpreted in terms of a state of the enzyme (redox or structural) generated only during turnover which is needed to initially bind CO to the cofactor. Under high electron flux conditions, new EPR inflections are observed at g = 5.78, 5.15 and g = 1.95, 1.81 and tentatively assigned to S = 3/2 and 1/2 states of the CO-bound cofactor and 1 equiv of oxidized P cluster, respectively. Sudden removal of CO from the environment results in the slow decay (>10 min) of both the hi-CO signal and CO inhibition of acetylene reduction activity. The use of ethylene glycol to quench enzymatic activity strongly inhibits the decay of hi-CO (in the presence of CO) and the subsequent decay of lo-CO (after removal of CO) but does not prevent the reversible interconversion hi-CO left and right arrow lo-CO + CO.

High throughput purification of combinatorial libraries.

This article summarizes recent advances at MDS Panlabs which provide for the high-throughput preparative HPLC purification of our Optiverse screening library. Topics covered include unique methods for the preparation, purification, characterization, and plating of purified screening compounds. Also discussed are procedures for data tracking as samples proceed through the purification process.