Translational control plays an important role in the adaptive heat-shock response of Streptomyces coelicolor.

Stress-induced adaptations require multiple levels of regulation in all organisms to repair cellular damage. In the present study we evaluated the genome-wide transcriptional and translational changes following heat stress exposure in the soil-dwelling model actinomycete bacterium, Streptomyces coelicolor. The combined analysis revealed an unprecedented level of translational control of gene expression, deduced through polysome profiling, in addition to transcriptional changes. Our data show little correlation between the transcriptome and 'translatome'; while an obvious downward trend in genome wide transcription was observed, polysome associated transcripts following heat-shock showed an opposite upward trend. A handful of key protein players, including the major molecular chaperones and proteases were highly induced at both the transcriptional and translational level following heat-shock, a phenomenon known as 'potentiation'. Many other transcripts encoding cold-shock proteins, ABC-transporter systems, multiple transcription factors were more highly polysome-associated following heat stress; interestingly, these protein families were not induced at the transcriptional level and therefore were not previously identified as part of the stress response. Thus, stress coping mechanisms at the level of gene expression in this bacterium go well beyond the induction of a relatively small number of molecular chaperones and proteases in order to ensure cellular survival at non-physiological temperatures.

Correcting for link loss in causal network inference caused by regulator interference.

MOTIVATION: There are a number of algorithms to infer causal regulatory networks from time series (gene expression) data. Here we analyse the phenomena of regulator interference, where regulators with similar dynamics mutually suppress both the probability of regulating a target and the associated link strength; for instance, interference between two identical strong regulators reduces link probabilities by ∼50%. RESULTS: We construct a robust method to define an interference-corrected causal network based on an analysis of the conditional link probabilities that recovers links lost through interference. On a large real network (Streptomyces coelicolor, phosphate depletion), we demonstrate that significant interference can occur between regulators with a correlation as low as 0.865, losing an estimated 34% of links by interference. However, levels of interference cannot be predicted from the correlation between regulators alone and are data specific. Validating against known networks, we show that high numbers of functional links are lost by regulator interference. Performance against other methods on DREAM4 data is excellent. AVAILABILITY AND IMPLEMENTATION: The method is implemented in R and is publicly available as the NIACS package at http://www2.warwick.ac.uk/fac/sci/systemsbiology/research/software.

Metabolic switches and adaptations deduced from the proteomes of Streptomyces coelicolor wild type and phoP mutant grown in batch culture.

Bacteria in the genus Streptomyces are soil-dwelling oligotrophs and important producers of secondary metabolites. Previously, we showed that global messenger RNA expression was subject to a series of metabolic and regulatory switches during the lifetime of a fermentor batch culture of Streptomyces coelicolor M145. Here we analyze the proteome from eight time points from the same fermentor culture and, because phosphate availability is an important regulator of secondary metabolite production, compare this to the proteome of a similar time course from an S. coelicolor mutant, INB201 (ΔphoP), defective in the control of phosphate utilization. The proteomes provide a detailed view of enzymes involved in central carbon and nitrogen metabolism. Trends in protein expression over the time courses were deduced from a protein abundance index, which also revealed the importance of stress pathway proteins in both cultures. As expected, the ΔphoP mutant was deficient in expression of PhoP-dependent genes, and several putatively compensatory metabolic and regulatory pathways for phosphate scavenging were detected. Notably there is a succession of switches that coordinately induce the production of enzymes for five different secondary metabolite biosynthesis pathways over the course of the batch cultures.

Extracting regulator activity profiles by integration of de novo motifs and expression data: characterizing key regulators of nutrient depletion responses in Streptomyces coelicolor.

Determining transcriptional regulator activities is a major focus of systems biology, providing key insight into regulatory mechanisms and co-regulators. For organisms such as Escherichia coli, transcriptional regulator binding site data can be integrated with expression data to infer transcriptional regulator activities. However, for most organisms there is only sparse data on their transcriptional regulators, while their associated binding motifs are largely unknown. Here, we address the challenge of inferring activities of unknown regulators by generating de novo (binding) motifs and integrating with expression data. We identify a number of key regulators active in the metabolic switch, including PhoP with its associated directed repeat PHO box, candidate motifs for two SARPs, a CRP family regulator, an iron response regulator and that for LexA. Experimental validation for some of our predictions was obtained using gel-shift assays. Our analysis is applicable to any organism for which there is a reasonable amount of complementary expression data and for which motifs (either over represented or evolutionary conserved) can be identified in the genome.

The PII protein GlnK is a pleiotropic regulator for morphological differentiation and secondary metabolism in Streptomyces coelicolor.

GlnK is an important nitrogen sensor protein in Streptomyces coelicolor. Deletion of glnK results in a medium-dependent failure of aerial mycelium and spore formation and loss of antibiotic production. Thus, GlnK is not only a regulator of nitrogen metabolism but also of morphological differentiation and secondary metabolite production. Through a comparative transcriptomic approach between the S. coelicolor wild-type and a S. coelicolor glnK mutant strain, 142 genes were identified that are differentially regulated in both strains. Among these are genes of the ram and rag operon, which are involved in S. coelicolor morphogenesis, as well as genes involved in gas vesicle biosynthesis and ectoine biosynthesis. Surprisingly, no relevant nitrogen genes were found to be differentially regulated, revealing that GlnK is not an important nitrogen sensor under the tested conditions.

Metabolic modeling and analysis of the metabolic switch in Streptomyces coelicolor.

BACKGROUND: The transition from exponential to stationary phase in Streptomyces coelicolor is accompanied by a major metabolic switch and results in a strong activation of secondary metabolism. Here we have explored the underlying reorganization of the metabolome by combining computational predictions based on constraint-based modeling and detailed transcriptomics time course observations. RESULTS: We reconstructed the stoichiometric matrix of S. coelicolor, including the major antibiotic biosynthesis pathways, and performed flux balance analysis to predict flux changes that occur when the cell switches from biomass to antibiotic production. We defined the model input based on observed fermenter culture data and used a dynamically varying objective function to represent the metabolic switch. The predicted fluxes of many genes show highly significant correlation to the time series of the corresponding gene expression data. Individual mispredictions identify novel links between antibiotic production and primary metabolism. CONCLUSION: Our results show the usefulness of constraint-based modeling for providing a detailed interpretation of time course gene expression data.

The dynamic architecture of the metabolic switch in Streptomyces coelicolor.

BACKGROUND: During the lifetime of a fermenter culture, the soil bacterium S. coelicolor undergoes a major metabolic switch from exponential growth to antibiotic production. We have studied gene expression patterns during this switch, using a specifically designed Affymetrix genechip and a high-resolution time-series of fermenter-grown samples. RESULTS: Surprisingly, we find that the metabolic switch actually consists of multiple finely orchestrated switching events. Strongly coherent clusters of genes show drastic changes in gene expression already many hours before the classically defined transition phase where the switch from primary to secondary metabolism was expected. The main switch in gene expression takes only 2 hours, and changes in antibiotic biosynthesis genes are delayed relative to the metabolic rearrangements. Furthermore, global variation in morphogenesis genes indicates an involvement of cell differentiation pathways in the decision phase leading up to the commitment to antibiotic biosynthesis. CONCLUSIONS: Our study provides the first detailed insights into the complex sequence of early regulatory events during and preceding the major metabolic switch in S. coelicolor, which will form the starting point for future attempts at engineering antibiotic production in a biotechnological setting.

Phosphate acquisition components of the Myxococcus xanthus Pho regulon are regulated by both phosphate availability and development.

In many organisms, phosphatase expression and phosphate (P) uptake are coordinately regulated by the Pho regulon. In Myxococcus xanthus P limitation initiates multicellular development, a process associated with changes in phosphatase expression. We sought here to characterize the link between P acquisition and development in this bacterium, an organism capable of preying upon other microorganisms as a sole nutrient source. M. xanthus seems to possess no significant internal P stores, as reducing the P concentration to less than 10 microM retarded growth within one doubling time. Pyrophosphate, polyphosphate, and glyceraldehyde-3-phosphate could support growth as sole P sources, although many other P-containing biomolecules could not (including nucleic acids and phospholipids). Several Pho regulon promoters were found to be highly active during vegetative growth, and P limitation specifically induced pstSCAB, AcPA1, and pho3 promoter activity and repressed pit expression. Enhanced pstSCAB and pho3 promoter activities in a phoP4 mutant (in the presence of high and low concentrations of P) suggested that PhoP4 acts as a repressor of these genes. However, in a phoP4 background, the activities of pstSCAB remained P regulated, suggesting that there is additional regulation by a P-sensitive factor. Initiation of multicellular development caused immediate down-regulation of Pho regulon genes and caused pstSCAB and pho3 promoter activities to become P insensitive. Hence, P acquisition components of the M. xanthus Pho regulon are regulated by both P availability and development, with developmental down-regulation overriding up-regulation by P limitation. These observations suggest that when development is initiated, subsequent changes in P availability become irrelevant to the population, which presumably has sufficient intrinsic P to ensure completion of the developmental program.

Protein-protein interactions between two-component system transmitter and receiver domains of Myxococcus xanthus.

We present a novel dataset assessing the specificity of protein-protein interactions between 69 transmitter and receiver domains from two-component system (TCS)-signalling pathways. TCS require a conserved protein-protein interaction between partner transmitter and receiver domains for signal transduction. The complex prokaryote Myxococcus xanthus possesses an unusually large number of TCS genes, many of which have no obvious interaction partners. Interactions between TCS domains of M. xanthus were assessed using a yeast two-hybrid assay, in which domains were expressed as both bait and prey translational fusions. LacZ production was monitored as an indicator of protein-protein interaction, and the strength of interactions classified as weak, medium or strong. Two-hundred and fifty-five transmitter-receiver domain interactions were observed (46 strong), allowing identification of potential signalling partners for individual M. xanthus TCS proteins. In addition, the dataset provides interesting 'meta' information. For instance, many strong interactions were identified between different transmitter domain pairs (34) and receiver domain pairs (23), suggesting a surprisingly large degree of heterodimerisation of these domains. Proteins in our dataset that exhibited similar 'profiles' of interactions, often shared a similar biological function, suggesting that interaction profiles can provide information on biological function, even considering sets of homologous domains.

The structure/function relationship of a dual-substrate (betaalpha)8-isomerase.

Two structures of phosphoribosyl isomerase A (PriA) from Streptomyces coelicolor, involved in both histidine and tryptophan biosynthesis, were solved at 1.8A resolution. A closed conformer was obtained, which represents the first complete structure of PriA, revealing hitherto unnoticed molecular interactions and the occurrence of conformational changes. Inspection of these conformers, including ligand-docking simulations, allowed identification of residues involved in substrate recognition, chemical catalysis and conformational changes. These predictions were validated by mutagenesis and functional analysis. Arg19 and Ser81 were shown to play critical roles within the carboxyl and amino phosphate-binding sites, respectively; the catalytic residues Asp11 and Asp130 are responsible for both activities; and Thr166 and Asp171, which make an unusual contact, are likely to elicit the conformational changes needed for adopting the active site architectures. This represents the first report of the structure/function relationship of this (betaalpha)8-isomerase.

Reductions in Resting Blood Pressure in Young Adults When Isometric Exercise Is Performed Whilst Walking.

Aerobic and isometric training have been shown to reduce resting blood pressure, but simultaneous aerobic and isometric training have not been studied. The purpose of this study was to compare the changes in resting systolic (SBP), diastolic (DBP), and mean arterial blood pressure (MAP) after 6 weeks of either (i) simultaneous walking and isometric handgrip exercise (WHG), (ii) walking (WLK), (iii) isometric handgrip exercise (IHG), or control (CON). Forty-eight healthy sedentary participants (age 20.7 ± 1.7 yrs, mass 67.2 ± 10.2 kg, height 176.7 ± 1.2 cm, male n = 26, and female n = 22) were randomly allocated, to one of four groups (n = 12 in each). Training was performed 4 × week-1 and involved either treadmill walking for 30 minutes (WLK), handgrip exercise 3 × 10 s at 20% MVC (IHG), or both performed simultaneously (WHG). Resting SBP, DBP, and MAP were recorded at rest, before and after the 6-week study period. Reductions in resting blood pressure were significantly greater in the simultaneous walking and handgrip group than any other group. These results show that simultaneous walking and handgrip training may have summative effects on reductions in resting blood pressure.

Molecular annotation of ketol-acid reductoisomerases from Streptomyces reveals a novel amino acid biosynthesis interlock mediated by enzyme promiscuity.

The 6-phosphogluconate dehydrogenase superfamily oxidize and reduce a wide range of substrates, making their functional annotation challenging. Ketol-acid reductoisomerase (KARI), encoded by the ilvC gene in branched-chain amino acids biosynthesis, is a promiscuous reductase enzyme within this superfamily. Here, we obtain steady-state enzyme kinetic parameters for 10 IlvC homologues from the genera Streptomyces and Corynebacterium, upon eight selected chemically diverse substrates, including some not normally recognized by enzymes of this superfamily. This biochemical data suggested a Streptomyces biosynthetic interlock between proline and the branched-chain amino acids, mediated by enzyme substrate promiscuity, which was confirmed via mutagenesis and complementation analyses of the proC, ilvC1 and ilvC2 genes in Streptomyces coelicolor. Moreover, both ilvC orthologues and paralogues were analysed, such that the relationship between gene duplication and functional diversification could be explored. The KARI paralogues present in S. coelicolor and Streptomyces lividans, despite their conserved high sequence identity (97%), were shown to be more promiscuous, suggesting a recent functional diversification. In contrast, the KARI paralogue from Streptomyces viridifaciens showed selectivity towards the synthesis of valine precursors, explaining its recruitment within the biosynthetic gene cluster of valanimycin. These results allowed us to assess substrate promiscuity indices as a tool to annotate new molecular functions with metabolic implications.

Optimized submerged batch fermentation strategy for systems scale studies of metabolic switching in Streptomyces coelicolor A3(2).

BACKGROUND: Systems biology approaches to study metabolic switching in Streptomyces coelicolor A3(2) depend on cultivation conditions ensuring high reproducibility and distinct phases of culture growth and secondary metabolite production. In addition, biomass concentrations must be sufficiently high to allow for extensive time-series sampling before occurrence of a given nutrient depletion for transition triggering. The present study describes for the first time the development of a dedicated optimized submerged batch fermentation strategy as the basis for highly time-resolved systems biology studies of metabolic switching in S. coelicolor A3(2). RESULTS: By a step-wise approach, cultivation conditions and two fully defined cultivation media were developed and evaluated using strain M145 of S. coelicolor A3(2), providing a high degree of cultivation reproducibility and enabling reliable studies of the effect of phosphate depletion and L-glutamate depletion on the metabolic transition to antibiotic production phase. Interestingly, both of the two carbon sources provided, D-glucose and L-glutamate, were found to be necessary in order to maintain high growth rates and prevent secondary metabolite production before nutrient depletion. Comparative analysis of batch cultivations with (i) both L-glutamate and D-glucose in excess, (ii) L-glutamate depletion and D-glucose in excess, (iii) L-glutamate as the sole source of carbon and (iv) D-glucose as the sole source of carbon, reveal a complex interplay of the two carbon sources in the bacterium's central carbon metabolism. CONCLUSIONS: The present study presents for the first time a dedicated cultivation strategy fulfilling the requirements for systems biology studies of metabolic switching in S. coelicolor A3(2). Key results from labelling and cultivation experiments on either or both of the two carbon sources provided indicate that in the presence of D-glucose, L-glutamate was the preferred carbon source, while D-glucose alone appeared incapable of maintaining culture growth, likely due to a metabolic bottleneck at the oxidation of pyruvate to acetyl-CoA.

Multicopy proC in Streptomyces coelicolor A3(2) elicits a transient production of prodiginines, while proC deletion does not yield a proline auxotroph.

The last step of proline biosynthesis is typically catalysed by the enzyme Δ(1)-pyrroline-5-carboxylate reductase, encoded by the proC gene. Complete genome sequencing of Streptomyces coelicolor, a soil-dwelling Gram-positive bacterium that uses proline as a precursor for synthesis of prodiginine, revealed a single copy of this gene. Unexpectedly, disruption of this proC homologue (Sco3337) in S. coelicolor M145 yielded a prototrophic strain, yet the reductase activity of Sco3337 was confirmed by complementation of an Escherichia coli proC mutant. Multicopy proC within different genetic contexts elicited a transient production of prodiginines, which showed differential production kinetics of the two most common forms of this natural product produced by S. coelicolor, i.e. streptorubin B (cyclic) and undecylprodigiosin (linear). The metabolic and evolutionary implications of these observations are discussed.

An approach to pathway reconstruction using whole genome metabolic models and sensitive sequence searching.

Metabolic models have the potential to impact on genome annotation and on the interpretation of gene expression and other high throughput genome data. The genome of Streptomyces coelicolor genome has been sequenced and some 30% of the open reading frames (ORFs) lack any functional annotation. A recently constructed metabolic network model for S. coelicolor highlights biochemical functions which should exist to make the metabolic model complete and consistent. These include 205 reactions for which no ORF is associated. Here we combine protein functional predictions for the unannotated open reading frames in the genome with \'missing but expected\' functions inferred from the metabolic model. The approach allows function predictions to be evaluated in the context of the biochemical pathway reconstruction, and feed back iteratively into the metabolic model. We describe the approach and discuss a few illustrative examples.

The response regulator PhoP4 is required for late developmental events in Myxococcus xanthus.

Phosphate regulation is complex in the developmental prokaryote Myxococcus xanthus, and requires at least four two-component systems (TCSs). Here, the identification and characterization of a member of one TCS, designated PhoP4, is reported. phoP4 insertion and in-frame deletion strains caused spore viability to be decreased by nearly two orders of magnitude, and reduced all three development-specific phosphatase activities by 80-90 % under phosphate-limiting conditions. Microarray and quantitative PCR analyses demonstrated that PhoP4 is also required for appropriate expression of the predicted pstSCAB-phoU operon of inorganic phosphate assimilation genes. Unlike the case for the other three M. xanthus Pho TCSs, the chromosomal region around phoP4 does not contain a partner histidine kinase gene. Yeast two-hybrid analyses reveal that PhoP4 interacts reciprocally with PhoR2, the histidine kinase of the Pho2 TCS; however, the existence of certain phenotypic differences between phoP4 and phoR2 mutants suggests that PhoP4 interacts with another, as-yet unidentified, histidine kinase.

Occurrence of a putative ancient-like isomerase involved in histidine and tryptophan biosynthesis.

We report the occurrence of an isomerase with a putative (betaalpha)(8)-barrel structure involved in both histidine and trypto-phan biosynthesis in Streptomyces coelicolor A3(2) and Mycobacterium tuberculosis HR37Rv. Deletion of a hisA homologue (SCO2050) putatively encoding N'-[(5'-phosphoribosyl)-formimino]-5 amino-imidazole-4-carboxamide ribonucleotide isomerase from the chromosome of S. coelicolor A3(2) generated a double auxotrophic mutant for histidine and tryptophan. The bifunctional gene SCO2050 and its orthologue Rv1603 from M. tuberculosis complemented both hisA and trpF mutants of Escherichia coli. Expression of the E. coli trpF gene in the S. coelicolor mutant only complemented the tryptophan auxo-trophy, and the hisA gene only complemented the histidine auxotrophy. The discovery of this enzyme, which has a broad-substrate specificity, has implications for the evolution of metabolic pathways and may prove to be important for understanding the evolution of the (betaalpha)(8)-barrels.

Development of a synthetic minimal medium for Listeria monocytogenes.

A defined solid and liquid minimal medium, HTM, which contained methionine and cysteine as the sole amino acids, was developed for Listeria monocytogenes. Complex broth-grown L. monocytogenes had to adapt to HTM by inducing amino acid biosyntheis. HTM is the simplest minimal medium available for growth of L. monocytogenes.

Light-induced carotenogenesis in Myxococcus xanthus: functional characterization of the ECF sigma factor CarQ and antisigma factor CarR.

Illumination of dark-grown Myxococcus xanthus with blue light leads to the induction of carotenoid synthesis. Central to this response is the activation of the light-inducible promoter, PcarQRS, and the transcription of three downstream genes, carQ, carR and carS. Sequence analysis predicted that CarQ is a member of the ECF (extracytoplasmic function) subfamily of RNA polymerase sigma factors, and that CarR is an inner membrane protein. Genetic analysis strongly implied that CarR is an antisigma factor that sequesters CarQ in a transcriptionally inactive complex. Using in vitro transcription run-off assays, we present biochemical evidence that CarQ functions as a bacterial sigma factor and is responsible for transcription initiation at PcarQRS. Similar experiments using the crtI promoter failed to implicate CarQ in direct transcription of the crtI gene. Experiments using the yeast two-hybrid system demonstrated a protein-protein interaction between CarQ and CarR, providing evidence of a CarQ-CarR complex. The yeast two-hybrid system data also indicated that CarR is capable of oligomerization. Fractionation of M. xanthus membranes with the detergent sarkosyl showed that CarR was associated with the inner membrane. Furthermore, CarR was found to be unstable in illuminated stationary phase cells, providing a possible mechanism by which the CarR-CarQ complex is disrupted.

Expression, purification and preliminary crystallographic analysis of phosphoribosyl isomerase (PriA) from Streptomyces coelicolor.

The priA gene encoding the enzyme phosphoribosyl isomerase from Streptomyces coelicolor, a novel bifunctional enzyme involved in both histidine and tryptophan biosynthesis, was heterologously expressed and purified in Escherichia coli as an N-terminal His-tag fusion. The purified recombinant enzyme was crystallized using the hanging-drop method in 1.50 M ammonium sulfate and 100 mM sodium citrate pH 4.8. Crystals were obtained of up to 0.05 x 0.05 x 0.3 mm in size. A full data set to 2 A resolution was collected at the ESRF beamline ID14-1 and space group P3(1,2)21 was assigned, with unit-cell parameters a = 65.1, c = 104.7 A.