A Dynamic Genome-Scale Model Identifies Metabolic Pathways Associated with Cold Tolerance in Saccharomyces kudriavzevii

Saccharomyces kudriavzevii is a cold-tolerant species identified as a good alternative for industrial winemaking. Although S. kudriavzevii has never been found in winemaking, its co-occurrence with Saccharomyces cerevisiae in Mediterranean oaks is well documented. This sympatric association is believed to be possible due to the different growth temperatures of the two yeast species. However, the mechanisms behind the cold tolerance of S. kudriavzevii are not well understood. In this work, we propose the use of a dynamic genome-scale model to compare the metabolic routes used by S. kudriavzevii at two temperatures, 25°C and 12°C, to decipher pathways relevant to cold tolerance. The model successfully recovered the dynamics of biomass and external metabolites and allowed us to link the observed phenotype with exact intracellular pathways. The model predicted fluxes that are consistent with previous findings, but it also led to novel results which we further confirmed with intracellular metabolomics and transcriptomic data. The proposed model (along with the corresponding code) provides a comprehensive picture of the mechanisms of cold tolerance that occur within S. kudriavzevii. The proposed strategy offers a systematic approach to explore microbial diversity from extracellular fermentation data at low temperatures. IMPORTANCE Nonconventional yeasts promise to provide new metabolic pathways for producing industrially relevant compounds and tolerating specific stressors such as cold temperatures. The mechanisms behind the cold tolerance of S. kudriavzevii or its sympatric relationship with S. cerevisiae in Mediterranean oaks are not well understood. This study proposes a dynamic genome-scale model to investigate metabolic pathways relevant to cold tolerance. The predictions of the model would indicate the ability of S. kudriavzevii to produce assimilable nitrogen sources from extracellular proteins present in its natural niche. These predictions were further confirmed with metabolomics and transcriptomic data. This finding suggests that not only the different growth temperature preferences but also this proteolytic activity may contribute to the sympatric association with S. cerevisiae. Further exploration of these natural adaptations could lead to novel engineering targets for the biotechnological industry.

Dynamic genome-scale modeling of Saccharomyces cerevisiae unravels mechanisms for ester formation during alcoholic fermentation

Fermentation employing Saccharomyces cerevisiae has produced alcoholic beverages and bread for millennia. More recently, S. cerevisiae has been used to manufacture specific metabolites for the food, pharmaceutical, and cosmetic industries. Among the most important of these metabolites are compounds associated with desirable aromas and flavors, including higher alcohols and esters. Although the physiology of yeast has been well-studied, its metabolic modulation leading to aroma production in relevant industrial scenarios such as winemaking is still unclear. Here we ask what are the underlying metabolic mechanisms that explain the conserved and varying behavior of different yeasts regarding aroma formation under enological conditions? We employed dynamic flux balance analysis (dFBA) to answer this key question using the latest genome-scale metabolic model (GEM) of S. cerevisiae. The model revealed several conserved mechanisms among wine yeasts, for example, acetate ester formation is dependent on intracellular metabolic acetyl-CoA/CoA levels, and the formation of ethyl esters facilitates the removal of toxic fatty acids from cells using CoA. Species-specific mechanisms were also found, such as a preference for the shikimate pathway leading to more 2-phenylethanol production in the Opale strain as well as strain behavior varying notably during the carbohydrate accumulation phase and carbohydrate accumulation inducing redox restrictions during a later cell growth phase for strain Uvaferm. In conclusion, our new metabolic model of yeast under enological conditions revealed key metabolic mechanisms in wine yeasts, which will aid future research strategies to optimize their behavior in industrial settings.

Modelling the physiological status of yeast during wine fermentation enables the prediction of secondary metabolism

Saccharomyces non-cerevisiae yeasts are gaining momentum in wine fermentation due to their potential to reduce ethanol content and achieve attractive aroma profiles. However, the design of the fermentation process for new species requires intensive experimentation. The use of mechanistic models could automate process design, yet to date, most fermentation models have focused on primary metabolism. Therefore, these models do not provide insight into the production of secondary metabolites essential for wine quality, such as aromas. In this work, we formulate a continuous model that accounts for the physiological status of yeast, that is, exponential growth, growth under nitrogen starvation and transition to stationary or decay phases. To do so, we assumed that nitrogen starvation is associated with carbohydrate accumulation and the induction of a set of transcriptional changes associated with the stationary phase. The model accurately described the dynamics of time series data for biomass and primary and secondary metabolites obtained for various yeast species in single culture fermentations. We also used the proposed model to explore different process designs, showing how the addition of nitrogen could affect the aromatic profile of wine. This study underlines the potential of incorporating yeast physiology into batch fermentation modelling and provides a new means of automating process design.

First-principles calculations of bulk, surface and interfacial phases and properties of silicon graphite composites as anode materials for lithium ion batteries

The high energy density offered by silicon along with its mineralogical abundance in the earth's crust, make silicon a very promising material for lithium-ion-battery anodes. Despite these potential advantages, graphitic carbon is still the state of the art due to its high conductivity and structural stability upon electrochemical cycling. Composite materials combine the advantages of silicon and graphitic carbon, making them promising materials for the next generation of anodes. However, successfully implementing them in electric vehicles and electronic devices depends on an understanding of the phase, surface and interface properties related to their performance and lifetime. To this end we employ electronic structure calculations to investigate crystalline silicon-graphite surfaces and grain boundaries exhibiting various orientations and degrees of lithiation. We observe a linear relationship between the mixing enthalpies and volumes of both Li-Si and Li-C systems, which results in an empirical relationship between the voltage and the volume expansion of both anode materials. Assuming thermodynamic equilibrium, we find that the lithiation of graphite only commences after Li_xSi has been lithiated to x = 2.5. Furthermore, we find that lithium ions stabilize silicon surfaces, but are unlikely to adsorb on graphite surfaces. Finally, lithium ions stabilize silicon-graphite interfaces, increasing the likelihood of adhesion as core@shell over yolk@shell configurations with increasing degree of lithiation. These observations explain how lithium might accelerate the crystallization of silicon-graphite composites and the formation of smaller nanoparticles with improved performance.

IT governance enablers in relation to IoT implementation: a systematic literature review

Purpose

The purpose of this paper is to obtain a list of recommendations addressed by the information technology (IT) governance enablers in relation to IoT implementation. The reason behind this it is the lack of information about these instances which could the organizations to be more effective when implementing IoT.

Design/methodology/approach

The objectives will be obtained using the methodology – systematic literature review.

Findings

During the research, a list of recommendations was created on each IT governance enabler in relation to IoT implementation, showing the flaws that exist at the literature level for each enabler.

Originality/value

The state of art of this research is a creation of a list of recommendations according to IT governance enablers to be applied on an IoT implementation.

Saccharomyces cerevisiae and S. kudriavzevii Synthetic Wine Fermentation Performance Dissected by Predictive Modeling

Wineries face unprecedented challenges due to new market demands and climate change effects on wine quality. New yeast starters including non-conventional Saccharomyces species, such as S. kudriavzevii, may contribute to deal with some of these challenges. The design of new fermentations using non-conventional yeasts requires an improved understanding of the physiology and metabolism of these cells. Dynamic modeling brings the potential of exploring the most relevant mechanisms and designing optimal processes more systematically. In this work we explore mechanisms by means of a model selection, reduction and cross-validation pipeline which enables to dissect the most relevant fermentation features for the species under consideration, Saccharomyces cerevisiae T73 and Saccharomyces kudriavzevii CR85. The pipeline involved the comparison of a collection of models which incorporate several alternative mechanisms with emphasis on the inhibitory effects due to temperature and ethanol. We focused on defining a minimal model with the minimum number of parameters, to maximize the identifiability and the quality of cross-validation. The selected model was then used to highlight differences in behavior between species. The analysis of model parameters would indicate that the specific growth rate and the transport of hexoses at initial times are higher for S. cervisiae T73 while S. kudriavzevii CR85 diverts more flux for glycerol production and cellular maintenance. As a result, the fermentations with S. kudriavzevii CR85 are typically slower; produce less ethanol but higher glycerol. Finally, we also explored optimal initial inoculation and process temperature to find the best compromise between final product characteristics and fermentation duration. Results reveal that the production of glycerol is distinctive in S. kudriavzevii CR85, it was not possible to achieve the same production of glycerol with S. cervisiae T73 in any of the conditions tested. This result brings the idea that the optimal design of mixed cultures may have an enormous potential for the improvement of final wine quality.

Genome-wide scans between two honeybee populations reveal putative signatures of human-mediated selection

Human‐mediated selection has left signatures in the genomes of many domesticated animals, including the European dark honeybee, Apis mellifera mellifera, which has been selected by apiculturists for centuries. Using whole‐genome sequence information, we investigated selection signatures in spatially separated honeybee subpopulations (Switzerland, n = 39 and France, n = 17). Three different test statistics were calculated in windows of 2 kb (fixation index, cross‐population extended haplotype homozygosity and cross‐population composite likelihood ratio) and combined into a recently developed composite selection score. Applying a stringent false discovery rate of 0.01, we identified six significant selective sweeps distributed across five chromosomes covering eight genes. These genes are associated with multiple molecular and biological functions, including regulation of transcription, receptor binding and signal transduction. Of particular interest is a selection signature on chromosome 1, which corresponds to the WNT4 gene, the family of which is conserved across the animal kingdom with a variety of functions. In Drosophila melanogaster, WNT4 alleles have been associated with differential wing, cross vein and abdominal phenotypes. Defining phenotypic characteristics of different Apis mellifera ssp., which are typically used as selection criteria, include colour and wing venation pattern. This signal is therefore likely to be a good candidate for human mediated‐selection arising from different applied breeding practices in the two managed populations.

Data-driven reverse engineering of signaling pathways using ensembles of dynamic models

Despite significant efforts and remarkable progress, the inference of signaling networks from experimental data remains very challenging. The problem is particularly difficult when the objective is to obtain a dynamic model capable of predicting the effect of novel perturbations not considered during model training. The problem is ill-posed due to the nonlinear nature of these systems, the fact that only a fraction of the involved proteins and their post-translational modifications can be measured, and limitations on the technologies used for growing cells in vitro, perturbing them, and measuring their variations. As a consequence, there is a pervasive lack of identifiability. To overcome these issues, we present a methodology called SELDOM (enSEmbLe of Dynamic lOgic-based Models), which builds an ensemble of logic-based dynamic models, trains them to experimental data, and combines their individual simulations into an ensemble prediction. It also includes a model reduction step to prune spurious interactions and mitigate overfitting. SELDOM is a data-driven method, in the sense that it does not require any prior knowledge of the system: the interaction networks that act as scaffolds for the dynamic models are inferred from data using mutual information. We have tested SELDOM on a number of experimental and in silico signal transduction case-studies, including the recent HPN-DREAM breast cancer challenge. We found that its performance is highly competitive compared to state-of-the-art methods for the purpose of recovering network topology. More importantly, the utility of SELDOM goes beyond basic network inference (i.e. uncovering static interaction networks): it builds dynamic (based on ordinary differential equation) models, which can be used for mechanistic interpretations and reliable dynamic predictions in new experimental conditions (i.e. not used in the training). For this task, SELDOM's ensemble prediction is not only consistently better than predictions from individual models, but also often outperforms the state of the art represented by the methods used in the HPN-DREAM challenge.

SYNBADm: a tool for optimization-based automated design of synthetic gene circuits

MOTIVATION

The design of de novo circuits with predefined performance specifications is a challenging problem in Synthetic Biology. Computational models and tools have proved to be crucial for a successful wet lab implementation. Natural gene circuits are complex, subject to evolutionary tradeoffs and playing multiple roles. However, most synthetic designs implemented to date are simple and perform a single task. As the field progresses, advanced computational tools are needed in order to handle greater levels of circuit complexity in a more flexible way and considering multiple design criteria.

RESULTS

This works presents SYNBADm (SYNthetic Biology Automated optimal Design in Matlab), a software toolbox for the automatic optimal design of gene circuits with targeted functions from libraries of components. This tool makes use of global optimization to simultaneously search the space of structures and kinetic parameters. SYNBADm can efficiently handle high levels of circuit complexity and can consider multiple design criteria through multiobjective optimization. Further, it provides flexible design capabilities, i.e. the user can define the modeling framework, library of components and target performance function(s).

AVAILABILITY AND IMPLEMENTATION

SYNBADm runs under the popular MATLAB computational environment, and is available under GPLv3 license at https://sites.google.com/site/synbadm CONTACT: ireneotero@iim.csic.es or julio@iim.csic.es.

AMIGO2, a toolbox for dynamic modeling, optimization and control in systems biology

Motivation: Many problems of interest in dynamic modeling and control of biological systems can be posed as non-linear optimization problems subject to algebraic and dynamic constraints. In the context of modeling, this is the case of, e.g. parameter estimation, optimal experimental design and dynamic flux balance analysis. In the context of control, model-based metabolic engineering or drug dose optimization problems can be formulated as (multi-objective) optimal control problems. Finding a solution to those problems is a very challenging task which requires advanced numerical methods.

Results: This work presents the AMIGO2 toolbox: the first multiplatform software tool that automatizes the solution of all those problems, offering a suite of state-of-the-art (multi-objective) global optimizers and advanced simulation approaches.

Availability and Implementation: The toolbox and its documentation are available at: sites.google.com/site/amigo2toolbox.

Contact:ebalsa@iim.csic.es

Supplementary information:Supplementary data are available at Bioinformatics online.

Reverse engineering of logic-based differential equation models using a mixed-integer dynamic optimization approach

MOTIVATION

Systems biology models can be used to test new hypotheses formulated on the basis of previous knowledge or new experimental data, contradictory with a previously existing model. New hypotheses often come in the shape of a set of possible regulatory mechanisms. This search is usually not limited to finding a single regulation link, but rather a combination of links subject to great uncertainty or no information about the kinetic parameters.

RESULTS

In this work, we combine a logic-based formalism, to describe all the possible regulatory structures for a given dynamic model of a pathway, with mixed-integer dynamic optimization (MIDO). This framework aims to simultaneously identify the regulatory structure (represented by binary parameters) and the real-valued parameters that are consistent with the available experimental data, resulting in a logic-based differential equation model. The alternative to this would be to perform real-valued parameter estimation for each possible model structure, which is not tractable for models of the size presented in this work. The performance of the method presented here is illustrated with several case studies: a synthetic pathway problem of signaling regulation, a two-component signal transduction pathway in bacterial homeostasis, and a signaling network in liver cancer cells.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONTACT

julio@iim.csic.es or saezrodriguez@ebi.ac.uk.

MEIGO: an open-source software suite based on metaheuristics for global optimization in systems biology and bioinformatics

Background

Optimization is the key to solving many problems in computational biology. Global optimization methods, which provide a robust methodology, and metaheuristics in particular have proven to be the most efficient methods for many applications. Despite their utility, there is a limited availability of metaheuristic tools.

Results

We present MEIGO, an R and Matlab optimization toolbox (also available in Python via a wrapper of the R version), that implements metaheuristics capable of solving diverse problems arising in systems biology and bioinformatics. The toolbox includes the enhanced scatter search method (eSS) for continuous nonlinear programming (cNLP) and mixed-integer programming (MINLP) problems, and variable neighborhood search (VNS) for Integer Programming (IP) problems. Additionally, the R version includes BayesFit for parameter estimation by Bayesian inference. The eSS and VNS methods can be run on a single-thread or in parallel using a cooperative strategy. The code is supplied under GPLv3 and is available at http://www.iim.csic.es/~gingproc/meigo.html. Documentation and examples are included. The R package has been submitted to BioConductor. We evaluate MEIGO against optimization benchmarks, and illustrate its applicability to a series of case studies in bioinformatics and systems biology where it outperforms other state-of-the-art methods.

Conclusions

MEIGO provides a free, open-source platform for optimization that can be applied to multiple domains of systems biology and bioinformatics. It includes efficient state of the art metaheuristics, and its open and modular structure allows the addition of further methods.

CellNOptR: a flexible toolkit to train protein signaling networks to data using multiple logic formalisms

Background

Cells process signals using complex and dynamic networks. Studying how this is performed in a context and cell type specific way is essential to understand signaling both in physiological and diseased situations. Context-specific medium/high throughput proteomic data measured upon perturbation is now relatively easy to obtain but formalisms that can take advantage of these features to build models of signaling are still comparatively scarce.

Results

Here we present CellNOptR, an open-source R software package for building predictive logic models of signaling networks by training networks derived from prior knowledge to signaling (typically phosphoproteomic) data. CellNOptR features different logic formalisms, from Boolean models to differential equations, in a common framework. These different logic model representations accommodate state and time values with increasing levels of detail. We provide in addition an interface via Cytoscape (CytoCopteR) to facilitate use and integration with Cytoscape network-based capabilities.

Conclusions

Models generated with this pipeline have two key features. First, they are constrained by prior knowledge about the network but trained to data. They are therefore context and cell line specific, which results in enhanced predictive and mechanistic insights. Second, they can be built using different logic formalisms depending on the richness of the available data. Models built with CellNOptR are useful tools to understand how signals are processed by cells and how this is altered in disease. They can be used to predict the effect of perturbations (individual or in combinations), and potentially to engineer therapies that have differential effects/side effects depending on the cell type or context.

State-time spectrum of signal transduction logic models

Despite the current wealth of high-throughput data, our understanding of signal transduction is still incomplete. Mathematical modeling can be a tool to gain an insight into such processes. Detailed biochemical modeling provides deep understanding, but does not scale well above relatively a few proteins. In contrast, logic modeling can be used where the biochemical knowledge of the system is sparse and, because it is parameter free (or, at most, uses relatively a few parameters), it scales well to large networks that can be derived by manual curation or retrieved from public databases. Here, we present an overview of logic modeling formalisms in the context of training logic models to data, and specifically the different approaches to modeling qualitative to quantitative data (state) and dynamics (time) of signal transduction. We use a toy model of signal transduction to illustrate how different logic formalisms (Boolean, fuzzy logic and differential equations) treat state and time. Different formalisms allow for different features of the data to be captured, at the cost of extra requirements in terms of computational power and data quality and quantity. Through this demonstration, the assumptions behind each formalism are discussed, as well as their advantages and disadvantages and possible future developments.

Model human behavior: don't constrain it!

Healthcare processes are complex and highly variable from day to day. Healthcare process execution can be affected by any participant in a process, including clinicians, the patient, and the patient's family, as well as environmental factors such as clinician, staff, facility and equipment availability, and patient clinical status. However, only a few solutions exist that enable computer support for a process to address the full complexity and variability of healthcare processes. We have re-conceptualized workflow and developed an innovative process representation and execution framework based on concepts from software engineering, machine learning, complexity, and database management. This new framework frees processes to track human behavior, thereby releasing us from the constraints of past methods. Our approach is also serving as a new architecture for software systems.

Colorimetric detection of eukaryotic gene expression with DNA-derivatized gold nanoparticles

Thiol-linked DNA-gold nanoparticles were used in a novel colorimetric method to detect the presence of specific mRNA from a total RNA extract of yeast cells. The method allowed detection of expression of the FSY1 gene that encodes a specific fructose/H+ symporter in Saccharomyces bayanus PYCC 4565. FSY1 is strongly expressed when the yeast is grown in fructose as the sole carbon source, while cells cultivated in glucose as the sole carbon source repress gene expression. The presence of FSY1 mRNA is detected based on color change of a sample containing total RNA extracted from the organism and gold nanoparticles derivatized with a 15-mer of complementary single stranded DNA upon addition of NaCl. If FSY1 mRNA is present, the solution remains pink, changing to blue-purple in the absence of FSY1 mRNA. Direct detection of specific expression was possible from only 0.3 microg of unamplified total RNA without any further enhancement. This novel method is inexpensive, very easy to perform as no amplification or signal enhancement steps are necessary and takes less than 15 min to develop after total RNA extraction. No temperature control is necessary and color change can be easily detected visually.

Mass spectrometric and thermodynamic studies reveal the role of water molecules in complexes formed between SH2 domains and tyrosyl phosphopeptides

BACKGROUND

SH2 domains have a fundamental role in signal transduction. These domains interact with proteins containing phosphorylated tyrosine residues and, in doing so, mediate the interactions of proteins involved in tyrosine kinase signalling. The issue of specificity in SH2 domain interactions is therefore of great interest in terms of understanding tyrosine kinase signal-transduction pathways and in the discovery of drugs to inhibit them. Water molecules are found at the interfaces of many complexes, however, to date little attention has been paid to their role in dictating specificity.

RESULTS

Here we use a combination of nanoflow electrospray ionization mass spectrometry (ESI-MS), isothermal titration calorimetry and structural data to investigate the effect of water molecules in complexes formed between the SH2 domain of tyrosine kinase Src and tyrosyl phosphopeptides. Binding studies have been performed using a series of different peptides that were selected to allow changes in the water content at the complex interface and demonstrate changes in specificity. ESI-MS enables quantification of the number of water molecules that interact with a higher affinity than those generally found solvating the biomolecular complex.

CONCLUSIONS

Comparing the interactions of different peptides, we show that an intricate network of water molecules have a key role in dictating specificity. The use of mass spectrometry to quantify tightly bound water molecules may prove of general use in structural biology, where an independent determination of the water molecules associated with a structure would be advantageous. Furthermore, the ability to assess whether given water molecules are important in high-affinity binding could make this method a precious tool in drug design.