Francesca Di Patti

Enhanced stochastic oscillations in autocatalytic reactions

We study a simplified scheme of k coupled autocatalytic reactions, previously introduced by Togashi and Kaneko. The role of stochastic fluctuations is elucidated through the use of the van Kampen system-size expansion and the results compared with direct stochastic simulations. Regular temporal oscillations are predicted to occur for the concentration of the various chemical constituents, with an enhanced amplitude resulting from a resonance which is induced by the intrinsic graininess of the system. The associated power spectra are determined and have a different form depending on the number of chemical constituents k . We make detailed comparisons in the two cases k=4 and k=8 . Agreement between the theoretical and numerical results for the power spectrum is good in both cases. The resulting spectrum is especially interesting in the k=8 system, since it has two peaks, which the system-size expansion is still able to reproduce accurately.

Spatial model of autocatalytic reactions.

Biological cells with all of their surface structure and complex interior stripped away are essentially vesicles--membranes composed of lipid bilayers which form closed sacs. Vesicles are thought to be relevant as models of primitive protocells, and they could have provided the ideal environment for prebiotic reactions to occur. In this paper, we investigate the stochastic dynamics of a set of autocatalytic reactions, within a spatially bounded domain, so as to mimic a primordial cell. The discreteness of the constituents of the autocatalytic reactions gives rise to large sustained oscillations even when the number of constituents is quite large. These oscillations are spatiotemporal in nature, unlike those found in previous studies, which consisted only of temporal oscillations. We speculate that these oscillations may have a role in seeding membrane instabilities which lead to vesicle division. In this way synchronization could be achieved between protocell growth and the reproduction rate of the constituents (the protogenetic material) in simple protocells.

Deterministic and stochastic aspects of VEGF-A production and the cooperative behavior of tumoral cell colony.

A model is proposed to study the process of hypoxia-induced angiogenesis in cancer cells. The model accounts for the role played by the vascular endothelial growth factor (VEGF)-A in regulating the oxygen intake. VEGF-A is dynamically controlled by the HIF-1α concentration. If not degraded, HIF-1α can bind to the subunit termed HIF-1β and so experience translocation to the nucleus, to exert its proper transcriptional activity. The delicate balance between these opposing tendencies translates into the emergence of distinct macroscopic behaviors in terms of the associated molecular concentrations that we here trace back to normoxia, hypoxia and death regimes. These aspects are firstly analyzed with reference to the ideal mean-field scenario. Stochastic fluctuations are also briefly discussed and shown to seed a cooperative interaction among cellular units, competing for the same oxygen reservoir.

Multiple-scale theory of topology-driven patterns on directed networks

Dynamical processes on networks are currently being considered in different domains of cross-disciplinary interest. Reaction-diffusion systems hosted on directed graphs are in particular relevant for their widespread applications, from computer networks to traffic systems. Due to the peculiar spectrum of the discrete Laplacian operator, homogeneous fixed points can turn unstable, on a directed support, because of the topology of the network, a phenomenon which cannot be induced on undirected graphs. A linear analysis can be performed to single out the conditions that underly the instability. The complete characterization of the patterns, which are eventually attained beyond the linear regime of exponential growth, calls instead for a full nonlinear treatment. By performing a multiple time scale perturbative calculation, we here derive an effective equation for the nonlinear evolution of the amplitude of the most unstable mode, close to the threshold of criticality. This is a Stuart-Landau equation the complex coefficients of which appear to depend on the topological features of the embedding directed graph. The theory proves adequate versus simulations, as confirmed by operating with a paradigmatic reaction-diffusion model.

Can a microscopic stochastic model explain the emergence of pain cycles in patients

A stochastic model is introduced here to investigate the molecular mechanisms which trigger the perception of pain. The action of analgesic drug compounds is discussed in a dynamical context, where the competition with inactive species is explicitly accounted for. Finite size effects inevitably perturb the mean-field dynamics: oscillations in the amount of bound receptors are spontaneously manifested, driven by the noise which is intrinsic to the system under scrutiny. These effects are investigated both numerically, via stochastic simulations, and analytically, through a large size expansion. The claim that our findings could provide a consistent interpretative framework for explaining the emergence of cyclic behaviors in response to analgesic treatments is substantiated.

Genome-Wide Analysis of Promoters: Clustering by Alignment and Analysis of Regular Patterns

In this paper we perform a genome-wide analysis of H. sapiens promoters. To this aim, we developed and combined two mathematical methods that allow us to (i) classify promoters into groups characterized by specific global structural features, and (ii) recover, in full generality, any regular sequence in the different classes of promoters. One of the main findings of this analysis is that H. sapiens promoters can be classified into three main groups. Two of them are distinguished by the prevalence of weak or strong nucleotides and are characterized by short compositionally biased sequences, while the most frequent regular sequences in the third group are strongly correlated with transposons. Taking advantage of the generality of these mathematical procedures, we have compared the promoter database of H. sapiens with those of other species. We have found that the above-mentioned features characterize also the evolutionary content appearing in mammalian promoters, at variance with ancestral species in the phylogenetic tree, that exhibit a definitely lower level of differentiation among promoters.

Topological stabilization for synchronized dynamics on networks

Abstract A general scheme is proposed and tested to control the symmetry breaking instability of a homogeneous solution of a spatially extended multispecies model, defined on a network. The inherent discreteness of the space makes it possible to act on the topology of the inter-nodes contacts to achieve the desired degree of stabilization, without altering the dynamical parameters of the model. Both symmetric and asymmetric couplings are considered. In this latter setting the web of contacts is assumed to be balanced, for the homogeneous equilibrium to exist. The performance of the proposed method are assessed, assuming the Complex Ginzburg-Landau equation as a reference model. In this case, the implemented control allows one to stabilize the synchronous limit cycle, hence time-dependent, uniform solution. A system of coupled real Ginzburg-Landau equations is also investigated to obtain the topological stabilization of a homogeneous and constant fixed point.

Optimal search strategies on complex multi-linked networks.

In this paper we consider the problem of optimal search strategies on multi-linked networks, i.e. graphs whose nodes are endowed with several independent sets of links. We focus preliminarily on agents randomly hopping along the links of a graph, with the additional possibility of performing non-local hops to randomly chosen nodes with a given probability. We show that an optimal combination of the two jump rules exists that maximises the efficiency of target search, the optimum reflecting the topology of the network. We then generalize our results to multi-linked networks with an arbitrary number of mutually interfering link sets.

System size expansion for systems with an absorbing state

The well-known van Kampen system size expansion, while of rather general applicability, is shown to fail to reproduce some qualitative features of the time evolution for systems with an absorbing state, apart from a transient initial time interval. We generalize the van Kampen ansatz by introducing a new prescription leading to non-Gaussian fluctuations around the absorbing state. The two expansion predictions are explicitly compared for the infinite range voter model with speciation as a paradigmatic model with an absorbing state. The new expansion, both for a finite size system in the large time limit and at finite time in the large size limit, converges to the exact solution as obtained in a numerical implementation using the Gillespie algorithm. Furthermore, the predicted lifetime distribution is shown to have the correct asymptotic behavior.

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