Source author record

Enzo Marinari

Enzo Marinari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Supervised perceptron learning vs unsupervised Hebbian unlearning: Approaching optimal memory retrieval in Hopfield-like networks

The Hebbian unlearning algorithm, i.e. an unsupervised local procedure used to improve the retrieval properties in Hopfield-like neural networks, is numerically compared to a supervised algorithm to train a linear symmetric perceptron. We analyze the stability of the stored memories: basins of attraction obtained by the Hebbian unlearning technique are found to be comparable in size to those obtained in the symmetric perceptron, while the two algorithms are found to converge in the same region of Gardner's space of interactions, having followed similar learning paths. A geometric interpretation of Hebbian unlearning is proposed to explain its optimal performances. Because the Hopfield model is also a prototypical model of disordered magnetic system, it might be possible to translate our results to other models of interest for memory storage in materials.

preprint2019arXiv

Competing endogenous RNA crosstalk at system level

microRNAs (miRNAs) regulate gene expression at post-transcriptional level by repressing target RNA molecules. Competition to bind miRNAs tends in turn to correlate their targets, establishing effective RNA-RNA interactions that can influence expression levels, buffer fluctuations and promote signal propagation. Such a potential has been characterized mathematically for small motifs both at steady state and \red{away from stationarity}. Experimental evidence, on the other hand, suggests that competing endogenous RNA (ceRNA) crosstalk is rather weak. Extended miRNA-RNA networks could however favour the integration of many crosstalk interactions, leading to significant large-scale effects in spite of the weakness of individual links. To clarify the extent to which crosstalk is sustained by the miRNA interactome, we have studied its emergent systemic features in silico in large-scale miRNA-RNA network reconstructions. We show that, although generically weak, system-level crosstalk patterns (i) are enhanced by transcriptional heterogeneities, (ii) can achieve high-intensity even for RNAs that are not co-regulated, (iii) are robust to variability in transcription rates, and (iv) are significantly non-local, i.e. correlate weakly with miRNA-RNA interaction parameters. Furthermore, RNA levels are generically more stable when crosstalk is strongest. As some of these features appear to be encoded in the network's topology, crosstalk may functionally be favoured by natural selection. These results suggest that, besides their repressive role, miRNAs mediate a weak but resilient and context-independent network of cross-regulatory interactions that interconnect the transcriptome, stabilize expression levels and support system-level responses.

preprint2019arXiv

Passive advection of fractional Brownian motion by random layered flows

We study statistical properties of the process $Y(t)$ of a passive advection by quenched random layered flows in situations when the inter-layer transfer is governed by a fractional Brownian motion $X(t)$ with the Hurst index $H \in (0,1)$. We show that the disorder-averaged mean-squared displacement of the passive advection grows in the large time $t$ limit in proportion to $t^{2 - H}$, which defines a family of anomalous super-diffusions. We evaluate the disorder-averaged Wigner-Ville spectrum of the advection process $Y(t)$ and demonstrate that it has a rather unusual power-law form $1/f^{3 - H}$ with a characteristic exponent which exceed the value $2$. Our results also suggest that sample-to-sample fluctuations of the spectrum can be very important.

preprint2016arXiv

Constrained Allocation Flux Balance Analysis

New experimental results on bacterial growth inspire a novel top-down approach to study cell metabolism, combining mass balance and proteomic constraints to extend and complement Flux Balance Analysis. We introduce here Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic costs associated to growth are accounted for in an effective way through a single additional genome-wide constraint. Its roots lie in the experimentally observed pattern of proteome allocation for metabolic functions, allowing to bridge regulation and metabolism in a transparent way under the principle of growth-rate maximization. We provide a simple method to solve CAFBA efficiently and propose an "ensemble averaging" procedure to account for unknown protein costs. Applying this approach to modeling E. coli metabolism, we find that, as the growth rate increases, CAFBA solutions cross over from respiratory, growth-yield maximizing states (preferred at slow growth) to fermentative states with carbon overflow (preferred at fast growth). In addition, CAFBA allows for quantitatively accurate predictions on the rate of acetate excretion and growth yield based on only 3 parameters determined by empirical growth laws.

preprint2016arXiv

Noise Processing by MicroRNA-Mediated Circuits: the Incoherent Feed-Forward Loop, Revisited

The intrinsic stochasticity of gene expression is usually mitigated in higher eukaryotes by post-transcriptional regulation channels that stabilise the output layer, most notably protein levels. The discovery of small non-coding RNAs (miRNAs) in specific motifs of the genetic regulatory network has led to identifying noise buffering as the possible key function they exert in regulation. Recent in vitro} and in silico studies have corroborated this hypothesis. It is however also known that miRNA-mediated noise reduction is hampered by transcriptional bursting in simple topologies. Here, using stochastic simulations validated by analytical calculations based on van Kampen's expansion, we revisit the noise-buffering capacity of the miRNA-mediated Incoherent Feed Forward Loop (IFFL), a small module that is widespread in the gene regulatory networks of higher eukaryotes, in order to account for the effects of intermittency in the transcriptional activity of the modulator gene. We show that bursting considerably alters the circuit's ability to control static protein noise. By comparing with other regulatory architectures, we find that direct transcriptional regulation significantly outperforms the IFFL in a broad range of kinetic parameters. This suggests that, under pulsatile inputs, static noise reduction may be less important than dynamical aspects of noise and information processing in characterising the performance of regulatory elements.

preprint2016arXiv

Probing the limits to microRNA-mediated control of gene expression

According to the `ceRNA hypothesis', microRNAs (miRNAs) may act as mediators of an effective positive interaction between long coding or non-coding RNA molecules, carrying significant potential implications for a variety of biological processes. Here, inspired by recent work providing a quantitative description of small regulatory elements as information-conveying channels, we characterize the effectiveness of miRNA-mediated regulation in terms of the optimal information flow achievable between modulator (transcription factors) and target nodes (long RNAs). Our findings show that, while a sufficiently large degree of target derepression is needed to activate miRNA-mediated transmission, (a) in case of differential mechanisms of complex processing and/or transcriptional capabilities, regulation by a post-transcriptional miRNA-channel can outperform that achieved through direct transcriptional control; moreover, (b) in the presence of large populations of weakly interacting miRNA molecules the extra noise coming from titration disappears, allowing the miRNA-channel to process information as effectively as the direct channel. These observations establish the limits of miRNA-mediated post-transcriptional cross-talk and suggest that, besides providing a degree of noise buffering, this type of control may be effectively employed in cells both as a failsafe mechanism and as a preferential fine tuner of gene expression, pointing to the specific situations in which each of these functionalities is maximized.

preprint2015arXiv

Cross correlations of the American baby names

The quantitative description of cultural evolution is a challenging task. The most difficult part of the problem is probably to find the appropriate measurable quantities that can make more quantitative such evasive concepts as, for example, dynamics of cultural movements, behavior patterns and traditions of the people. A strategy to tackle this issue is to observe particular features of human activities, i.e. cultural traits, such as names given to newborns. We study the names of babies born in the United States of America from 1910 to 2012. Our analysis shows that groups of different correlated states naturally emerge in different epochs, and we are able to follow and decrypt their evolution. While these groups of states are stable across many decades, a sudden reorganization occurs in the last part of the twentieth century. We think that this kind of quantitative analysis can be possibly extended to other cultural traits: although databases covering more than one century (as the one we used) are rare, the cultural evolution on shorter time scales can be studied thanks to the fact that many human activities are usually recorded in the present digital era.

preprint2015arXiv

Quantitative constraint-based computational model of tumor-to-stroma coupling via lactate shuttle

Cancer cells utilize large amounts of ATP to sustain growth, relying primarily on non-oxidative, fermentative pathways for its production. In many types of cancers this leads, even in the presence of oxygen, to the secretion of carbon equivalents (usually in the form of lactate) in the cell's surroundings, a feature known as the Warburg effect. While the molecular basis of this phenomenon are still to be elucidated, it is clear that the spilling of energy resources contributes to creating a peculiar microenvironment for tumors, possibly characterized by a degree of toxicity. This suggests that mechanisms for recycling the fermentation products (e.g. a lactate shuttle) may be active, effectively inducing a mutually beneficial metabolic coupling between aberrant and non-aberrant cells. Here we analyze this scenario through a large-scale in silico metabolic model of interacting human cells. By going beyond the cell-autonomous description, we show that elementary physico-chemical constraints indeed favor the establishment of such a coupling under very broad conditions. The characterization we obtained by tuning the aberrant cell's demand for ATP, amino-acids and fatty acids and/or the imbalance in nutrient partitioning provides quantitative support to the idea that synergistic multi-cell effects play a central role in cancer sustainment.

preprint2014arXiv

Spontaneous energy-barrier formation in an entropy-driven glassy dynamics

The description of activated relaxation of glassy systems in the multidimensional configurational space is a long-standing open problem. We develop a phenomenological description of the out-of-equilibrium dynamics of a model with a rough potential energy landscape and we analyse it both numerically and analytically. The model provides an example of dynamics where typical relaxation channels go over finite potential energy barriers despite the presence of less-energy-demanding escaping paths in configurational space; we expect this phenomenon to be also relevant in the thermally activated regime of realistic models of glass-formers. In this case, we found that typical dynamical paths episodically reach an high fixed threshold energy unexpectedly giving rise to a simple thermally activated aging phenomenology. In order to unveil this peculiar aging behavior we introduce a novel description of the dynamics in terms of spontaneously emerging dynamical basins.

preprint2013arXiv

Counting and correcting thermodynamically infeasible flux cycles in genome-scale metabolic networks

Thermodynamics constrains the flow of matter in a reaction network to occur through routes along which the Gibbs energy decreases, implying that viable steady-state flux patterns should be void of closed reaction cycles. Identifying and removing cycles in large reaction networks can unfortunately be a highly challenging task from a computational viewpoint. We propose here a method that accomplishes it by combining a relaxation algorithm and a Monte Carlo procedure to detect loops, with ad hoc rules (discussed in detail) to eliminate them. As test cases, we tackle (a) the problem of identifying infeasible cycles in the E. coli metabolic network and (b) the problem of correcting thermodynamic infeasibilities in the Flux-Balance-Analysis solutions for 15 human cell-type specific metabolic networks. Results for (a) are compared with previous analyses of the same issue, while results for (b) are weighed against alternative methods to retrieve thermodynamically viable flux patterns based on minimizing specific global quantities. Our method on one hand outperforms previous techniques and, on the other, corrects loopy solutions to Flux Balance Analysis. As a byproduct, it also turns out to be able to reveal possible inconsistencies in model reconstructions.

preprint2013arXiv

Energy metabolism and glutamate-glutamine cycle in the brain: a stoichiometric modeling perspective

The energetics of cerebral activity critically relies on the functional and metabolic interactions between neurons and astrocytes. Important open questions include the relation between neuronal versus astrocytic energy demand, glucose uptake and intercellular lactate transfer, as well as their dependence on the level of activity. We have developed a large-scale, constraint-based network model of the metabolic partnership between astrocytes and glutamatergic neurons that allows for a quantitative appraisal of the extent to which stoichiometry alone drives the energetics of the system. We find that the velocity of the glutamate-glutamine cycle ($V_{cyc}$) explains part of the uncoupling between glucose and oxygen utilization at increasing $V_{cyc}$ levels. Thus, we are able to characterize different activation states in terms of the tissue oxygen-glucose index (OGI). Calculations show that glucose is taken up and metabolized according to cellular energy requirements, and that partitioning of the sugar between different cell types is not significantly affected by $V_{cyc}$. Furthermore, both the direction and magnitude of the lactate shuttle between neurons and astrocytes turn out to depend on the relative cell glucose uptake while being roughly independent of $V_{cyc}$. These findings suggest that, in absence of ad hoc activity-related constraints on neuronal and astrocytic metabolism, the glutamate-glutamine cycle does not control the relative energy demand of neurons and astrocytes, and hence their glucose uptake and lactate exchange.

preprint2013arXiv

MicroRNAs as a selective channel of communication between competing RNAs: a steady-state theory

It has recently been suggested that the competition for a finite pool of microRNAs (miRNA) gives rise to effective interactions among their common targets (competing endogenous RNAs or ceRNAs) that could prove to be crucial for post-transcriptional regulation (PTR). We have studied a minimal model of PTR where the emergence and the nature of such interactions can be characterized in detail at steady state. Sensitivity analysis shows that binding free energies and repression mechanisms are the key ingredients for the cross-talk between ceRNAs to arise. Interactions emerge in specific ranges of repression values, can be symmetrical (one ceRNA influences another and vice-versa) or asymmetrical (one ceRNA influences another but not the reverse) and may be highly selective, while possibly limited by noise. In addition, we show that non-trivial correlations among ceRNAs can emerge in experimental readouts due to transcriptional fluctuations even in absence of miRNA-mediated cross-talk.

preprint2013arXiv

RNA-based regulation: dynamics and response to perturbations of competing RNAs

The observation that, through a titration mechanism, microRNAs (miRNAs) can act as mediators of effective interactions among their common targets (competing endogenous RNAs or ceRNAs) has brought forward the idea ('ceRNA hypothesis') that RNAs can regulate each other in extended 'cross-talk' networks. Such an ability might play a major role in post-transcriptional regulation (PTR) in shaping a cell's protein repertoire. Recent work focusing on the emergent properties of the cross-talk networks has emphasized the high flexibility and selectivity that may be achieved at stationarity. On the other hand, dynamical aspects, possibly crucial on the relevant time scales, are far less clear. We have carried out a dynamical study of the ceRNA hypothesis on a model of PTR. Sensitivity analysis shows that ceRNA cross-talk is dynamically extended, i.e. it may take place on time scales shorter than those required to achieve stationairity even in cases where no cross-talk occurs in the steady state, and is possibly amplified. Besides, in case of large, transfection-like perturbations the system may develop strongly non-linear, threshold response. Finally, we show that the ceRNA effect provides a very efficient way for a cell to achieve fast positive shifts in the level of a ceRNA when necessary. These results indicate that competition for miRNAs may indeed provide an elementary mechanism to achieve system-level regulatory effects on the transcriptome over physiologically relevant time scales.

preprint2012arXiv

A scalable algorithm to explore the Gibbs energy landscape of genome-scale metabolic networks

The integration of various types of genomic data into predictive models of biological networks is one of the main challenges currently faced by computational biology. Constraint-based models in particular play a key role in the attempt to obtain a quantitative understanding of cellular metabolism at genome scale. In essence, their goal is to frame the metabolic capabilities of an organism based on minimal assumptions that describe the steady states of the underlying reaction network via suitable stoichiometric constraints, specifically mass balance and energy balance (i.e. thermodynamic feasibility). The implementation of these requirements to generate viable configurations of reaction fluxes and/or to test given flux profiles for thermodynamic feasibility can however prove to be computationally intensive. We propose here a fast and scalable stoichiometry-based method to explore the Gibbs energy landscape of a biochemical network at steady state. The method is applied to the problem of reconstructing the Gibbs energy landscape underlying metabolic activity in the human red blood cell, and to that of identifying and removing thermodynamically infeasible reaction cycles in the Escherichia coli metabolic network (iAF1260). In the former case, we produce consistent predictions for chemical potentials (or log-concentrations) of intracellular metabolites; in the latter, we identify a restricted set of loops (23 in total) in the periplasmic and cytoplasmic core as the origin of thermodynamic infeasibility in a large sample ($10^6$) of flux configurations generated randomly and compatibly with the prior information available on reaction reversibility.

preprint2012arXiv

Von Neumann's growth model: statistical mechanics and biological applications

We review recent work on the statistical mechanics of Von Neumann's growth model and discuss its application to cellular metabolic networks. In this context, we present a detailed analysis of the physiological scenario underlying optimality a la Von Neumann in the metabolism of the bacterium Escherichia coli, showing that optimal solutions are characterized by a considerable microscopic flexibility accompanied by a robust emergent picture for the key physiological functions. This suggests that the ideas behind optimal economic growth in Von Neumann's model can be helpful in uncovering functional organization principles of cell energetics.

preprint2011arXiv

Computing fluxes and chemical potential distributions in biochemical networks: energy balance analysis of the human red blood cell

The analysis of non-equilibrium steady states of biochemical reaction networks relies on finding the configurations of fluxes and chemical potentials satisfying stoichiometric (mass balance) and thermodynamic (energy balance) constraints. Efficient methods to explore such states are crucial to predict reaction directionality, calculate physiologic ranges of variability, estimate correlations, and reconstruct the overall energy balance of the network from the underlying molecular processes. While different techniques for sampling the space generated by mass balance constraints are currently available, thermodynamics is generically harder to incorporate. Here we introduce a method to sample the free energy landscape of a reaction network at steady state. In its most general form, it allows to calculate distributions of fluxes and concentrations starting from trial functions that may contain prior biochemical information. We apply our method to the human red blood cell's metabolic network, whose space of mass-balanced flux states has been sampled extensively in recent years. Specifically, we profile its thermodynamically feasible flux configurations, characterizing in detail how fluctuations of fluxes and potentials are correlated. Based on this, we derive the cell's energy balance in terms of entropy production, chemical work done and thermodynamic efficiency.

preprint2010arXiv

Large random correlations in individual mean field spin glass samples

We argue that complex systems must possess long range correlations and illustrate this idea on the example of the mean field spin glass model. Defined on the complete graph, this model has no genuine concept of distance, but the long range character of correlations is translated into a broad distribution of the spin-spin correlation coefficients for almost all realizations of the random couplings. When we sample the whole phase space we find that this distribution is so broad indeed that at low temperatures it essentially becomes uniform, with all possible correlation values appearing with the same probability. The distribution of correlations inside a single phase space valley is also studied and found to be much narrower.

preprint2009arXiv

Intrinsic limitations of inverse inference in the pairwise Ising spin glass

We analyze the limits inherent to the inverse reconstruction of a pairwise Ising spin glass based on susceptibility propagation. We establish the conditions under which the susceptibility propagation algorithm is able to reconstruct the characteristics of the network given first- and second-order local observables, evaluate eventual errors due to various types of noise in the originally observed data, and discuss the scaling of the problem with the number of degrees of freedom.

preprint2008arXiv

Ranking by Loops: a new approach to categorization

We introduce Loop Ranking, a new ranking measure based on the detection of closed paths, which can be computed in an efficient way. We analyze it with respect to several ranking measures which have been proposed in the past, and are widely used to capture the relative importance of the vertices in complex networks. We argue that Loop Ranking is a very appropriate measure to quantify the role of both vertices and edges in the network traffic.

preprint2007arXiv

Inferring DNA sequences from mechanical unzipping data: the large-bandwidth case

The complementary strands of DNA molecules can be separated when stretched apart by a force; the unzipping signal is correlated to the base content of the sequence but is affected by thermal and instrumental noise. We consider here the ideal case where opening events are known to a very good time resolution (very large bandwidth), and study how the sequence can be reconstructed from the unzipping data. Our approach relies on the use of statistical Bayesian inference and of Viterbi decoding algorithm. Performances are studied numerically on Monte Carlo generated data, and analytically. We show how multiple unzippings of the same molecule may be exploited to improve the quality of the prediction, and calculate analytically the number of required unzippings as a function of the bandwidth, the sequence content, the elasticity parameters of the unzipped strands.

preprint2003arXiv

Localization of Denaturation Bubbles in Random DNA Sequences

We study the thermodynamic and dynamic behaviors of twist-induced denaturation bubbles in a long, stretched random sequence of DNA. The small bubbles associated with weak twist are delocalized. Above a threshold torque, the bubbles of several tens of bases or larger become preferentially localized to \AT-rich segments. In the localized regime, the bubbles exhibit ``aging'' and move around sub-diffusively with continuously varying dynamic exponents. These properties are derived using results of large-deviation theory together with scaling arguments, and are verified by Monte-Carlo simulations.

preprint1997arXiv

Numerical Simulations of Spin Glass Systems

We discuss the status of Monte Carlo simulations of (mainly finite dimensional) spin glass systems. After a short historical note and a brief theoretical introduction we start by discussing the (crucial) 3D case: the warm phase, the critical point and the cold phase, the ultrametric structure and the out of equilibrium dynamics. With the same style we discuss the cases of 4D and 2D. In a few appendices we give some details about the definition of states and about the tempering Monte Carlo approach.

preprint1995arXiv

Replica Theory and Large D Josephson Junction Hypercubic Models

We study the statistical mechanics of a $D$-dimensional array of Josephson junctions in presence of a magnetic field on a lattice of side $2$. In the high temperature region the thermodynamical properties can be computed in the limit $D \to \infty$. A conjectural form of the thermodynamic properties in the low temperature phase is obtained by assuming that they are the same of an appropriate spin glass system, based on quenched disordered couplings. Numerical simulations show that this conjecture is very accurate in one regime of the magnetic field, while it is probably slightly inaccurate in a second regime.

preprint1992arXiv

Simulated Tempering: A New Monte Carlo Scheme

We propose a new global optimization method ({\em Simulated Tempering}) for simulating effectively a system with a rough free energy landscape (i.e. many coexisting states) at finite non-zero temperature. This method is related to simulated annealing, but here the temperature becomes a dynamic variable, and the system is always kept at equilibrium. We analyze the method on the Random Field Ising Model, and we find a dramatic improvement over conventional Metropolis and cluster methods. We analyze and discuss the conditions under which the method has optimal performances.