The Lanthipeptide Synthetase-like Protein CA_C0082 Is an Effector of Agr Quorum Sensing in Clostridium acetobutylicum.

Lanthipeptide synthetases are present in all domains of life. They catalyze a crucial step during lanthipeptide biosynthesis by introducing thioether linkages during posttranslational peptide modification. Lanthipeptides have a wide range of functions, including antimicrobial and morphogenetic activities. Intriguingly, several Clostridium species contain lanthipeptide synthetase-like genes of the class II (lanM) family but lack other components of the lanthipeptide biosynthetic machinery. In all instances, these genes are located immediately downstream of putative agr quorum sensing operons. The physiological role and mode of action of the encoded LanM-like proteins remain uncertain as they lack conserved catalytic residues. Here we show for the industrial organism Clostridium acetobutylicum that the LanM-like protein CA_C0082 is not required for the production of active AgrD-derived signaling peptide but nevertheless acts as an effector of Agr quorum sensing. Expression of CA_C0082 was shown to be controlled by the Agr system and is a prerequisite for granulose (storage polymer) formation. The accumulation of granulose, in turn, was shown to be required for maximal spore formation but also to reduce early solvent formation. CA_C0082 and its putative homologs appear to be closely associated with Agr systems predicted to employ signaling peptides with six-membered ring structures and may represent a new subfamily of LanM-like proteins. This is the first time their contribution to bacterial Agr signaling has been described.

A genome-scale metabolic model of Cupriavidus necator H16 integrated with TraDIS and transcriptomic data reveals metabolic insights for biotechnological applications.

Exploiting biological processes to recycle renewable carbon into high value platform chemicals provides a sustainable and greener alternative to current reliance on petrochemicals. In this regard Cupriavidus necator H16 represents a particularly promising microbial chassis due to its ability to grow on a wide range of low-cost feedstocks, including the waste gas carbon dioxide, whilst also naturally producing large quantities of polyhydroxybutyrate (PHB) during nutrient-limited conditions. Understanding the complex metabolic behaviour of this bacterium is a prerequisite for the design of successful engineering strategies for optimising product yields. We present a genome-scale metabolic model (GSM) of C. necator H16 (denoted iCN1361), which is directly constructed from the BioCyc database to improve the readability and reusability of the model. After the initial automated construction, we have performed extensive curation and both theoretical and experimental validation. By carrying out a genome-wide essentiality screening using a Transposon-directed Insertion site Sequencing (TraDIS) approach, we showed that the model could predict gene knockout phenotypes with a high level of accuracy. Importantly, we indicate how experimental and computational predictions can be used to improve model structure and, thus, model accuracy as well as to evaluate potential false positives identified in the experiments. Finally, by integrating transcriptomics data with iCN1361 we create a condition-specific model, which, importantly, better reflects PHB production in C. necator H16. Observed changes in the omics data and in-silico-estimated alterations in fluxes were then used to predict the regulatory control of key cellular processes. The results presented demonstrate that iCN1361 is a valuable tool for unravelling the system-level metabolic behaviour of C. necator H16 and can provide useful insights for designing metabolic engineering strategies.

Toward Full-Stack In Silico Synthetic Biology: Integrating Model Specification, Simulation, Verification, and Biological Compilation.

We present the Infobiotics Workbench (IBW), a user-friendly, scalable, and integrated computational environment for the computer-aided design of synthetic biological systems. It supports an iterative workflow that begins with specification of the desired synthetic system, followed by simulation and verification of the system in high-performance environments and ending with the eventual compilation of the system specification into suitable genetic constructs. IBW integrates modeling, simulation, verification, and biocompilation features into a single software suite. This integration is achieved through a new domain-specific biological programming language, the Infobiotics Language (IBL), which tightly combines these different aspects of in silico synthetic biology into a full-stack integrated development environment. Unlike existing synthetic biology modeling or specification languages, IBL uniquely blends modeling, verification, and biocompilation statements into a single file. This allows biologists to incorporate design constraints within the specification file rather than using decoupled and independent formalisms for different in silico analyses. This novel approach offers seamless interoperability across different tools as well as compatibility with SBOL and SBML frameworks and removes the burden of doing manual translations for standalone applications. We demonstrate the features, usability, and effectiveness of IBW and IBL using well-established synthetic biological circuits.

Engineering improved ethylene production: Leveraging systems biology and adaptive laboratory evolution.

Ethylene is a small hydrocarbon gas widely used in the chemical industry. Annual worldwide production currently exceeds 150 million tons, producing considerable amounts of CO2 contributing to climate change. The need for a sustainable alternative is therefore imperative. Ethylene is natively produced by several different microorganisms, including Pseudomonas syringae pv. phaseolicola via a process catalyzed by the ethylene-forming enzyme (EFE), subsequent heterologous expression of EFE has led to ethylene production in non-native bacterial hosts including Escherichia coli and cyanobacteria. However, solubility of EFE and substrate availability remain rate-limiting steps in biological ethylene production. We employed a combination of genome-scale metabolic modelling, continuous fermentation, and protein evolution to enable the accelerated development of a high efficiency ethylene producing E. coli strain, yielding a 49-fold increase in production, the most significant improvement reported to date. Furthermore, we have clearly demonstrated that this increased yield resulted from metabolic adaptations that were uniquely linked to EFE (wild type versus mutant). Our findings provide a novel solution to deregulate metabolic bottlenecks in key pathways, which can be readily applied to address other engineering challenges.

Dual Predation by Bacteriophage and Bdellovibrio bacteriovorus Can Eradicate Escherichia coli Prey in Situations where Single Predation Cannot.

Bacteria are preyed upon by diverse microbial predators, including bacteriophage and predatory bacteria, such as Bdellovibrio bacteriovorus While bacteriophage are used as antimicrobial therapies in Eastern Europe and are being applied for compassionate use in the United States, predatory bacteria are only just beginning to reveal their potential therapeutic uses. However, predation by either predator type can falter due to different adaptations arising in the prey bacteria. When testing poultry farm wastewater for novel Bdellovibrio isolates on Escherichia coli prey lawns, individual composite plaques were isolated containing both an RTP (rosette-tailed-phage)-like-phage and a B. bacteriovorus strain and showing central prey lysis and halos of extra lysis. Combining the purified phage with a lab strain of B. bacteriovorus HD100 recapitulated haloed plaques and increased killing of the E. coli prey in liquid culture, showing an effective side-by-side action of these predators compared to their actions alone. Using approximate Bayesian computation to select the best fitting from a variety of different mathematical models demonstrated that the experimental data could be explained only by assuming the existence of three prey phenotypes: (i) sensitive to both predators, (ii) genetically resistant to phage only, and (iii) plastic resistant to B. bacteriovorus only. Although each predator reduces prey availability for the other, high phage numbers did not abolish B. bacteriovorus predation, so both predators are competent to coexist and are causing different selective pressures on the bacterial surface while, in tandem, controlling prey bacterial numbers efficiently. This suggests that combinatorial predator therapy could overcome problems of phage resistance.IMPORTANCE With increasing levels of antibiotic resistance, the development of alternative antibacterial therapies is urgently needed. Two potential alternatives are bacteriophage and predatory bacteria. Bacteriophage therapy has been used, but prey/host specificity and the rapid acquisition of bacterial resistance to bacteriophage are practical considerations. Predatory bacteria are of interest due to their broad Gram-negative bacterial prey range and the lack of simple resistance mechanisms. Here, a bacteriophage and a strain of Bdellovibrio bacteriovorus, preyed side by side on a population of E. coli, causing a significantly greater decrease in prey numbers than either alone. Such combinatorial predator therapy may have greater potential than individual predators since prey surface changes selected for by each predator do not protect prey against the other predator.

Gsmodutils: a python based framework for test-driven genome scale metabolic model development.

MOTIVATION: Genome scale metabolic models (GSMMs) are increasingly important for systems biology and metabolic engineering research as they are capable of simulating complex steady-state behaviour. Constraints based models of this form can include thousands of reactions and metabolites, with many crucial pathways that only become activated in specific simulation settings. However, despite their widespread use, power and the availability of tools to aid with the construction and analysis of large scale models, little methodology is suggested for their continued management. For example, when genome annotations are updated or new understanding regarding behaviour is discovered, models often need to be altered to reflect this. This is quickly becoming an issue for industrial systems and synthetic biotechnology applications, which require good quality reusable models integral to the design, build, test and learn cycle. RESULTS: As part of an ongoing effort to improve genome scale metabolic analysis, we have developed a test-driven development methodology for the continuous integration of validation data from different sources. Contributing to the open source technology based around COBRApy, we have developed the gsmodutils modelling framework placing an emphasis on test-driven design of models through defined test cases. Crucially, different conditions are configurable allowing users to examine how different designs or curation impact a wide range of system behaviours, minimizing error between model versions. AVAILABILITY AND IMPLEMENTATION: The software framework described within this paper is open source and freely available from http://github.com/SBRCNottingham/gsmodutils. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Measuring and modelling the response of Klebsiella pneumoniae KPC prey to Bdellovibrio bacteriovorus predation, in human serum and defined buffer.

In worldwide conditions of increasingly antibiotic-resistant hospital infections, it is important to research alternative therapies. Bdellovibrio bacteriovorus bacteria naturally prey on Gram-negative pathogens, including antibiotic-resistant strains and so B. bacteriovorus have been proposed as "living antibiotics" to combat antimicrobially-resistant pathogens. Predator-prey interactions are complex and can be altered by environmental components. To be effective B. bacteriovorus predation needs to work in human body fluids such as serum where predation dynamics may differ to that studied in laboratory media. Here we combine mathematical modelling and lab experimentation to investigate the predation of an important carbapenem-resistant human pathogen, Klebsiella pneumoniae, by B. bacteriovorus in human serum versus buffer. We show experimentally that B. bacteriovorus is able to reduce prey numbers in each environment, on different timescales. Our mathematical model captures the underlying dynamics of the experimentation, including an initial predation-delay at the predator-prey-serum interface. Our research shows differences between predation in buffer and serum and highlights both the potential and limitations of B. bacteriovorus acting therapeutically against K. pneumoniae in serum, informing future research into the medicinal behaviours and dosing of this living antibacterial.

A new accuracy measure based on bounded relative error for time series forecasting.

Many accuracy measures have been proposed in the past for time series forecasting comparisons. However, many of these measures suffer from one or more issues such as poor resistance to outliers and scale dependence. In this paper, while summarising commonly used accuracy measures, a special review is made on the symmetric mean absolute percentage error. Moreover, a new accuracy measure called the Unscaled Mean Bounded Relative Absolute Error (UMBRAE), which combines the best features of various alternative measures, is proposed to address the common issues of existing measures. A comparative evaluation on the proposed and related measures has been made with both synthetic and real-world data. The results indicate that the proposed measure, with user selectable benchmark, performs as well as or better than other measures on selected criteria. Though it has been commonly accepted that there is no single best accuracy measure, we suggest that UMBRAE could be a good choice to evaluate forecasting methods, especially for cases where measures based on geometric mean of relative errors, such as the geometric mean relative absolute error, are preferred.

Identification of Gene Modules Associated with Low Temperatures Response in Bambara Groundnut by Network-Based Analysis.

Bambara groundnut (Vigna subterranea (L.) Verdc.) is an African legume and is a promising underutilized crop with good seed nutritional values. Low temperature stress in a number of African countries at night, such as Botswana, can effect the growth and development of bambara groundnut, leading to losses in potential crop yield. Therefore, in this study we developed a computational pipeline to identify and analyze the genes and gene modules associated with low temperature stress responses in bambara groundnut using the cross-species microarray technique (as bambara groundnut has no microarray chip) coupled with network-based analysis. Analyses of the bambara groundnut transcriptome using cross-species gene expression data resulted in the identification of 375 and 659 differentially expressed genes (p<0.01) under the sub-optimal (23°C) and very sub-optimal (18°C) temperatures, respectively, of which 110 genes are commonly shared between the two stress conditions. The construction of a Highest Reciprocal Rank-based gene co-expression network, followed by its partition using a Heuristic Cluster Chiseling Algorithm resulted in 6 and 7 gene modules in sub-optimal and very sub-optimal temperature stresses being identified, respectively. Modules of sub-optimal temperature stress are principally enriched with carbohydrate and lipid metabolic processes, while most of the modules of very sub-optimal temperature stress are significantly enriched with responses to stimuli and various metabolic processes. Several transcription factors (from MYB, NAC, WRKY, WHIRLY & GATA classes) that may regulate the downstream genes involved in response to stimulus in order for the plant to withstand very sub-optimal temperature stress were highlighted. The identified gene modules could be useful in breeding for low-temperature stress tolerant bambara groundnut varieties.

The Infobiotics Workbench: an integrated in silico modelling platform for Systems and Synthetic Biology.

SUMMARY: The Infobiotics Workbench is an integrated software suite incorporating model specification, simulation, parameter optimization and model checking for Systems and Synthetic Biology. A modular model specification allows for straightforward creation of large-scale models containing many compartments and reactions. Models are simulated either using stochastic simulation or numerical integration, and visualized in time and space. Model parameters and structure can be optimized with evolutionary algorithms, and model properties calculated using probabilistic model checking. AVAILABILITY: Source code and binaries for Linux, Mac and Windows are available at http://www.infobiotics.org/infobiotics-workbench/; released under the GNU General Public License (GPL) version 3. CONTACT: Natalio.Krasnogor@nottingham.ac.uk.

A thermodynamic switch modulates abscisic acid receptor sensitivity.

Abscisic acid (ABA) is a key hormone regulating plant growth, development and the response to biotic and abiotic stress. ABA binding to pyrabactin resistance (PYR)/PYR1-like (PYL)/Regulatory Component of Abscisic acid Receptor (RCAR) intracellular receptors promotes the formation of stable complexes with certain protein phosphatases type 2C (PP2Cs), leading to the activation of ABA signalling. The PYR/PYL/RCAR family contains 14 genes in Arabidopsis and is currently the largest plant hormone receptor family known; however, it is unclear what functional differentiation exists among receptors. Here, we identify two distinct classes of receptors, dimeric and monomeric, with different intrinsic affinities for ABA and whose differential properties are determined by the oligomeric state of their apo forms. Moreover, we find a residue in PYR1, H60, that is variable between family members and plays a key role in determining oligomeric state. In silico modelling of the ABA activation pathway reveals that monomeric receptors have a competitive advantage for binding to ABA and PP2Cs. This work illustrates how receptor oligomerization can modulate hormonal responses and more generally, the sensitivity of a ligand-dependent signalling system.

Stochastic and deterministic multiscale models for systems biology: an auxin-transport case study.

BACKGROUND: Stochastic and asymptotic methods are powerful tools in developing multiscale systems biology models; however, little has been done in this context to compare the efficacy of these methods. The majority of current systems biology modelling research, including that of auxin transport, uses numerical simulations to study the behaviour of large systems of deterministic ordinary differential equations, with little consideration of alternative modelling frameworks. RESULTS: In this case study, we solve an auxin-transport model using analytical methods, deterministic numerical simulations and stochastic numerical simulations. Although the three approaches in general predict the same behaviour, the approaches provide different information that we use to gain distinct insights into the modelled biological system. We show in particular that the analytical approach readily provides straightforward mathematical expressions for the concentrations and transport speeds, while the stochastic simulations naturally provide information on the variability of the system. CONCLUSIONS: Our study provides a constructive comparison which highlights the advantages and disadvantages of each of the considered modelling approaches. This will prove helpful to researchers when weighing up which modelling approach to select. In addition, the paper goes some way to bridging the gap between these approaches, which in the future we hope will lead to integrative hybrid models.