Search

Quadratic Phase Coupling (QPC) serves as an essential statistical instrument for evaluating nonlinear synchronization within multivariate time series data, especially in signal processing and neuroscience fields. This study explores the precision of QPC detection using numerical estimates derived from cross-bicoherence and bivariate Granger causality within a straightforward, yet noisy, instantaneous multiplier model. It further assesses the impact of accidental statistically significant bifrequency interactions, introducing new metrics such as the ratio of bispectral quadratic phase coupling and the ratio of bivariate Granger causality quadratic phase coupling. Ratios nearing 1 signify a high degree of accuracy in detecting QPC. The coupling strength between interacting channels is identified as a key element that introduces nonlinearities, influencing the signal-to-noise ratio in the output channel. The model is tested across 59 experimental conditions of simulated recordings, with each condition evaluated against six coupling strength values, covering a wide range of carrier frequencies to examine a broad spectrum of scenarios. The findings demonstrate that the bispectral method outperforms bivariate Granger causality, particularly in identifying specific QPC under conditions of very weak couplings and in the presence of noise. The detection of specific QPC is crucial for neuroscience applications aimed at better understanding the temporal and spatial coordination between different brain regions., (© 2024. The Author(s).)

Lactococcus lactis is used as cell-factory and strain selections are regularly performed to improve production processes. When selection regimes only allow desired phenotypes to survive, for instance by using antibiotics to select for cells that do not grow in a specific condition, the presence of more resistant subpopulations with a wildtype genotype severely slows down the procedure. While the food grade organism L. lactis is not often exposed to antibiotics we characterized its response to ampicillin in more detail, to better understand emerging population heterogeneity and how this might affect strain selection procedures. Using growth-dependent viability assays we identified persister subpopulations in stationary and exponential phase. Growth-independent viability assays revealed a 100 times larger subpopulation that did not grow on plates or in liquid medium, but had an intact membrane and could maintain a pH gradient. Over one third of these cells restored their intracellular pH when we induced a temporary collapse, indicating that this subpopulation was metabolically active and in a viable but non-culturable state. Exposure of L. lactis MG1363 to ampicillin therefore results in a heterogeneous population response with different dormancy states. These dormant cells should be considered in survival-based strain selection procedures.

Fluorescent proteins (FPs) are widely used in many organisms, but are commonly characterised in vitro. However, the in vitro properties may poorly reflect in vivo performance. Therefore, we characterised 27 FPs in vivo using Saccharomyces cerevisiae as model organism. We linked the FPs via a T2A peptide to a control FP, producing equimolar expression of the 2 FPs from 1 plasmid. Using this strategy, we characterised the FPs for brightness, photostability, photochromicity and pH-sensitivity, achieving a comprehensive in vivo characterisation. Many FPs showed different in vivo properties compared to existing in vitro data. Additionally, various FPs were photochromic, which affects readouts due to complex bleaching kinetics. Finally, we codon optimized the best performing FPs for optimal expression in yeast, and found that codon-optimization alters FP characteristics. These FPs improve experimental signal readout, opening new experimental possibilities. Our results may guide future studies in yeast that employ fluorescent proteins.

In glaucoma, although axonal injury drives retinal ganglion cell (RGC) death, little is known about the underlying pathomechanisms. To provide new mechanistic insights and identify new biomarkers, we combined latest non-targeting metabolomics analyses to profile altered metabolites in the mouse whole retina 2, 4, and 7 days after optic nerve crush (NC). Ultra-high-performance liquid chromatography quadrupole time-of-flight mass spectrometry and liquid chromatography Fourier transform mass spectrometry covering wide spectrum of metabolites in combination highlighted 30 metabolites that changed its concentration after NC. The analysis displayed similar changes for purine nucleotide and glutathione as reported previously in another animal model of axonal injury and detected multiple metabolites that increased after the injury. After studying the specificity of the identified metabolites to RGCs in histological sections using imaging mass spectrometry, two metabolites, i.e., L-acetylcarnitine and phosphatidylcholine were increased not only preceding the peak of RGC death in the whole retina but also at the RGC layer (2.3-fold and 1.2-fold, respectively). These phospholipids propose novel mechanisms of RGC death and may serve as early biomarkers of axonal injury. The combinatory metabolomics analyses promise to illuminate pathomechanisms, reveal biomarkers, and allow the discovery of new therapeutic targets of glaucoma.

Many organisms have several similar transporters with different affinities for the same substrate. Typically, high-affinity transporters are expressed when substrate is scarce and low-affinity ones when it is abundant. The benefit of using low instead of high-affinity transporters remains unclear, especially when additional nutrient sensors are present. Here, we investigate two hypotheses. It was previously hypothesized that there is a trade-off between the affinity and the catalytic efficiency of transporters, and we find some but no definitive support for it. Additionally, we propose that for uptake by facilitated diffusion, at saturating substrate concentrations, lowering the affinity enhances the net uptake rate by reducing substrate efflux. As a consequence, there exists an optimal, external-substrate-concentration dependent transporter affinity. A computational model of Saccharomyces cerevisiae glycolysis shows that using the low affinity HXT3 transporter instead of the high affinity HXT6 enhances the steady-state flux by 36%. We tried to test this hypothesis with yeast strains expressing a single glucose transporter modified to have either a high or a low affinity. However, due to the intimate link between glucose perception and metabolism, direct experimental proof for this hypothesis remained inconclusive. Still, our theoretical results provide a novel reason for the presence of low-affinity transport systems.

The inherent stochasticity of molecular reactions prevents us from predicting the exact state of single-cells in a population. However, when a population grows at steady-state, the probability to observe a cell with particular combinations of properties is fixed. Here we validate and exploit existing theory on the statistics of single-cell growth in order to predict the probability of phenotypic characteristics such as cell-cycle times, volumes, accuracy of division and cell-age distributions, using real-time imaging data for Bacillus subtilis and Escherichia coli. Our results show that single-cell growth-statistics can accurately be predicted from a few basic measurements. These equations relate different phenotypic characteristics, and can therefore be used in consistency tests of experimental single-cell growth data and prediction of single-cell statistics. We also exploit these statistical relations in the development of a fast stochastic-simulation algorithm of single-cell growth and protein expression. This algorithm greatly reduces computational burden, by recovering the statistics of growing cell-populations from the simulation of only one of its lineages. Our approach is validated by comparison of simulations and experimental data. This work illustrates a methodology for the prediction, analysis and tests of consistency of single-cell growth and protein expression data from a few basic statistical principles.

Protein expression in a single cell depends on its global physiological state. Moreover, genetically-identical cells exhibit variability (noise) in protein expression, arising from the stochastic nature of biochemical processes, cell growth and division. While it is well understood how cellular growth rate influences mean protein expression, little is known about the relationship between growth rate and noise in protein expression. Here we quantify this relationship in Bacillus subtilis by a novel combination of experiments and theory. We measure the effects of promoter activity and growth rate on the expression of a fluorescent protein in single cells. We disentangle the observed protein expression noise into protein-specific and systemic contributions, using theory and variance decomposition. We find that noise in protein expression depends solely on mean expression levels, regardless of whether expression is set by promoter activity or growth rate, and that noise increases linearly with growth rate. Our results can aid studies of (synthetic) gene circuits of single cells and their condition dependence.

Protein expression is shaped by evolutionary processes that tune microbial fitness. The limited biosynthetic capacity of a cell constrains protein expression and forces the cell to carefully manage its protein economy. In a chemostat, the physiology of the cell feeds back on the growth conditions, hindering intuitive understanding of how changes in protein concentration affect fitness. Here, we aim to provide a theoretical framework that addresses the selective pressures and optimal evolutionary-strategies in the chemostat. We show that the optimal enzyme levels are the result of a trade-off between the cost of their production and the benefit of their catalytic function. We also show that deviations from optimal enzyme levels are directly related to selection coefficients. The maximal fitness strategy for an organism in the chemostat is to express a well-defined metabolic subsystem known as an elementary flux mode. Using a coarse-grained, kinetic model of Saccharomyces cerevisiae's metabolism and growth, we illustrate that the dynamics and outcome of evolution in a chemostat can be very counter-intuitive: Strictly-respiring and strictly-fermenting strains can evolve from a common ancestor. This work provides a theoretical framework that relates a kinetic, mechanistic view on metabolism with cellular physiology and evolutionary dynamics in the chemostat.

By clustering various ion channels and transporters, ankyrin-G (AnkG) configures the membrane-excitation platforms in neurons and cardiomyocytes. AnkG itself localizes to specific areas on the plasma membrane via s-palmitoylation of Cys. However, the structural mechanism by which AnkG anchors to the membrane is not understood. In this study, we solved the crystal structures of the reduced and oxidized forms of the AnkG s-palmitoylation domain and used multiple long-term coarse-grained molecular dynamics simulations to analyze their membrane association. Here we report that the membrane anchoring of AnkG was facilitated by s-palmitoylation, defining a stable binding interface on the lipid membrane, and that AnkG without s-palmitoylation also preferred to stay near the membrane but did not have a unique binding interface. This suggests that AnkG in the juxtamembrane region is primed to accept lipid modification at Cys, and once that happens AnkG constitutes a rigid structural base upon which a membrane-excitation platform can be assembled.

Gene circuits that control metabolism should restore metabolic functions upon environmental changes. Whether gene networks are capable of steering metabolism to optimal states is an open question. Here we present a method to identify such optimal gene networks. We show that metabolic network optimisation over a range of environments results in an input-output relationship for the gene network that guarantees optimal metabolic states. Optimal control is possible if the gene network can achieve this input-output relationship. We illustrate our approach with the best-studied regulatory network in yeast, the galactose network. We find that over the entire range of external galactose concentrations, the regulatory network is able to optimally steer galactose metabolism. Only a few gene network parameters affect this optimal regulation. The other parameters can be tuned independently for optimisation of other functions, such as fast and low-noise gene expression. This study highlights gene network plasticity, evolvability, and modular functionality.