Lunaemycins, New Cyclic Hexapeptide Antibiotics from the Cave Moonmilk-Dweller Streptomyces lunaelactis MM109T.

Streptomyces lunaelactis strains have been isolated from moonmilk deposits, which are calcium carbonate speleothems used for centuries in traditional medicine for their antimicrobial properties. Genome mining revealed that these strains are a remarkable example of a Streptomyces species with huge heterogeneity regarding their content in biosynthetic gene clusters (BGCs) for specialized metabolite production. BGC 28a is one of the cryptic BGCs that is only carried by a subgroup of S. lunaelactis strains for which in silico analysis predicted the production of nonribosomal peptide antibiotics containing the non-proteogenic amino acid piperazic acid (Piz). Comparative metabolomics of culture extracts of S. lunaelactis strains either holding or not holding BGC 28a combined with MS/MS-guided peptidogenomics and 1H/13C NMR allowed us to identify the cyclic hexapeptide with the amino acid sequence (D-Phe)-(L-HO-Ile)-(D-Piz)-(L-Piz)-(D-Piz)-(L-Piz), called lunaemycin A, as the main compound synthesized by BGC 28a. Molecular networking further identified 18 additional lunaemycins, with 14 of them having their structure elucidated by HRMS/MS. Antimicrobial assays demonstrated a significant bactericidal activity of lunaemycins against Gram-positive bacteria, including multi-drug resistant clinical isolates. Our work demonstrates how an accurate in silico analysis of a cryptic BGC can highly facilitate the identification, the structural elucidation, and the bioactivity of its associated specialized metabolites.

AURTHO: Autoregulation of transcription factors as facilitator of cis-acting element discovery.

Transcriptional regulation is key in bacteria for providing an adequate response in time and space to changing environmental conditions. However, despite decades of research, the binding sites and therefore the target genes and the function of most transcription factors (TFs) remain unknown. Filling this gap in knowledge through conventional methods represents a colossal task which we demonstrate here can be significantly facilitated by a widespread feature in transcriptional control: the autoregulation of TFs implying that the yet unknown transcription factor binding site (TFBS) is neighboring the TF itself. In this work, we describe the "AURTHO" methodology (AUtoregulation of oRTHOlogous transcription factors), consisting of analyzing upstream regions of orthologous TFs in order to uncover their associated TFBSs. AURTHO enabled the de novo identification of novel TFBSs with an unprecedented improvement in terms of quantity and reliability. DNA-protein interaction studies on a selection of candidate cis-acting elements yielded an >90 % success rate, demonstrating the efficacy of AURTHO at highlighting true TF-TFBS couples and confirming the identification in a near future of a plethora of TFBSs across all bacterial species.

On the Risks of Phylogeny-Based Strain Prioritization for Drug Discovery: Streptomyces lunaelactis as a Case Study.

Strain prioritization for drug discovery aims at excluding redundant strains of a collection in order to limit the repetitive identification of the same molecules. In this work, we wanted to estimate what can be unexploited in terms of the amount, diversity, and novelty of compounds if the search is focused on only one single representative strain of a species, taking Streptomyces lunaelactis as a model. For this purpose, we selected 18 S. lunaelactis strains taxonomically clustered with the archetype strain S. lunaelactis MM109T. Genome mining of all S. lunaelactis isolated from the same cave revealed that 54% of the 42 biosynthetic gene clusters (BGCs) are strain specific, and five BGCs are not present in the reference strain MM109T. In addition, even when a BGC is conserved in all strains such as the bag/fev cluster involved in bagremycin and ferroverdin production, the compounds produced highly differ between the strains and previously unreported compounds are not produced by the archetype MM109T. Moreover, metabolomic pattern analysis uncovered important profile heterogeneity, confirming that identical BGC predisposition between two strains does not automatically imply chemical uniformity. In conclusion, trying to avoid strain redundancy based on phylogeny and genome mining information alone can compromise the discovery of new natural products and might prevent the exploitation of the best naturally engineered producers of specific molecules.

A Single Biosynthetic Gene Cluster Is Responsible for the Production of Bagremycin Antibiotics and Ferroverdin Iron Chelators.

Biosynthetic gene clusters (BGCs) are organized groups of genes involved in the production of specialized metabolites. Typically, one BGC is responsible for the production of one or several similar compounds with bioactivities that usually only vary in terms of strength and/or specificity. Here we show that the previously described ferroverdins and bagremycins, which are families of metabolites with different bioactivities, are produced from the same BGC, whereby the fate of the biosynthetic pathway depends on iron availability. Under conditions of iron depletion, the monomeric bagremycins are formed, representing amino-aromatic antibiotics resulting from the condensation of 3-amino-4-hydroxybenzoic acid with p-vinylphenol. Conversely, when iron is abundantly available, the biosynthetic pathway additionally produces a molecule based on p-vinylphenyl-3-nitroso-4-hydroxybenzoate, which complexes iron to form the trimeric ferroverdins that have anticholesterol activity. Thus, our work shows a unique exception to the concept that BGCs should only produce a single family of molecules with one type of bioactivity and that in fact different bioactive molecules may be produced depending on the environmental conditions.IMPORTANCE Access to whole-genome sequences has exposed the general incidence of the so-called cryptic biosynthetic gene clusters (BGCs), thereby renewing their interest for natural product discovery. As a consequence, genome mining is the often first approach implemented to assess the potential of a microorganism for producing novel bioactive metabolites. By revealing a new level of complexity of natural product biosynthesis, we further illustrate the difficulty of estimation of the panel of molecules associated with a BGC based on genomic information alone. Indeed, we found that the same gene cluster is responsible for the production of compounds which differ in terms of structure and bioactivity. The production of these different compounds responds to different environmental triggers, which suggests that multiplication of culture conditions is essential for revealing the entire panel of molecules made by a single BGC.

Complete Genome Sequence of Streptomyces lunaelactis MM109T, Isolated from Cave Moonmilk Deposits.

Streptomyces lunaelactis MM109T is a ferroverdin A (anticholesterol) producer isolated from cave moonmilk deposits. The complete genome sequence of MM109T was obtained by combining Oxford Nanopore MinION and Illumina HiSeq and MiSeq technologies, revealing an 8.4-Mb linear chromosome and two plasmids, pSLUN1 (127,264 bp, linear) and pSLUN2 (46,827 bp, circular).

Isolation, Characterization, and Antibacterial Activity of Hard-to-Culture Actinobacteria from Cave Moonmilk Deposits.

Cave moonmilk deposits host an abundant and diverse actinobacterial population that has a great potential for producing novel natural bioactive compounds. In our previous attempt to isolate culturable moonmilk-dwelling Actinobacteria, only Streptomyces species were recovered, whereas a metagenetic study of the same deposits revealed a complex actinobacterial community including 46 actinobacterial genera in addition to streptomycetes. In this work, we applied the rehydration-centrifugation method to lessen the occurrence of filamentous species and tested a series of strategies to achieve the isolation of hard-to-culture and rare Actinobacteria from the moonmilk deposits of the cave "Grotte des Collemboles". From the "tips and tricks" that were tested, separate autoclaving of the components of the International Streptomyces Project (ISP) medium number 5 (ISP5) medium, prolonged incubation time, and dilution of the moonmilk suspension were found to most effectively improve colony forming units. Taxonomic analyses of the 40 isolates revealed new representatives of the Agromyces, Amycolatopsis, Kocuria, Micrococcus, Micromonospora, Nocardia, and Rhodococcus species, as well as additional new streptomycetes. The applied methodologies allowed the isolation of strains associated with both the least and most abundant moonmilk-dwelling actinobacterial operational taxonomic units. Finally, bioactivity screenings revealed that some isolates displayed high antibacterial activities, and genome mining uncovered a strong potential for the production of natural compounds.

Cracking the regulatory code of biosynthetic gene clusters as a strategy for natural product discovery.

The World Health Organization (WHO) describes antibiotic resistance as "one of the biggest threats to global health, food security, and development today", as the number of multi- and pan-resistant bacteria is rising dangerously. Acquired resistance phenomena also impair antifungals, antivirals, anti-cancer drug therapy, while herbicide resistance in weeds threatens the crop industry. On the positive side, it is likely that the chemical space of natural products goes far beyond what has currently been discovered. This idea is fueled by genome sequencing of microorganisms which unveiled numerous so-called cryptic biosynthetic gene clusters (BGCs), many of which are transcriptionally silent under laboratory culture conditions, and by the fact that most bacteria cannot yet be cultivated in the laboratory. However, brute force antibiotic discovery does not yield the same results as it did in the past, and researchers have had to develop creative strategies in order to unravel the hidden potential of microorganisms such as Streptomyces and other antibiotic-producing microorganisms. Identifying the cis elements and their corresponding transcription factors(s) involved in the control of BGCs through bioinformatic approaches is a promising strategy. Theoretically, we are a few 'clicks' away from unveiling the culturing conditions or genetic changes needed to activate the production of cryptic metabolites or increase the production yield of known compounds to make them economically viable. In this opinion article, we describe and illustrate the idea beyond 'cracking' the regulatory code for natural product discovery, by presenting a series of proofs of concept, and discuss what still should be achieved to increase the rate of success of this strategy.

Assessment of the Potential Role of Streptomyces in Cave Moonmilk Formation.

Moonmilk is a karstic speleothem mainly composed of fine calcium carbonate crystals (CaCO3) with different textures ranging from pasty to hard, in which the contribution of biotic rock-building processes is presumed to involve indigenous microorganisms. The real microbial input in the genesis of moonmilk is difficult to assess leading to controversial hypotheses explaining the origins and the mechanisms (biotic vs. abiotic) involved. In this work, we undertook a comprehensive approach in order to assess the potential role of filamentous bacteria, particularly a collection of moonmilk-originating Streptomyces, in the genesis of this speleothem. Scanning electron microscopy (SEM) confirmed that indigenous filamentous bacteria could indeed participate in moonmilk development by serving as nucleation sites for CaCO3 deposition. The metabolic activities involved in CaCO3 transformation were furthermore assessed in vitro among the collection of moonmilk Streptomyces, which revealed that peptides/amino acids ammonification, and to a lesser extend ureolysis, could be privileged metabolic pathways participating in carbonate precipitation by increasing the pH of the bacterial environment. Additionally, in silico search for the genes involved in biomineralization processes including ureolysis, dissimilatory nitrate reduction to ammonia, active calcium ion transport, and reversible hydration of CO2 allowed to identify genetic predispositions for carbonate precipitation in Streptomyces. Finally, their biomineralization abilities were confirmed by environmental SEM, which allowed to visualize the formation of abundant mineral deposits under laboratory conditions. Overall, our study provides novel evidences that filamentous Actinobacteria could be key protagonists in the genesis of moonmilk through a wide spectrum of biomineralization processes.

A Phenotypic and Genotypic Analysis of the Antimicrobial Potential of Cultivable Streptomyces Isolated from Cave Moonmilk Deposits.

Moonmilk speleothems of limestone caves host a rich microbiome, among which Actinobacteria represent one of the most abundant phyla. Ancient medical texts reported that moonmilk had therapeutical properties, thereby suggesting that its filamentous endemic actinobacterial population might be a source of natural products useful in human treatment. In this work, a screening approach was undertaken in order to isolate cultivable Actinobacteria from moonmilk of the Grotte des Collemboles in Belgium, to evaluate their taxonomic profile, and to assess their potential in biosynthesis of antimicrobials. Phylogenetic analysis revealed that all 78 isolates were exclusively affiliated to the genus Streptomyces and clustered into 31 distinct phylotypes displaying various pigmentation patterns and morphological features. Phylotype representatives were tested for antibacterial and antifungal activities and their genomes were mined for secondary metabolite biosynthetic genes coding for non-ribosomal peptide synthetases (NRPSs), and polyketide synthases (PKS). The moonmilk Streptomyces collection was found to display strong inhibitory activities against a wide range of reference organisms, as 94, 71, and 94% of the isolates inhibited or impaired the growth of Gram-positive, Gram-negative bacteria, and fungi, respectively. Interestingly, 90% of the cave strains induced strong growth suppression against the multi-drug resistant Rasamsonia argillacea, a causative agent of invasive mycosis in cystic fibrosis and chronic granulomatous diseases. No correlation was observed between the global antimicrobial activity of an individual strain and the number of NRPS and PKS genes predicted in its genome, suggesting that approaches for awakening cryptic metabolites biosynthesis should be applied to isolates with no antimicrobial phenotype. Overall, our work supports the common belief that moonmilk might effectively treat various infectious diseases thanks to the presence of a highly diverse population of prolific antimicrobial producing Streptomyces, and thus may indeed constitute a promising reservoir of potentially novel active natural compounds.

A coarse-grained simulation study of the structures, energetics, and dynamics of linear and circular DNA with its ions.

We have studied the structural, energetics, and dynamical properties of a variety of linear and circular DNA fragments using a solvent-mediated coarse-grained (CG) model of DNA with explicit ions recently developed by us [Naômé et al., J. Chem. Theory Comput., 2014, 10, 3541-3549]. We particularly examined the treatment of electrostatics and determined that a large cutoff is necessary to properly reproduce the DNA flexibility. Moreover, it is crucial to include long-ranged electrostatic interactions: a Particle Mesh Ewald scheme at low resolution is sufficient to avoid structural artifacts. We calculated the ring closure probabilities, as j-factors, for DNA fragments of different lengths from equilibrium, as well as restrained molecular dynamics (MD) simulations. The latter force integration method provided accurate results without model fitting. We generated topology and energy maps for DNA minicircles of various lengths and helical densities, at low and high ion concentrations. A general trend for structure compaction is observed, driven by an increase in writhing as the ionic concentration increases. Finally, we applied a reconstruction procedure to generate detailed molecular structures from the various superhelical conformations generated by the CG MD of the DNA minicircles. These pre-equilibrated reconstructed atomistic structures can serve as starting material for atomistic simulations.

Systematic hierarchical coarse-graining with the inverse Monte Carlo method.

We outline our coarse-graining strategy for linking micro- and mesoscales of soft matter and biological systems. The method is based on effective pairwise interaction potentials obtained in detailed ab initio or classical atomistic Molecular Dynamics (MD) simulations, which can be used in simulations at less accurate level after scaling up the size. The effective potentials are obtained by applying the inverse Monte Carlo (IMC) method [A. P. Lyubartsev and A. Laaksonen, Phys. Rev. E 52(4), 3730-3737 (1995)] on a chosen subset of degrees of freedom described in terms of radial distribution functions. An in-house software package MagiC is developed to obtain the effective potentials for arbitrary molecular systems. In this work we compute effective potentials to model DNA-protein interactions (bacterial LiaR regulator bound to a 26 base pairs DNA fragment) at physiological salt concentration at a coarse-grained (CG) level. Normally the IMC CG pair-potentials are used directly as look-up tables but here we have fitted them to five Gaussians and a repulsive wall. Results show stable association between DNA and the model protein as well as similar position fluctuation profile.

A Solvent-Mediated Coarse-Grained Model of DNA Derived with the Systematic Newton Inversion Method.

We present a new class of coarse-grained (CG) force fields (FFs) for B-DNA with explicit ions suited for large-scale mesoscale simulations at microsecond-micrometer scale using a wide spectrum of particle simulation methods from molecular dynamics to dissipative particle dynamics. The effective solvent-mediated pairwise interactions making up the FFs are obtained by inverting radial distribution functions and other particle-particle distributions obtained from all-atom simulations of numbers of octadecamer DNA fragments from the Ascona B-DNA library. The inverse Monte Carlo (IMC) method, later known as Newton inversion (NI) (Lyubartsev, A. P.; Laaksonen, A. Phys. Rev. E, 1995, 52, 3730-3737), was used together with the iterative Boltzmann inversion (IBI) scheme to compute the effective CG potentials. We show that this systematic structure-based approach is capable of providing converged potentials that accurately reproduce the structural features of the underlying atomistic system within a few percents of relative difference. We also show that a simple one-site-per-nucleotide model with 10 intramolecular pair interaction potentials is able to reproduce key features of DNA, for example, the persistence length and its dependence on the ionic concentration, experimentally determined around 50 nm at physiological salt concentration.

Molecular dynamics simulation of 8-oxoguanine containing DNA fragments reveals altered hydration and ion binding patterns.

We present results of molecular dynamics (MD) simulations of a double-stranded DNA fragment of sequence 5'-d(CGCGAGTTCGCG)-3' and 3'-d(GCGCTCAAGCGC)-5' in its native form and altered with guanine in the G6/C19 base pair replaced by 8-oxoguanine (8oxoG). DNA was simulated in solution with Na(+) counterions and explicit water molecules. Trajectories were analyzed for radial distribution functions (RDFs), spatial distribution functions (SDFs), residence times, and occupancies, for both water and ions. Structural parameters of DNA fragments were also considered. Particularly, we report sensible differences in the statistical and dynamical behavior of water and ions between the healthy and lesion containing models. Structural features such as sugar-phosphate backbone torsion angles also significantly differ. We conclude that 8oxoG noticeably modifies its close environment and that the observed changes in the vicinity of the lesion site might help the enzymes in charge of cleansing 8oxoG from human DNA to locate the aberration in its intrahelical context. This work supports the existence of a determinant early recognition mode of 8oxoG for hOGG1.