Monomeric α-synuclein activates the plasma membrane calcium pump.

Alpha-synuclein (aSN) is a membrane-associated and intrinsically disordered protein, well known for pathological aggregation in neurodegeneration. However, the physiological function of aSN is disputed. Pull-down experiments have pointed to plasma membrane Ca2+ -ATPase (PMCA) as a potential interaction partner. From proximity ligation assays, we find that aSN and PMCA colocalize at neuronal synapses, and we show that calcium expulsion is activated by aSN and PMCA. We further show that soluble, monomeric aSN activates PMCA at par with calmodulin, but independent of the autoinhibitory domain of PMCA, and highly dependent on acidic phospholipids and membrane-anchoring properties of aSN. On PMCA, the key site is mapped to the acidic lipid-binding site, located within a disordered PMCA-specific loop connecting the cytosolic A domain and transmembrane segment 3. Our studies point toward a novel physiological role of monomeric aSN as a stimulator of calcium clearance in neurons through activation of PMCA.

Data integration strategies for whole-cell modeling.

Data makes the world go round-and high quality data is a prerequisite for precise models, especially for whole-cell models (WCM). Data for WCM must be reusable, contain information about the exact experimental background, and should-in its entirety-cover all relevant processes in the cell. Here, we review basic requirements to data for WCM and strategies how to combine them. As a species-specific resource, we introduce the Yeast Cell Model Data Base (YCMDB) to illustrate requirements and solutions. We discuss recent standards for data as well as for computational models including the modeling process as data to be reported. We outline strategies for constructions of WCM despite their inherent complexity.

Enhancement of essential cofactors for in vivo biocatalysis.

A scarcity of cofactors, necessary metabolites or substrates for in vivo enzymatic reactions, is among the major barriers for product synthesis in metabolically engineered cells. This work compares our recently developed cofactor-boosting strategy, which uses xylose reductase (XR) and lactose to increase the intracellular levels of reduced or oxidized nicotinamide adenine dinucleotide (phosphate) (NAD(P)H), adenosine triphosphate (ATP) and acetyl coenzymeA (acetyl-CoA), with other previously reported methods. We demonstrated that the XR/lactose approach enhances levels of sugar alcohols and sugar phosphates, which leads to elevated levels of crucial cofactors required by specific metabolic pathways. The patterns of cofactor enhancement are not uniform and depend upon the specific pathway components that are overexpressed. We term this model the "user-pool" model. Here, we investigated metabolite alteration in the fatty-alcohol-producing system in the presence of XR/lactose within an early time frame (5 min after the bioconversion started). All metabolite data were analyzed using untargeted metabolomics. We found that the XR/lactose system could improve fatty-alcohol production as early as 5 min after the bioconversion started. The enhancement of key cofactors and intermediates, such as hexitol, NAD(P)H, ATP, 3-phosphoglycerate, acetyl-CoA, 6-phosphogluconate (6-PG) and glutathione, was consistent with those previously reported on a longer time scale (after 1 h). However, measurements performed at the early time reported here showed detectable differences in metabolite enhancement patterns, such as those of ATP, NADPH, acetyl-CoA and glutathione. These data could serve as a basis for future analysis of metabolic flux alteration by the XR/lactose system. Comparative analysis of the cofactor enhancement by XR and other methods suggests that XR/lactose can serve as a simple tool to increase levels of various cofactors for microbial cell factories.

kboolnet: a toolkit for the verification, validation, and visualization of reaction-contingency (rxncon) models.

BACKGROUND: Computational models of cell signaling networks are extremely useful tools for the exploration of underlying system behavior and prediction of response to various perturbations. By representing signaling cascades as executable Boolean networks, the previously developed rxncon ("reaction-contingency") formalism and associated Python package enable accurate and scalable modeling of signal transduction even in large (thousands of components) biological systems. The models are split into reactions, which generate states, and contingencies, that impinge on reactions; this avoids the so-called "combinatorial explosion" of system size. Boolean description of the biological system compensates for the poor availability of kinetic parameters which are necessary for quantitative models. Unfortunately, few tools are available to support rxncon model development, especially for large, intricate systems. RESULTS: We present the kboolnet toolkit ( https://github.com/Kufalab-UCSD/kboolnet , complete documentation at https://github.com/Kufalab-UCSD/kboolnet/wiki ), an R package and a set of scripts that seamlessly integrate with the python-based rxncon software and collectively provide a complete workflow for the verification, validation, and visualization of rxncon models. The verification script VerifyModel.R checks for responsiveness to repeated stimulations as well as consistency of steady state behavior. The validation scripts TruthTable.R, SensitivityAnalysis.R, and ScoreNet.R provide various readouts for the comparison of model predictions to experimental data. In particular, ScoreNet.R compares model predictions to a cloud-stored MIDAS-format experimental database to provide a numerical score for tracking model accuracy. Finally, the visualization scripts allow for graphical representations of model topology and behavior. The entire kboolnet toolkit is cloud-enabled, allowing for easy collaborative development; most scripts also allow for the extraction and analysis of individual user-defined "modules". CONCLUSION: The kboolnet toolkit provides a modular, cloud-enabled workflow for the development of rxncon models, as well as their verification, validation, and visualization. This will enable the creation of larger, more comprehensive, and more rigorous models of cell signaling using the rxncon formalism in the future.

Limitation of sucrose biosynthesis shapes carbon partitioning during plant cold acclimation.

Cold acclimation is a multigenic process by which many plant species increase their freezing tolerance. Stabilization of photosynthesis and carbohydrate metabolism plays a crucial role in cold acclimation. To study regulation of primary and secondary metabolism during cold acclimation of Arabidopsis thaliana, metabolic mutants with deficiencies in either starch or flavonoid metabolism were exposed to 4°C. Photosynthesis was determined together with amounts of carbohydrates, anthocyanins, organic acids and enzyme activities of the central carbohydrate metabolism. Starch deficiency was found to significantly delay soluble sugar accumulation during cold acclimation, while starch overaccumulation did not affect accumulation dynamics but resulted in lower total amounts of \sucrose and glucose. Anthocyanin amounts were lowered in both starch deficient and overaccumulating mutants. Vice versa, flavonoid deficiency did not result in a changed starch amount, which suggested a unidirectional signalling link between starch and flavonoid metabolism. Mathematical modelling of carbon metabolism indicated kinetics of sucrose biosynthesis to be limiting for carbon partitioning in leaf tissue during cold exposure. Together with cold-induced dynamics of citrate, fumarate and malate amounts, this provided evidence for a central role of sucrose phosphate synthase activity in carbon partitioning between biosynthetic and dissimilatory pathways which stabilizes photosynthesis and metabolism at low temperature.

Identification of periodic attractors in Boolean networks using a priori information.

Boolean networks (BNs) have been developed to describe various biological processes, which requires analysis of attractors, the long-term stable states. While many methods have been proposed to detection and enumeration of attractors, there are no methods which have been demonstrated to be theoretically better than the naive method and be practically used for large biological BNs. Here, we present a novel method to calculate attractors based on a priori information, which works much and verifiably faster than the naive method. We apply the method to two BNs which differ in size, modeling formalism, and biological scope. Despite these differences, the method presented here provides a powerful tool for the analysis of both networks. First, our analysis of a BN studying the effect of the microenvironment during angiogenesis shows that the previously defined microenvironments inducing the specialized phalanx behavior in endothelial cells (ECs) additionally induce stalk behavior. We obtain this result from an extended network version which was previously not analyzed. Second, we were able to heuristically detect attractors in a cell cycle control network formalized as a bipartite Boolean model (bBM) with 3158 nodes. These attractors are directly interpretable in terms of genotype-to-phenotype relationships, allowing network validation equivalent to an in silico mutagenesis screen. Our approach contributes to the development of scalable analysis methods required for whole-cell modeling efforts.

Acclimation in plants - the Green Hub consortium.

Acclimation is the capacity to adapt to environmental changes within the lifetime of an individual. This ability allows plants to cope with the continuous variation in ambient conditions to which they are exposed as sessile organisms. Because environmental changes and extremes are becoming even more pronounced due to the current period of climate change, enhancing the efficacy of plant acclimation is a promising strategy for mitigating the consequences of global warming on crop yields. At the cellular level, the chloroplast plays a central role in many acclimation responses, acting both as a sensor of environmental change and as a target of cellular acclimation responses. In this Perspective article, we outline the activities of the Green Hub consortium funded by the German Science Foundation. The main aim of this research collaboration is to understand and strategically modify the cellular networks that mediate plant acclimation to adverse environments, employing Arabidopsis, tobacco (Nicotiana tabacum) and Chlamydomonas as model organisms. These efforts will contribute to 'smart breeding' methods designed to create crop plants with improved acclimation properties. To this end, the model oilseed crop Camelina sativa is being used to test modulators of acclimation for their potential to enhance crop yield under adverse environmental conditions. Here we highlight the current state of research on the role of gene expression, metabolism and signalling in acclimation, with a focus on chloroplast-related processes. In addition, further approaches to uncovering acclimation mechanisms derived from systems and computational biology, as well as adaptive laboratory evolution with photosynthetic microbes, are highlighted.

A yeast cell cycle model integrating stress, signaling, and physiology.

The cell division cycle in eukaryotic cells is a series of highly coordinated molecular interactions that ensure that cell growth, duplication of genetic material, and actual cell division are precisely orchestrated to give rise to two viable progeny cells. Moreover, the cell cycle machinery is responsible for incorporating information about external cues or internal processes that the cell must keep track of to ensure a coordinated, timely progression of all related processes. This is most pronounced in multicellular organisms, but also a cardinal feature in model organisms such as baker's yeast. The complex and integrative behavior is difficult to grasp and requires mathematical modeling to fully understand the quantitative interplay of the single components within the entire system. Here, we present a self-oscillating mathematical model of the yeast cell cycle that comprises all major cyclins and their main regulators. Furthermore, it accounts for the regulation of the cell cycle machinery by a series of external stimuli such as mating pheromones and changes in osmotic pressure or nutrient quality. We demonstrate how the external perturbations modify the dynamics of cell cycle components and how the cell cycle resumes after adaptation to or relief from stress.

Sperm migration in the genital tract-In silico experiments identify key factors for reproductive success.

Sperm migration in the female genital tract controls sperm selection and, therefore, reproductive success as male gametes are conditioned for fertilization while their number is dramatically reduced. Mechanisms underlying sperm migration are mostly unknown, since in vivo investigations are mostly unfeasible for ethical or practical reasons. By presenting a spatio-temporal model of the mammalian female genital tract combined with agent-based description of sperm motion and interaction as well as parameterizing it with bovine data, we offer an alternative possibility for studying sperm migration in silico. The model incorporates genital tract geometry as well as biophysical principles of sperm motion observed in vitro such as positive rheotaxis and thigmotaxis. This model for sperm migration from vagina to oviducts was successfully tested against in vivo data from literature. We found that physical sperm characteristics such as velocity and directional stability as well as sperm-fluid interactions and wall alignment are critical for success, i.e. sperms reaching the oviducts. Therefore, we propose that these identified sperm parameters should be considered in detail for conditioning sperm in artificial selection procedures since the natural processes are normally bypassed in reproductive in vitro technologies. The tremendous impact of mucus flow to support sperm accumulation in the oviduct highlights the importance of a species-specific optimum time window for artificial insemination regarding ovulation. Predictions from our extendable in silico experimental system will improve assisted reproduction in humans, endangered species, and livestock.

Loss of hepatic Mboat7 leads to liver fibrosis.

OBJECTIVE: The rs641738C>T variant located near the membrane-bound O-acyltransferase domain containing 7 (MBOAT7) locus is associated with fibrosis in liver diseases, including non-alcoholic fatty liver disease (NAFLD), alcohol-related liver disease, hepatitis B and C. We aim to understand the mechanism by which the rs641738C>T variant contributes to pathogenesis of NAFLD. DESIGN: Mice with hepatocyte-specific deletion of MBOAT7 (Mboat7Δhep) were generated and livers were characterised by histology, flow cytometry, qPCR, RNA sequencing and lipidomics. We analysed the association of rs641738C>T genotype with liver inflammation and fibrosis in 846 NAFLD patients and obtained genotype-specific liver lipidomes from 280 human biopsies. RESULTS: Allelic imbalance analysis of heterozygous human liver samples pointed to lower expression of the MBOAT7 transcript on the rs641738C>T haplotype. Mboat7Δhep mice showed spontaneous steatosis characterised by increased hepatic cholesterol ester content after 10 weeks. After 6 weeks on a high fat, methionine-low, choline-deficient diet, mice developed increased hepatic fibrosis as measured by picrosirius staining (p<0.05), hydroxyproline content (p<0.05) and transcriptomics, while the inflammatory cell populations and inflammatory mediators were minimally affected. In a human biopsied NAFLD cohort, MBOAT7 rs641738C>T was associated with fibrosis (p=0.004) independent of the presence of histological inflammation. Liver lipidomes of Mboat7Δhep mice and human rs641738TT carriers with fibrosis showed increased total lysophosphatidylinositol levels. The altered lysophosphatidylinositol and phosphatidylinositol subspecies in MBOAT7Δhep livers and human rs641738TT carriers were similar. CONCLUSION: Mboat7 deficiency in mice and human points to an inflammation-independent pathway of liver fibrosis that may be mediated by lipid signalling and a potentially targetable treatment option in NAFLD.

Nonalcoholic fatty liver disease stratification by liver lipidomics.

Nonalcoholic fatty liver disease (NAFLD) is a common metabolic dysfunction leading to hepatic steatosis. However, NAFLD's global impact on the liver lipidome is poorly understood. Using high-resolution shotgun mass spectrometry, we quantified the molar abundance of 316 species from 22 major lipid classes in liver biopsies of 365 patients, including nonsteatotic patients with normal or excessive weight, patients diagnosed with NAFL (nonalcoholic fatty liver) or NASH (nonalcoholic steatohepatitis), and patients bearing common mutations of NAFLD-related protein factors. We confirmed the progressive accumulation of di- and triacylglycerols and cholesteryl esters in the liver of NAFL and NASH patients, while the bulk composition of glycerophospho- and sphingolipids remained unchanged. Further stratification by biclustering analysis identified sphingomyelin species comprising n24:2 fatty acid moieties as membrane lipid markers of NAFLD. Normalized relative abundance of sphingomyelins SM 43:3;2 and SM 43:1;2 containing n24:2 and n24:0 fatty acid moieties, respectively, showed opposite trends during NAFLD progression and distinguished NAFL and NASH lipidomes from the lipidome of nonsteatotic livers. Together with several glycerophospholipids containing a C22:6 fatty acid moiety, these lipids serve as markers of early and advanced stages of NAFL.

Identification of small non-coding RNAs responsive to GUN1 and GUN5 related retrograde signals in Arabidopsis thaliana.

Chloroplast perturbations activate retrograde signalling pathways, causing dynamic changes of gene expression. Besides transcriptional control of gene expression, different classes of small non-coding RNAs (sRNAs) act in gene expression control, but comprehensive analyses regarding their role in retrograde signalling are lacking. We performed sRNA profiling in response to norflurazon (NF), which provokes retrograde signals, in Arabidopsis thaliana wild type (WT) and the two retrograde signalling mutants gun1 and gun5. The RNA samples were also used for mRNA and long non-coding RNA profiling to link altered sRNA levels to changes in the expression of their cognate target RNAs. We identified 122 sRNAs from all known sRNA classes that were responsive to NF in the WT. Strikingly, 142 and 213 sRNAs were found to be differentially regulated in both mutants, indicating a retrograde control of these sRNAs. Concomitant with the changes in sRNA expression, we detected about 1500 differentially expressed mRNAs in the NF-treated WT and around 900 and 1400 mRNAs that were differentially regulated in the gun1 and gun5 mutants, with a high proportion (~30%) of genes encoding plastid proteins. Furthermore, around 20% of predicted miRNA targets code for plastid-localised proteins. Among the sRNA-target pairs, we identified pairs with an anticorrelated expression as well pairs showing other expressional relations, pointing to a role of sRNAs in balancing transcriptional changes upon retrograde signals. Based on the comprehensive changes in sRNA expression, we assume a considerable impact of sRNAs in retrograde-dependent transcriptional changes to adjust plastidic and nuclear gene expression.

Deep Hidden Physics Modeling of Cell Signaling Networks.

According to the WHO, cancer is the second most common cause of death worldwide. The social and economic damage caused by cancer is high and rising. In Europe, the annual direct medical expenses alone amount to more than €129 billion. This results in an urgent need for new and sustainable therapeutics, which has currently not been met by the pharmaceutical industry; only 3.4% of cancer drugs entering Phase I clinical trials get to market. Phosphorylation sites are parts of the core machinery of kinase signaling networks, which are known to be dysfunctional in all types of cancer. Indeed, kinases are the second most common drug target yet. However, these inhibitors block all functions of a protein, and they commonly lead to the development of resistance and increased toxicity. To facilitate global and mechanistic modeling of cancer and clinically relevant cell signaling networks, the community will have to develop sophisticated data-driven deep-learning and mechanistic computational models that generate in silico probabilistic predictions of molecular signaling network rearrangements causally implicated in cancer.

Modeling indicates degradation of mRNA and protein as a potential regulation mechanisms during cold acclimation.

Plants are constantly exposed to temperature fluctuations, which have direct effects on all cellular reactions because temperature influences reaction likelihood and speed. Chloroplasts are crucial to temperature acclimation responses of plants, due to their photosynthetic reactions whose products play a central role in plant metabolism. Consequently, chloroplasts serve as sensors of temperature changes and are simultaneously major targets of temperature acclimation. The core subunits of the complexes involved in the light reactions of photosynthesis are encoded in the chloroplast. As a result, it is assumed that temperature acclimation in plants requires regulatory responses in chloroplast gene expression and protein turnover. We conducted western blot experiments to assess changes in the accumulation of two photosynthetic complexes (PSII, and Cytb6f complex) and the ATP synthase in tobacco plants over two days of acclimation to low temperature. Surprisingly, the concentration of proteins within the chloroplast varied negligibly compared to controls. To explain this observation, we used a simplified Ordinary Differential Equation (ODE) model of transcription, translation, mRNA degradation and protein degradation to explain how the protein concentration can be kept constant. This model takes into account temperature effects on these processes. Through simulations of the ODE model, we show that mRNA and protein degradation are possible targets for control during temperature acclimation. Our model provides a basis for future directions in research and the analysis of future results.

Moonlighting proteins - an approach to systematize the concept.

Moonlighting refers to a protein with at least two unrelated, mechanistically different functions. As a concept, moonlighting describes a large and diverse group of proteins which have been discovered in a multitude of organisms. As of today, a systematized view on these proteins is missing. Here, we propose a classification of moonlighting proteins by two classifiers. We use the function of the protein as a first classifier: activating - activating (Type I), activating - inhibiting (Type II), inhibiting - activating (Type III) and inhibiting - inhibiting (Type IV). To further specify the type of moonlighting protein, we used a second classifier based on the character of the factor that switches the function of the protein: external factor affecting the protein (Type A), change in the first pathway (Type B), change in the second pathway (Type C), equal competition between both pathways (Type D). Using a small two-pathway model we simulated these types of moonlighting proteins to elucidate possible behaviors of the types of moonlighting proteins. We find that, using the results of our simulations, we can classify the behavior of the moonlighting types into Blinker, Splitter andSwitch.

Chemical Reaction Networks Possess Intrinsic, Temperature-Dependent Functionality.

Temperature influences the life of many organisms in various ways. A great number of organisms live under conditions where their ability to adapt to changes in temperature can be vital and largely determines their fitness. Understanding the mechanisms and principles underlying this ability to adapt can be of great advantage, for example, to improve growth conditions for crops and increase their yield. In times of imminent, increasing climate change, this becomes even more important in order to find strategies and help crops cope with these fundamental changes. There is intense research in the field of acclimation that comprises fluctuations of various environmental conditions, but most acclimation research focuses on regulatory effects and the observation of gene expression changes within the examined organism. As thermodynamic effects are a direct consequence of temperature changes, these should necessarily be considered in this field of research but are often neglected. Additionally, compensated effects might be missed even though they are equally important for the organism, since they do not cause observable changes, but rather counteract them. In this work, using a systems biology approach, we demonstrate that even simple network motifs can exhibit temperature-dependent functional features resulting from the interplay of network structure and the distribution of activation energies over the involved reactions. The demonstrated functional features are (i) the reversal of fluxes within a linear pathway, (ii) a thermo-selective branched pathway with different flux modes and (iii) the increased flux towards carbohydrates in a minimal Calvin cycle that was designed to demonstrate temperature compensation within reaction networks. Comparing a system's response to either temperature changes or changes in enzyme activity we also dissect the influence of thermodynamic changes versus genetic regulation. By this, we expand the scope of thermodynamic modelling of biochemical processes by addressing further possibilities and effects, following established mathematical descriptions of biophysical properties.

Automated generation of bacterial resource allocation models.

Resource Balance Analysis (RBA) is a computational method based on resource allocation, which performs accurate quantitative predictions of whole-cell states (i.e. growth rate, metabolic fluxes, abundances of molecular machines including enzymes) across growth conditions. We present an integrated workflow of RBA together with the Python package RBApy. RBApy builds bacterial RBA models from annotated genome-scale metabolic models by adding descriptions of cellular processes relevant for growth and maintenance. The package includes functions for model simulation and calibration and for interfacing to Escher maps and Proteomaps for visualization. We demonstrate that RBApy faithfully reproduces results obtained by a hand-curated and experimentally validated RBA model for Bacillus subtilis. We also present a calibrated RBA model of Escherichia coli generated from scratch, which obtained excellent fits to measured flux values and enzyme abundances. RBApy makes whole-cell modelling accessible for a wide range of bacterial wild-type and engineered strains, as illustrated with a CO2-fixing Escherichia coli strain. AVAILABILITY: RBApy is available at /https://github.com/SysBioInra/RBApy, under the licence GNU GPL version 3, and runs on Linux, Mac and Windows distributions.

Spatial modeling of the membrane-cytosolic interface in protein kinase signal transduction.

The spatial architecture of signaling pathways and the interaction with cell size and morphology are complex, but little understood. With the advances of single cell imaging and single cell biology, it becomes crucial to understand intracellular processes in time and space. Activation of cell surface receptors often triggers a signaling cascade including the activation of membrane-attached and cytosolic signaling components, which eventually transmit the signal to the cell nucleus. Signaling proteins can form steep gradients in the cytosol, which cause strong cell size dependence. We show that the kinetics at the membrane-cytosolic interface and the ratio of cell membrane area to the enclosed cytosolic volume change the behavior of signaling cascades significantly. We suggest an estimate of average concentration for arbitrary cell shapes depending on the cell volume and cell surface area. The normalized variance, known from image analysis, is suggested as an alternative measure to quantify the deviation from the average concentration. A mathematical analysis of signal transduction in time and space is presented, providing analytical solutions for different spatial arrangements of linear signaling cascades. Quantification of signaling time scales reveals that signal propagation is faster at the membrane than at the nucleus, while this time difference decreases with the number of signaling components in the cytosol. Our investigations are complemented by numerical simulations of non-linear cascades with feedback and asymmetric cell shapes. We conclude that intracellular signal propagation is highly dependent on cell geometry and, thereby, conveys information on cell size and shape to the nucleus.

Differential T cell response against BK virus regulatory and structural antigens: A viral dynamics modelling approach.

BK virus (BKV) associated nephropathy affects 1-10% of kidney transplant recipients, leading to graft failure in about 50% of cases. Immune responses against different BKV antigens have been shown to have a prognostic value for disease development. Data currently suggest that the structural antigens and regulatory antigens of BKV might each trigger a different mode of action of the immune response. To study the influence of different modes of action of the cellular immune response on BKV clearance dynamics, we have analysed the kinetics of BKV plasma load and anti-BKV T cell response (Elispot) in six patients with BKV associated nephropathy using ODE modelling. The results show that only a small number of hypotheses on the mode of action are compatible with the empirical data. The hypothesis with the highest empirical support is that structural antigens trigger blocking of virus production from infected cells, whereas regulatory antigens trigger an acceleration of death of infected cells. These differential modes of action could be important for our understanding of BKV resolution, as according to the hypothesis, only regulatory antigens would trigger a fast and continuous clearance of the viral load. Other hypotheses showed a lower degree of empirical support, but could potentially explain the clearing mechanisms of individual patients. Our results highlight the heterogeneity of the dynamics, including the delay between immune response against structural versus regulatory antigens, and its relevance for BKV clearance. Our modelling approach is the first that studies the process of BKV clearance by bringing together viral and immune kinetics and can provide a framework for personalised hypotheses generation on the interrelations between cellular immunity and viral dynamics.

Information Thermodynamics for Time Series of Signal-Response Models.

The entropy production in stochastic dynamical systems is linked to the structure of their causal representation in terms of Bayesian networks. Such a connection was formalized for bipartite (or multipartite) systems with an integral fluctuation theorem in [Phys. Rev. Lett. 111, 180603 (2013)]. Here we introduce the information thermodynamics for time series, that are non-bipartite in general, and we show that the link between irreversibility and information can only result from an incomplete causal representation. In particular, we consider a backward transfer entropy lower bound to the conditional time series irreversibility that is induced by the absence of feedback in signal-response models. We study such a relation in a linear signal-response model providing analytical solutions, and in a nonlinear biological model of receptor-ligand systems where the time series irreversibility measures the signaling efficiency.