Source author record

Andrea De Martino

Andrea De Martino 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

20works
8topics
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

20 published item(s)

preprint2022arXiv

Relationship between fitness and heterogeneity in exponentially growing microbial populations

Despite major environmental and genetic differences, microbial metabolic networks are known to generate consistent physiological outcomes across vastly different organisms. This remarkable robustness suggests that, at least in bacteria, metabolic activity may be guided by universal principles. The constrained optimization of evolutionarily-motivated objective functions like the growth rate has emerged as the key theoretical assumption for the study of bacterial metabolism. While conceptually and practically useful in many situations, the idea that certain functions are optimized is hard to validate in data. Moreover, it is not always clear how optimality can be reconciled with the high degree of single-cell variability observed in experiments within microbial populations. To shed light on these issues, we develop an inverse modeling framework that connects the fitness of a population of cells (represented by the mean single-cell growth rate) to the underlying metabolic variability through the Maximum-Entropy inference of the distribution of metabolic phenotypes from data. While no clear objective function emerges, we find that, as the medium gets richer, the fitness and inferred variability for Escherichia coli populations follow and slowly approach the theoretically optimal bound defined by minimal reduction of variability at given fitness. These results suggest that bacterial metabolism may be crucially shaped by a population-level trade-off between growth and heterogeneity.

preprint2022arXiv

Statistical-physics approaches to RNA molecules, families and networks

This contribution focuses on the fascinating RNA molecule, its sequence-dependent folding driven by base-pairing interactions, the interplay between these interactions and natural evolution, and its multiple regulatory roles. The four of us have dug into these topics using the tools and the spirit of the statistical physics of disordered systems, and in particular the concept of a disordered (energy/fitness) landscape. After an introduction to RNA molecules and the perspectives they open not only in evolutionary and synthetic biology but also in medicine, we will introduce the important notions of energy and fitness landscapes for these molecules. In Section III we will review some models and algorithms for RNA sequence-to-secondary-structure mapping. Section IV discusses how the secondary-structure energy landscape can be derived from unzipping data. Section V deals with the inference of RNA structure from evolutionary sequence data sampled in different organisms. This will shift the focus from the `sequence-to-structure' mapping described in Section III to a `sequence-to-function' landscape that can be inferred from laboratory evolutionary data on DNA aptamers. Finally, in Section VI, we shall discuss the rich theoretical picture linking networks of interacting RNA molecules to the organization of robust, systemic regulatory programs. Along this path, we will therefore explore phenomena across multiple scales in space, number of molecules and time, showing how the biological complexity of the RNA world can be captured by the unifying concepts of statistical physics.

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.

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

Growth against entropy in bacterial metabolism: the phenotypic trade-off behind empirical growth rate distributions in E. coli

The solution space of genome-scale models of cellular metabolism provides a map between physically viable flux configurations and cellular metabolic phenotypes described, at the most basic level, by the corresponding growth rates. By sampling the solution space of E. coli's metabolic network, we show that empirical growth rate distributions recently obtained in experiments at single-cell resolution can be explained in terms of a trade-off between the higher fitness of fast-growing phenotypes and the higher entropy of slow-growing ones. Based on this, we propose a minimal model for the evolution of a large bacterial population that captures this trade-off. The scaling relationships observed in experiments encode, in such frameworks, for the same distance from the maximum achievable growth rate, the same degree of growth rate maximization, and/or the same rate of phenotypic change. Being grounded on genome-scale metabolic network reconstructions, these results allow for multiple implications and extensions in spite of the underlying conceptual simplicity.

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

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

Boolean constraint satisfaction problems for reaction networks

We define and study a class of (random) Boolean constraint satisfaction problems representing minimal feasibility constraints for networks of chemical reactions. The constraints we consider encode, respectively, for hard mass-balance conditions (where the consumption and production fluxes of each chemical species are matched) and for soft mass-balance conditions (where a net production of compounds is in principle allowed). We solve these constraint satisfaction problems under the Bethe approximation and derive the corresponding Belief Propagation equations, that involve 8 different messages. The statistical properties of ensembles of random problems are studied via the population dynamics methods. By varying a chemical potential attached to the activity of reactions, we find first order transitions and strong hysteresis, suggesting a non-trivial structure in the space of feasible solutions.

preprint2014arXiv

Inferring metabolic phenotypes from the exometabolome through a thermodynamic variational principle

Networks of biochemical reactions, like cellular metabolic networks, are kept in non-equilibrium steady states by the exchange fluxes connecting them to the environment. In most cases, feasible flux configurations can be derived from minimal mass-balance assumptions upon prescribing in- and out-take fluxes. Here we consider the problem of inferring intracellular flux patterns from extracellular metabolite levels. Resorting to a thermodynamic out of equilibrium variational principle to describe the network at steady state, we show that the switch from fermentative to oxidative phenotypes in cells can be characterized in terms of the glucose, lactate, oxygen and carbon dioxide concentrations. Results obtained for an exactly solvable toy model are fully recovered for a large scale reconstruction of human catabolism. Finally we argue that, in spite of the many approximations involved in the theory, available data for several human cell types are well described by the predicted phenotypic map of the problem.

preprint2013arXiv

A Novel Methodology to Estimate Metabolic Flux Distributions in Constraint-Based Models

Quite generally, constraint-based metabolic flux analysis describes the space of viable flux configurations for a metabolic network as a high-dimensional polytope defined by the linear constraints that enforce the balancing of production and consumption fluxes for each chemical species in the system. In some cases, the complexity of the solution space can be reduced by performing an additional optimization, while in other cases, knowing the range of variability of fluxes over the polytope provides a sufficient characterization of the allowed configurations. There are cases, however, in which the thorough information encoded in the individual distributions of viable fluxes over the polytope is required. Obtaining such distributions is known to be a highly challenging computational task when the dimensionality of the polytope is sufficiently large, and the problem of developing cost-effective {\it ad hoc} algorithms has recently seen a major surge of interest. Here, we propose a method that allows us to perform the required computation heuristically in a time scaling {\it linearly} with the number of reactions in the network, overcoming some limitations of similar techniques employed in recent years. As a case study, we apply it to the analysis of the human red blood cell metabolic network, whose solution space can be sampled by different exact techniques, like Hit-and-Run Monte Carlo (scaling roughly like the third power of the system size). Remarkably accurate estimates for the true distributions of viable reaction fluxes are obtained, suggesting that, although further improvements are desirable, our method enhances our ability to analyze the space of allowed configurations for large biochemical reaction networks.

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.

preprint2013arXiv

Searching for feasible stationary states in reaction networks by solving a Boolean constraint satisfaction problem

We analyze the solutions, on single network instances, of a recently introduced class of constraint-satisfaction problems (CSPs), describing feasible steady states of chemical reaction networks. First, we show that the CSPs generalize the scheme known as Network Expansion, which is recovered in a specific limit. Next, a full statistical mechanics characterization (including the phase diagram and a discussion of physical origin of the phase transitions) for Network Expansion is obtained. Finally, we provide a message-passing algorithm to solve the original CSPs in the most general form.

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

Reaction networks as systems for resource allocation: A variational principle for their non-equilibrium steady states

Within a fully microscopic setting, we derive a variational principle for the non-equilibrium steady states of chemical reaction networks, valid for time-scales over which chemical potentials can be taken to be slowly varying: at stationarity the system minimizes a global function of the reaction fluxes with the form of a Hopfield Hamiltonian with Hebbian couplings, that is explicitly seen to correspond to the rate of decay of entropy production over time. Guided by this analogy, we show that reaction networks can be formally re-cast as systems of interacting reactions that optimize the use of the available compounds by competing for substrates, akin to agents competing for a limited resource in an optimal allocation problem. As an illustration, we analyze the scenario that emerges in two simple cases: that of toy (random) reaction networks and that of a metabolic network model of the human red blood cell.

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.