Multifractality in spin glasses.

We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the computation of the singularity spectrum and its corresponding Legendre transform. A comparatively small number of site pairs controls the average correlation that governs the response to a magnetic field. We explain how this scenario of dramatic fluctuations (at length scales smaller than the coherence length) can be reconciled with the smooth, self-averaging behavior that has long been considered to describe spin-glass dynamics.

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 the disordered magnetic system, it might be possible to translate our results to other models of interest for memory storage in materials.

The Mpemba effect in spin glasses is a persistent memory effect.

The Mpemba effect occurs when a hot system cools faster than an initially colder one, when both are refrigerated in the same thermal reservoir. Using the custom-built supercomputer Janus II, we study the Mpemba effect in spin glasses and show that it is a nonequilibrium process, governed by the coherence length ξ of the system. The effect occurs when the bath temperature lies in the glassy phase, but it is not necessary for the thermal protocol to cross the critical temperature. In fact, the Mpemba effect follows from a strong relationship between the internal energy and ξ that turns out to be a sure-tell sign of being in the glassy phase. Thus, the Mpemba effect presents itself as an intriguing avenue for the experimental study of the coherence length in supercooled liquids and other glass formers.

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 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.

A statics-dynamics equivalence through the fluctuation-dissipation ratio provides a window into the spin-glass phase from nonequilibrium measurements.

We have performed a very accurate computation of the nonequilibrium fluctuation-dissipation ratio for the 3D Edwards-Anderson Ising spin glass, by means of large-scale simulations on the special-purpose computers Janus and Janus II. This ratio (computed for finite times on very large, effectively infinite, systems) is compared with the equilibrium probability distribution of the spin overlap for finite sizes. Our main result is a quantitative statics-dynamics dictionary, which could allow the experimental exploration of important features of the spin-glass phase without requiring uncontrollable extrapolations to infinite times or system sizes.

Forgetting Memories and Their Attractiveness.

We study numerically the memory that forgets, introduced in 1986 by Parisi by bounding the synaptic strength, with a mechanism that avoids confusion; allows remembering the pattern learned more recently; and has a physiologically very well-defined meaning. We analyze a number of features of this learning for a finite number of neurons and finite number of patterns. We discuss how the system behaves in the large but finite N limit. We analyze the basin of attraction of the patterns that have been learned, and we show that it is exponentially small in the age of the pattern.

Cross-correlations of 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, behavioral 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 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. Although these groups of states are stable across many decades, a sudden reorganization occurs in the last part of the 20th century. We unambiguously demonstrate that cultural evolution of society can be observed and quantified by looking at cultural traits. We think that this kind of quantitative analysis can be possibly extended to other cultural traits: Although databases covering more than one century (such as the one we used) are rare, the cultural evolution on shorter timescales can be studied due to the fact that many human activities are usually recorded in the present digital era.

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.

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.

Thermodynamic glass transition in a spin glass without time-reversal symmetry.

Spin glasses are a longstanding model for the sluggish dynamics that appear at the glass transition. However, spin glasses differ from structural glasses in a crucial feature: they enjoy a time reversal symmetry. This symmetry can be broken by applying an external magnetic field, but embarrassingly little is known about the critical behavior of a spin glass in a field. In this context, the space dimension is crucial. Simulations are easier to interpret in a large number of dimensions, but one must work below the upper critical dimension (i.e., in d < 6) in order for results to have relevance for experiments. Here we show conclusive evidence for the presence of a phase transition in a four-dimensional spin glass in a field. Two ingredients were crucial for this achievement: massive numerical simulations were carried out on the Janus special-purpose computer, and a new and powerful finite-size scaling method.

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 (i.e., the ceRNA hypothesis) that RNAs can regulate each other in extended cross-talk networks. Such an ability might play a major role in posttranscriptional regulation to shape 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 timescales, are far less clear. We have carried out a dynamical study of the ceRNA hypothesis on a model of posttranscriptional regulation. Sensitivity analysis shows that ceRNA cross-talk is dynamically extended, i.e., it may take place on timescales shorter than those required to achieve stationarity even in cases where no cross-talk occurs in the steady state, and is possibly amplified. In addition, in the case of large, transfection-like perturbations, the system may develop a strongly nonlinear, 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 timescales.

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 posttranscriptional regulation. We have studied a minimal model of posttranscriptional regulation 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 nontrivial correlations among ceRNAs can emerge in experimental readouts due to transcriptional fluctuations even in the absence of miRNA-mediated cross-talk.

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.

Identifying essential genes in Escherichia coli from a metabolic optimization principle.

Understanding the organization of reaction fluxes in cellular metabolism from the stoichiometry and the topology of the underlying biochemical network is a central issue in systems biology. In this task, it is important to devise reasonable approximation schemes that rely on the stoichiometric data only, because full-scale kinetic approaches are computationally affordable only for small networks (e.g., red blood cells, approximately 50 reactions). Methods commonly used are based on finding the stationary flux configurations that satisfy mass-balance conditions for metabolites, often coupling them to local optimization rules (e.g., maximization of biomass production) to reduce the size of the solution space to a single point. Such methods have been widely applied and have proven able to reproduce experimental findings for relatively simple organisms in specific conditions. Here, we define and study a constraint-based model of cellular metabolism where neither mass balance nor flux stationarity are postulated and where the relevant flux configurations optimize the global growth of the system. In the case of Escherichia coli, steady flux states are recovered as solutions, although mass-balance conditions are violated for some metabolites, implying a nonzero net production of the latter. Such solutions furthermore turn out to provide the correct statistics of fluxes for the bacterium E. coli in different environments and compare well with the available experimental evidence on individual fluxes. Conserved metabolic pools play a key role in determining growth rate and flux variability. Finally, we are able to connect phenomenological gene essentiality with "frozen" fluxes (i.e., fluxes with smaller allowed variability) in E. coli metabolism.

Spectral fingerprints of nonequilibrium dynamics: The case of a Brownian gyrator.

The same system can exhibit a completely different dynamical behavior when it evolves in equilibrium conditions or when it is driven out-of-equilibrium by, e.g., connecting some of its components to heat baths kept at different temperatures. Here we concentrate on an analytically solvable and experimentally relevant model of such a system-the so-called Brownian gyrator-a two-dimensional nanomachine that performs a systematic, on average, rotation around the origin under nonequilibrium conditions, while no net rotation takes place under equilibrium ones. On this example, we discuss a question whether it is possible to distinguish between two types of a behavior judging not upon the statistical properties of the trajectories of components but rather upon their respective spectral densities. The latter are widely used to characterize diverse dynamical systems and are routinely calculated from the data using standard built-in packages. From such a perspective, we inquire whether the power spectral densities possess some "fingerprint" properties specific to the behavior in nonequilibrium. We show that indeed one can conclusively distinguish between equilibrium and nonequilibrium dynamics by analyzing the cross-correlations between the spectral densities of both components in the short frequency limit, or from the spectral densities of both components evaluated at zero frequency. Our analytical predictions, corroborated by experimental and numerical results, open a new direction for the analysis of a nonequilibrium dynamics.

Recognition capabilities of a Hopfield model with auxiliary hidden neurons.

We study the recognition capabilities of the Hopfield model with auxiliary hidden layers, which emerge naturally upon a Hubbard-Stratonovich transformation. We show that the recognition capabilities of such a model at zero temperature outperform those of the original Hopfield model, due to a substantial increase of the storage capacity and the lack of a naturally defined basin of attraction. The modified model does not fall abruptly into the regime of complete confusion when memory load exceeds a sharp threshold. This latter circumstance, together with an increase of the storage capacity, renders such a modified Hopfield model a promising candidate for further research, with possible diverse applications.

A yield-cost tradeoff governs Escherichia coli's decision between fermentation and respiration in carbon-limited growth.

Living cells react to changes in growth conditions by re-shaping their proteome. This accounts for different stress-response strategies, both specific (i.e., aimed at increasing the availability of stress-mitigating proteins) and systemic (such as large-scale changes in the use of metabolic pathways aimed at a more efficient exploitation of resources). Proteome re-allocation can, however, imply significant biosynthetic costs. Whether and how such costs impact the growth performance are largely open problems. Focusing on carbon-limited E. coli growth, we integrate genome-scale modeling and proteomic data to address these questions at quantitative level. After deriving a simple formula linking growth rate, carbon intake, and biosynthetic costs, we show that optimal growth results from the tradeoff between yield maximization and protein burden minimization. Empirical data confirm that E. coli growth is indeed close to Pareto-optimal over a broad range of growth rates. Moreover, we establish that, while most of the intaken carbon is diverted into biomass precursors, the efficiency of ATP synthesis is the key driver of the yield-cost tradeoff. These findings provide a quantitative perspective on carbon overflow, the origin of growth laws and the multidimensional optimality of E. coli metabolism.

Finding long cycles in graphs.

We analyze the problem of discovering long cycles inside a graph. We propose and test two algorithms for this task. The first one is based on recent advances in statistical mechanics and relies on a message passing procedure. The second follows a more standard Monte Carlo Markov chain strategy. Special attention is devoted to Hamiltonian cycles of (nonregular) random graphs of minimal connectivity equal to 3.

ceRNA crosstalk stabilizes protein expression and affects the correlation pattern of interacting proteins.

Gene expression is a noisy process and several mechanisms, both transcriptional and post-transcriptional, can stabilize protein levels in cells. Much work has focused on the role of miRNAs, showing in particular that miRNA-mediated regulation can buffer expression noise for lowly expressed genes. Here, using in silico simulations and mathematical modeling, we demonstrate that miRNAs can exert a much broader influence on protein levels by orchestrating competition-induced crosstalk between mRNAs. Most notably, we find that miRNA-mediated cross-talk (i) can stabilize protein levels across the full range of gene expression rates, and (ii) modifies the correlation pattern of co-regulated interacting proteins, changing the sign of correlations from negative to positive. The latter feature may constitute a potentially robust signature of the existence of RNA crosstalk induced by endogenous competition for miRNAs in standard cellular conditions.

Sample-to-sample fluctuations of power spectrum of a random motion in a periodic Sinai model.

The Sinai model of a tracer diffusing in a quenched Brownian potential is a much-studied problem exhibiting a logarithmically slow anomalous diffusion due to the growth of energy barriers with the system size. However, if the potential is random but periodic, the regime of anomalous diffusion crosses over to one of normal diffusion once a tracer has diffused over a few periods of the system. Here we consider a system in which the potential is given by a Brownian bridge on a finite interval (0,L) and then periodically repeated over the whole real line and study the power spectrum S(f) of the diffusive process x(t) in such a potential. We show that for most of realizations of x(t) in a given realization of the potential, the low-frequency behavior is S(f)∼A/f^{2}, i.e., the same as for standard Brownian motion, and the amplitude A is a disorder-dependent random variable with a finite support. Focusing on the statistical properties of this random variable, we determine the moments of A of arbitrary, negative, or positive order k and demonstrate that they exhibit a multifractal dependence on k and a rather unusual dependence on the temperature and on the periodicity L, which are supported by atypical realizations of the periodic disorder. We finally show that the distribution of A has a log-normal left tail and exhibits an essential singularity close to the right edge of the support, which is related to the Lifshitz singularity. Our findings are based both on analytic results and on extensive numerical simulations of the process x(t).