Exploration and Retrieval of Virus-Related Molecular Data Using ExTaxsI: The Monkeypox Use Case

Retrieval and visualization of biological data are essential for understanding complex systems. With the increasing volume of data generated from high-throughput sequencing technologies, effective and optimized data visualization tools have become indispensable. This is particularly relevant in the COVID-19 postpandemic period, where understanding the diversity and interactions of microbial communities (i.e., viral and bacterial) constitutes an important asset to develop and plan suitable interventions.In this chapter, we show the usage and the potentials of ExTaxsI (Exploring Taxonomy Information) tool to retrieve viral biodiversity data stored in National Center for Biotechnology Information (NCBI) databases and create the related visualization. In addition, by integrating different functions and modules, the tool generates relevant types of visualization plots to facilitate the exploration of microbial biodiversity communities useful to deep dive into ecological and taxonomic relationships among different species and identify potential significant targets.Using the Monkeypox virus as a case study, this work points out significant perspectives on biological data visualization, which can be used to gain insights into the ecology, evolution, and pathogenesis of viruses. Accordingly, we show the potentiality of ExTaxsI to organize and describe the available/downloaded data in an easy, simple, and interpretable way allowing the user to interact dynamically with the visualization plots through specific filters, zoom, and explore functions.

Mobile app to perform anonymized longitudinal studies in the context of COVID-19 adverse drug reaction monitoring, leveraging the citizenship engagement

BACKGROUND

Over the last few years, increasingly studies focused on the development of mobile apps as complementary tools to existing pharmacovigilance traditional surveillance systems for improving and facilitating adverse drug reactions reporting.

OBJECTIVE

In this study, we evaluated the potentiality of a new mobile app (vaxEffect@UniMiB) to perform longitudinal studies while preserving the anonymity of the respondents. We applied it to monitor the adverse drug reactions during COVID-19 vaccination campaign in a sample of Italian population.

METHODS

We administered vaxEffect@UniMiB to a convenience sample of academic subjects vaccinated at Milano-Bicocca University hub for COVID-19 during the Italian national vaccination campaign. vaxEffect@UniMiB was developed for both Android and iOS devices. The mobile app asks users to send their medical history and, upon every vaccine administration, their vaccination data and the adverse reactions that occurred within seven days after the vaccination, allowing the follow of reactions dynamic for each respondent. The app sends data over the web to an application server. The web server, along with receiving all user data, saves them in a SQL database server, and reminds patients to submit vaccine and adverse reactions' data by push notifications sent to the mobile app through Firebase Cloud Messaging. On initial startup of the app, a unique user identifier was generated for each respondent, so that its anonymity is completely ensured, while enabling longitudinal studies.

RESULTS

A total of 3712 people have been vaccinated during the first vaccination wave. A total of 2733 respondents between the ages of 19 and 80, coming from the University of Milano-Bicocca and the Politecnico of Milan, participated in the survey. Overall, we collected the information about vaccination and adverse reactions to the first vaccine dose for 2226 subjects (60.0% of vaccinated), to the second dose for 1610 subjects (43.3%), and, in a non-sponsored fashion, to the third dose for 169 individuals.

CONCLUSIONS

vaxEffect@UniMiB revealed to be the first attempt in performing longitudinal studies to monitor the same subject over time in terms of the reported ADRs after each vaccine administration, while guaranteeing at the same time complete anonymity of the subjects. A series of aspects contributed to a positive involvement from people in using this application to report their ADRs to vaccination: ease of use, availability from multiple platforms, anonymity of all the survey participants and protection of the submitted data and the healthcare workers' support.

CLINICALTRIAL

INTEGRATE: Model-based multi-omics data integration to characterize multi-level metabolic regulation

Metabolism is directly and indirectly fine-tuned by a complex web of interacting regulatory mechanisms that fall into two major classes. On the one hand, the expression level of the catalyzing enzyme sets the maximal theoretical flux level (i.e., the net rate of the reaction) for each enzyme-controlled reaction. On the other hand, metabolic regulation controls the metabolic flux through the interactions of metabolites (substrates, cofactors, allosteric modulators) with the responsible enzyme. High-throughput data, such as metabolomics and transcriptomics data, if analyzed separately, do not accurately characterize the hierarchical regulation of metabolism outlined above. They must be integrated to disassemble the interdependence between different regulatory layers controlling metabolism. To this aim, we propose INTEGRATE, a computational pipeline that integrates metabolomics and transcriptomics data, using constraint-based stoichiometric metabolic models as a scaffold. We compute differential reaction expression from transcriptomics data and use constraint-based modeling to predict if the differential expression of metabolic enzymes directly originates differences in metabolic fluxes. In parallel, we use metabolomics to predict how differences in substrate availability translate into differences in metabolic fluxes. We discriminate fluxes regulated at the metabolic and/or gene expression level by intersecting these two output datasets. We demonstrate the pipeline using a set of immortalized normal and cancer breast cell lines. In a clinical setting, knowing the regulatory level at which a given metabolic reaction is controlled will be valuable to inform targeted, truly personalized therapies in cancer patients.

ExTaxsI: an exploration tool of biodiversity molecular data

BACKGROUND

The increasing availability of multi-omics data is leading to regularly revised estimates of existing biodiversity data. In particular, the molecular data enable novel species to be characterized and the information linked to those already observed to be increased with new genomics data. For this reason, the management and visualization of existing molecular data, and their related metadata, through the implementation of easy-to-use IT tools have become a key point to design future research. The more users are able to access biodiversity-related information, the greater the ability of the scientific community to expand its knowledge in this area.

RESULTS

In this article we focus on the development of ExTaxsI (Exploring Taxonomy Information), an IT tool that can retrieve biodiversity data stored in NCBI databases and provide a simple and explorable visualization. We use 3 case studies to show how an efficient organization of the available data can lead to obtaining new information that is fundamental as a starting point for new research. Using this approach highlights the limits in the distribution of data availability, a key factor to consider in the experimental design phase of broad-spectrum studies such as metagenomics.

CONCLUSIONS

ExTaxsI can easily retrieve molecular data and its metadata with an explorable visualization, with the aim of helping researchers to improve experimental designs and highlight the main gaps in the coverage of available data.

Accelerated global sensitivity analysis of genome-wide constraint-based metabolic models

Background

Genome-wide reconstructions of metabolism opened the way to thorough investigations of cell metabolism for health care and industrial purposes. However, the predictions offered by Flux Balance Analysis (FBA) can be strongly affected by the choice of flux boundaries, with particular regard to the flux of reactions that sink nutrients into the system. To mitigate possible errors introduced by a poor selection of such boundaries, a rational approach suggests to focus the modeling efforts on the pivotal ones.

Methods

In this work, we present a methodology for the automatic identification of the key fluxes in genome-wide constraint-based models, by means of variance-based sensitivity analysis. The goal is to identify the parameters for which a small perturbation entails a large variation of the model outcomes, also referred to as sensitive parameters. Due to the high number of FBA simulations that are necessary to assess sensitivity coefficients on genome-wide models, our method exploits a master-slave methodology that distributes the computation on massively multi-core architectures. We performed the following steps: (1) we determined the putative parameterizations of the genome-wide metabolic constraint-based model, using Saltelli’s method; (2) we applied FBA to each parameterized model, distributing the massive amount of calculations over multiple nodes by means of MPI; (3) we then recollected and exploited the results of all FBA runs to assess a global sensitivity analysis.

Results

We show a proof-of-concept of our approach on latest genome-wide reconstructions of human metabolism Recon2.2 and Recon3D. We report that most sensitive parameters are mainly associated with the intake of essential amino acids in Recon2.2, whereas in Recon 3D they are associated largely with phospholipids. We also illustrate that in most cases there is a significant contribution of higher order effects.

Conclusion

Our results indicate that interaction effects between different model parameters exist, which should be taken into account especially at the stage of calibration of genome-wide models, supporting the importance of a global strategy of sensitivity analysis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04002-0.

MaREA4Galaxy: Metabolic reaction enrichment analysis and visualization of RNA-seq data within Galaxy

We present MaREA4Galaxy, a user-friendly tool that allows a user to characterize and to graphically compare groups of samples with different transcriptional regulation of metabolism, as estimated from cross-sectional RNA-seq data. The tool is available as plug-in for the widely-used Galaxy platform for comparative genomics and bioinformatics analyses. MaREA4Galaxy combines three modules. The Expression2RAS module, which, for each reaction of a specified set, computes a Reaction Activity Score (RAS) as a function of the expression level of genes encoding for the associated enzyme. The MaREA (Metabolic Reaction Enrichment Analysis) module that allows to highlight significant differences in reaction activities between specified groups of samples. The Clustering module which employs the RAS computed before as a metric for unsupervised clustering of samples into distinct metabolic subgroups; the Clustering tool provides different clustering techniques and implements standard methods to evaluate the goodness of the results.

Integration of single-cell RNA-seq data into population models to characterize cancer metabolism

Metabolic reprogramming is a general feature of cancer cells. Regrettably, the comprehensive quantification of metabolites in biological specimens does not promptly translate into knowledge on the utilization of metabolic pathways. By estimating fluxes across metabolic pathways, computational models hold the promise to bridge this gap between data and biological functionality. These models currently portray the average behavior of cell populations however, masking the inherent heterogeneity that is part and parcel of tumorigenesis as much as drug resistance. To remove this limitation, we propose single-cell Flux Balance Analysis (scFBA) as a computational framework to translate single-cell transcriptomes into single-cell fluxomes. We show that the integration of single-cell RNA-seq profiles of cells derived from lung adenocarcinoma and breast cancer patients into a multi-scale stoichiometric model of a cancer cell population: significantly 1) reduces the space of feasible single-cell fluxomes; 2) allows to identify clusters of cells with different growth rates within the population; 3) points out the possible metabolic interactions among cells via exchange of metabolites. The scFBA suite of MATLAB functions is available at https://github.com/BIMIB-DISCo/scFBA, as well as the case study datasets.

Cytotoxicity of chemotherapeutic agents and resistance to targeted treatments are the main reasons why cancer is still one of the top causes of death. As tumor cells are intrinsically resistant to therapies that target signaling pathways, targeting the metabolic hallmarks of cancer holds promise for more incisive treatments. Regrettably, the heterogeneity of cancer metabolism hinders the identification of effective treatments. To fully uncover the metabolic heterogeneity within tumors, characterization of metabolic programs (metabolic flux distributions) at the single-cell level is required. To fill the gap between current technologies for genomics and future technologies for fluxomics, both at the single-cell and the genome-wide scale, we propose to integrate cancer data from: 1) single-cell transcriptomics and 2) bulk metabolomics, into a multi-scale stoichiometric model, to deliver for the first time metabolic fluxomes at the single-cell level. To this end, we introduce a new paradigm for flux balance analysis and data integration in cancer metabolism to: 1) characterize metabolic heterogeneity, not only at the inter-, but also at the intra-tumor level 2) identify the metabolic interactions between cancer populations, whose role in resistance to metabolic treatments has been recently recognized 3) predict the collective response to drug targeting of metabolism.

Emerging ensembles of kinetic parameters to characterize observed metabolic phenotypes

Background

Determining the value of kinetic constants for a metabolic system in the exact physiological conditions is an extremely hard task. However, this kind of information is of pivotal relevance to effectively simulate a biological phenomenon as complex as metabolism.

Results

To overcome this issue, we propose to investigate emerging properties of ensembles of sets of kinetic constants leading to the biological readout observed in different experimental conditions. To this aim, we exploit information retrievable from constraint-based analyses (i.e. metabolic flux distributions at steady state) with the goal to generate feasible values for kinetic constants exploiting the mass action law. The sets retrieved from the previous step will be used to parametrize a mechanistic model whose simulation will be performed to reconstruct the dynamics of the system (until reaching the metabolic steady state) for each experimental condition. Every parametrization that is in accordance with the expected metabolic phenotype is collected in an ensemble whose features are analyzed to determine the emergence of properties of a phenotype. In this work we apply the proposed approach to identify ensembles of kinetic parameters for five metabolic phenotypes of E. Coli, by analyzing five different experimental conditions associated with the ECC2comp model recently published by Hädicke and collaborators.

Conclusions

Our results suggest that the parameter values of just few reactions are responsible for the emergence of a metabolic phenotype. Notably, in contrast with constraint-based approaches such as Flux Balance Analysis, the methodology used in this paper does not require to assume that metabolism is optimizing towards a specific goal.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2181-7) contains supplementary material, which is available to authorized users.

popFBA: tackling intratumour heterogeneity with Flux Balance Analysis

Motivation

Intratumour heterogeneity poses many challenges to the treatment of cancer. Unfortunately, the transcriptional and metabolic information retrieved by currently available computational and experimental techniques portrays the average behaviour of intermixed and heterogeneous cell subpopulations within a given tumour. Emerging single-cell genomic analyses are nonetheless unable to characterize the interactions among cancer subpopulations. In this study, we propose popFBA, an extension to classic Flux Balance Analysis, to explore how metabolic heterogeneity and cooperation phenomena affect the overall growth of cancer cell populations.

Results

We show how clones of a metabolic network of human central carbon metabolism, sharing the same stoichiometry and capacity constraints, may follow several different metabolic paths and cooperate to maximize the growth of the total population. We also introduce a method to explore the space of possible interactions, given some constraints on plasma supply of nutrients. We illustrate how alternative nutrients in plasma supply and/or a dishomogeneous distribution of oxygen provision may affect the landscape of heterogeneous phenotypes. We finally provide a technique to identify the most proliferative cells within the heterogeneous population.

Availability and implementation

the popFBA MATLAB function and the SBML model are available at https://github.com/BIMIB-DISCo/popFBA.

A metabolic core model elucidates how enhanced utilization of glucose and glutamine, with enhanced glutamine-dependent lactate production, promotes cancer cell growth: The WarburQ effect

Cancer cells share several metabolic traits, including aerobic production of lactate from glucose (Warburg effect), extensive glutamine utilization and impaired mitochondrial electron flow. It is still unclear how these metabolic rearrangements, which may involve different molecular events in different cells, contribute to a selective advantage for cancer cell proliferation. To ascertain which metabolic pathways are used to convert glucose and glutamine to balanced energy and biomass production, we performed systematic constraint-based simulations of a model of human central metabolism. Sampling of the feasible flux space allowed us to obtain a large number of randomly mutated cells simulated at different glutamine and glucose uptake rates. We observed that, in the limited subset of proliferating cells, most displayed fermentation of glucose to lactate in the presence of oxygen. At high utilization rates of glutamine, oxidative utilization of glucose was decreased, while the production of lactate from glutamine was enhanced. This emergent phenotype was observed only when the available carbon exceeded the amount that could be fully oxidized by the available oxygen. Under the latter conditions, standard Flux Balance Analysis indicated that: this metabolic pattern is optimal to maximize biomass and ATP production; it requires the activity of a branched TCA cycle, in which glutamine-dependent reductive carboxylation cooperates to the production of lipids and proteins; it is sustained by a variety of redox-controlled metabolic reactions. In a K-ras transformed cell line we experimentally assessed glutamine-induced metabolic changes. We validated computational results through an extension of Flux Balance Analysis that allows prediction of metabolite variations. Taken together these findings offer new understanding of the logic of the metabolic reprogramming that underlies cancer cell growth.

Hallmarks describing common key events in initiation, maintenance and progression of cancer have been identified. One hallmark deals with rewiring of metabolic reactions required to sustain enhanced cell proliferation. The availability of molecular, mechanistic models of cancer hallmarks will mightily improve optimized personal treatment and new drug discovery. Metabolism is the only hallmark for which it is currently possible to derive large scale mathematical models, which have predictive ability. In this paper, we exploit a constraint-based model of the core metabolism required for biomass conversion of the most relevant nutrients—glucose and glutamine—to clarify the logic of control of cancer metabolism. We newly report that, when available oxygen is not sufficient to fully oxidize available glucose and glutamine carbons–a situation compatible with that observed under normal oxygen conditions in human and in cancer cells growing in vitro—utilization of glutamine by reductive carboxylation and conversion of glucose and glutamine to lactate confer advantage for biomass production. Redox homeostasis can be maintained through the use of different alternative pathways. In conclusion, this paper offers a logic interpretation to the link between metabolic rewiring and enhanced proliferation, which may offer new approaches to targeted drug discovery and utilization.

Zooming-in on cancer metabolic rewiring with tissue specific constraint-based models

The metabolic rearrangements occurring in cancer cells can be effectively investigated with a Systems Biology approach supported by metabolic network modeling. We here present tissue-specific constraint-based core models for three different types of tumors (liver, breast and lung) that serve this purpose. The core models were extracted and manually curated from the corresponding genome-scale metabolic models in the Human Metabolic Atlas database with a focus on the pathways that are known to play a key role in cancer growth and proliferation. Along similar lines, we also reconstructed a core model from the original general human metabolic network to be used as a reference model. A comparative Flux Balance Analysis between the reference and the cancer models highlighted both a clear distinction between the two conditions and a heterogeneity within the three different cancer types in terms of metabolic flux distribution. These results emphasize the need for modeling approaches able to keep up with this tumoral heterogeneity in order to identify more suitable drug targets and develop effective treatments. According to this perspective, we identified key points able to reverse the tumoral phenotype toward the reference one or vice-versa.

Computational strategies for a system-level understanding of metabolism

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism. To this aim, a plethora of methods have been developed, yet it is generally hard to identify which computational strategy is most suited for the investigation of a specific aspect of metabolism. This review provides an up-to-date description of the computational methods available for the analysis of metabolic pathways, discussing their main advantages and drawbacks. In particular, attention is devoted to the identification of the appropriate scale and level of accuracy in the reconstruction of metabolic networks, and to the inference of model structure and parameters, especially when dealing with a shortage of experimental measurements. The choice of the proper computational methods to derive in silico data is then addressed, including topological analyses, constraint-based modeling and simulation of the system dynamics. A description of some computational approaches to gain new biological knowledge or to formulate hypotheses is finally provided.

cuTauLeaping: a GPU-powered tau-leaping stochastic simulator for massive parallel analyses of biological systems

Tau-leaping is a stochastic simulation algorithm that efficiently reconstructs the temporal evolution of biological systems, modeled according to the stochastic formulation of chemical kinetics. The analysis of dynamical properties of these systems in physiological and perturbed conditions usually requires the execution of a large number of simulations, leading to high computational costs. Since each simulation can be executed independently from the others, a massive parallelization of tau-leaping can bring to relevant reductions of the overall running time. The emerging field of General Purpose Graphic Processing Units (GPGPU) provides power-efficient high-performance computing at a relatively low cost. In this work we introduce cuTauLeaping, a stochastic simulator of biological systems that makes use of GPGPU computing to execute multiple parallel tau-leaping simulations, by fully exploiting the Nvidia's Fermi GPU architecture. We show how a considerable computational speedup is achieved on GPU by partitioning the execution of tau-leaping into multiple separated phases, and we describe how to avoid some implementation pitfalls related to the scarcity of memory resources on the GPU streaming multiprocessors. Our results show that cuTauLeaping largely outperforms the CPU-based tau-leaping implementation when the number of parallel simulations increases, with a break-even directly depending on the size of the biological system and on the complexity of its emergent dynamics. In particular, cuTauLeaping is exploited to investigate the probability distribution of bistable states in the Schlögl model, and to carry out a bidimensional parameter sweep analysis to study the oscillatory regimes in the Ras/cAMP/PKA pathway in S. cerevisiae.

In vivo and in silico analysis of PCNA ubiquitylation in the activation of the Post Replication Repair pathway in S. cerevisiae

BACKGROUND

The genome of living organisms is constantly exposed to several damaging agents that induce different types of DNA lesions, leading to cellular malfunctioning and onset of many diseases. To maintain genome stability, cells developed various repair and tolerance systems to counteract the effects of DNA damage. Here we focus on Post Replication Repair (PRR), the pathway involved in the bypass of DNA lesions induced by sunlight exposure and UV radiation. PRR acts through two different mechanisms, activated by mono- and poly-ubiquitylation of the DNA sliding clamp, called Proliferating Cell Nuclear Antigen (PCNA).

RESULTS

We developed a novel protocol to measure the time-course ratios between mono-, di- and tri-ubiquitylated PCNA isoforms on a single western blot, which were used as the wet readout for PRR events in wild type and mutant S. cerevisiae cells exposed to acute UV radiation doses. Stochastic simulations of PCNA ubiquitylation dynamics, performed by exploiting a novel mechanistic model of PRR, well fitted the experimental data at low UV doses, but evidenced divergent behaviors at high UV doses, thus driving the design of further experiments to verify new hypothesis on the functioning of PRR. The model predicted the existence of a UV dose threshold for the proper functioning of the PRR model, and highlighted an overlapping effect of Nucleotide Excision Repair (the pathway effectively responsible to clean the genome from UV lesions) on the dynamics of PCNA ubiquitylation in different phases of the cell cycle. In addition, we showed that ubiquitin concentration can affect the rate of PCNA ubiquitylation in PRR, offering a possible explanation to the DNA damage sensitivity of yeast strains lacking deubiquitylating enzymes.

CONCLUSIONS

We exploited an in vivo and in silico combinational approach to analyze for the first time in a Systems Biology context the events of PCNA ubiquitylation occurring in PRR in budding yeast cells. Our findings highlighted an intricate functional crosstalk between PRR and other events controlling genome stability, and evidenced that PRR is more complicated and still far less characterized than previously thought.

The role of feedback control mechanisms on the establishment of oscillatory regimes in the Ras/cAMP/PKA pathway in S. cerevisiae

In the yeast Saccharomyces cerevisiae, the Ras/cAMP/PKA pathway is involved in the regulation of cell growth and proliferation in response to nutritional sensing and stress conditions. The pathway is tightly regulated by multiple feedback loops, exerted by the protein kinase A (PKA) on a few pivotal components of the pathway. In this article, we investigate the dynamics of the second messenger cAMP by performing stochastic simulations and parameter sweep analysis of a mechanistic model of the Ras/cAMP/PKA pathway, to determine the effects that the modulation of these feedback mechanisms has on the establishment of stable oscillatory regimes. In particular, we start by studying the role of phosphodiesterases, the enzymes that catalyze the degradation of cAMP, which represent the major negative feedback in this pathway. Then, we show the results on cAMP oscillations when perturbing the amount of protein Cdc25 coupled with the alteration of the intracellular ratio of the guanine nucleotides (GTP/GDP), which are known to regulate the switch of the GTPase Ras protein. This multi-level regulation of the amplitude and frequency of oscillations in the Ras/cAMP/PKA pathway might act as a fine tuning mechanism for the downstream targets of PKA, as also recently evidenced by some experimental investigations on the nucleocytoplasmic shuttling of the transcription factor Msn2 in yeast cells.

Simulation of the Ras/cAMP/PKA pathway in budding yeast highlights the establishment of stable oscillatory states

In the yeast Saccharomyces cerevisiae, the Ras/cAMP/PKA pathway plays a major role in the regulation of metabolism, stress resistance and cell cycle progression. We extend here a mechanistic model of the Ras/cAMP/PKA pathway that we previously defined by describing the molecular interactions and post-translational modifications of proteins, and perform a computational analysis to investigate the dynamical behaviors of the components of this pathway, regulated by different control mechanisms. We carry out stochastic simulations to consider, in particular, the effect of the negative feedback loops on the activity of both Ira2 (a Ras-GAP) and Cdc25 (a Ras-GEF) proteins. Our results show that stable oscillatory regimes for the dynamics of cAMP can be obtained only through the activation of these feedback mechanisms, and when the amount of Cdc25 is within a specific range. In addition, we highlight that the levels of guanine nucleotides pools are able to regulate the pathway, by influencing the transition between stable steady states and oscillatory regimes.

Modeling and stochastic simulation of the Ras/cAMP/PKA pathway in the yeast Saccharomyces cerevisiae evidences a key regulatory function for intracellular guanine nucleotides pools

In the yeast Saccharomyces cerevisiae, the Ras/cAMP/PKA pathway is involved in the regulation of metabolism and cell cycle progression. The pathway is tightly regulated by several control mechanisms, as the feedback cycle ruled by the activity of phosphodiesterase. Here, we present a discrete mathematical model for the Ras/cAMP/PKA pathway that considers its principal cytoplasmic components and their mutual interactions. The tau-leaping algorithm is then used to perform stochastic simulations of the model. We investigate this system under various conditions, and we test how different values of several stochastic reaction constants affect the pathway behaviour. Finally, we show that the level of guanine nucleotides, GTP and GDP, could be relevant metabolic signals for the regulation of the whole pathway.

Modelling metapopulations with stochastic membrane systems

Metapopulations, or multi-patch systems, are models describing the interactions and the behavior of populations living in fragmented habitats. Dispersal, persistence and extinction are some of the characteristics of interest in ecological studies of metapopulations. In this paper, we propose a novel method to analyze metapopulations, which is based on a discrete and stochastic modelling framework in the area of Membrane Computing. New structural features of membrane systems, necessary to appropriately describe a multi-patch system, are introduced, such as the reduction of the maximal parallel consumption of objects, the spatial arrangement of membranes and the stochastic creation of objects. The role of the additional features, their meaning for a metapopulation model and the emergence of relevant behaviors are then investigated by means of stochastic simulations. Conclusive remarks and ideas for future research are finally presented.

Heterodyne near-field scattering: a technique for complex fluids

We implemented the heterodyne near-field scattering (HNFS) technique [Appl. Phys. Lett. 81, 4109 (2002)]], showing that it is a fairly valid alternative to traditional elastic low-angle light scattering and quite suitable for studying complex fluids such as colloidal systems. With respect to the original work, we adopted a different data reduction scheme, which allowed us to improve significantly the performance of the technique, at levels of sensitivity and accuracy much higher than those achievable with classical low-angle light scattering instrumentation. This method also relaxes the requirements on the optical/mechanical stability of the experimental setup and allows for a real time analysis. The HNFS technique has been tested by using calibrated colloidal particles and its capability of performing accurate particle sizing was ascertained on both monodisperse and bimodal particle distributions. Nonstationary samples, such as aggregating colloidal solutions, were profitably studied, and their kinetics quantitatively characterized.