Combination therapies proved to be a valuable strategy in the fight against cancer, thanks to their increased efficacy in inducing tumor cell death and in reducing tumor growth, metastatic potential, ...and the risk of developing drug resistance. The identification of effective combinations of drug targets generally relies on costly and time consuming processes based on
experiments. Here, we present a novel computational approach that, by integrating dynamic fuzzy modeling with multi-objective optimization, allows to efficiently identify novel combination cancer therapies, with a relevant saving in working time and costs. We tested this approach on a model of oncogenic K-ras cancer cells characterized by a marked Warburg effect. The computational approach was validated by its capability in finding out therapies already known in the literature for this type of cancer cell. More importantly, our results show that this method can suggest potential therapies consisting in a small number of molecular targets. In the model of oncogenic K-ras cancer cells, for instance, we identified combination of up to three targets, which affect different cellular pathways that are crucial for cancer proliferation and survival.
Several software tools for the simulation and analysis of biochemical reaction networks have been developed in the last decades; however, assessing and comparing their computational performance in ...executing the typical tasks of computational systems biology can be limited by the lack of a standardized benchmarking approach. To overcome these limitations, we propose here a novel tool, named SMGen, designed to automatically generate synthetic models of reaction networks that, by construction, are characterized by relevant features (e.g., system connectivity and reaction discreteness) and non-trivial emergent dynamics of real biochemical networks. The generation of synthetic models in SMGen is based on the definition of an undirected graph consisting of a single connected component that, generally, results in a computationally demanding task; to speed up the overall process, SMGen exploits a main–worker paradigm. SMGen is also provided with a user-friendly graphical user interface, which allows the user to easily set up all the parameters required to generate a set of synthetic models with any number of reactions and species. We analysed the computational performance of SMGen by generating batches of symmetric and asymmetric reaction-based models (RBMs) of increasing size, showing how a different number of reactions and/or species affects the generation time. Our results show that when the number of reactions is higher than the number of species, SMGen has to identify and correct a large number of errors during the creation process of the RBMs, a circumstance that increases the running time. Still, SMGen can generate synthetic models with hundreds of species and reactions in less than 7 s.
Calcium homeostasis and signaling processes in
, as well as in any eukaryotic organism, depend on various transporters and channels located on both the plasma and intracellular membranes. The ...activity of these proteins is regulated by a number of feedback mechanisms that act through the calmodulin-calcineurin pathway. When exposed to hypotonic shock (HTS), yeast cells respond with an increased cytosolic calcium transient, which seems to be conditioned by the opening of stretch-activated channels. To better understand the role of each channel and transporter involved in the generation and recovery of the calcium transient-and of their feedback regulations-we defined and analyzed a mathematical model of the calcium signaling response to HTS in yeast cells. The model was validated by comparing the simulation outcomes with calcium concentration variations before and during the HTS response, which were observed experimentally in both wild-type and mutant strains. Our results show that calcium normally enters the cell through the High Affinity Calcium influx System and mechanosensitive channels. The increase of the plasma membrane tension, caused by HTS, boosts the opening probability of mechanosensitive channels. This event causes a sudden calcium pulse that is rapidly dissipated by the activity of the vacuolar transporter Pmc1. According to model simulations, the role of another vacuolar transporter, Vcx1, is instead marginal, unless calcineurin is inhibited or removed. Our results also suggest that the mechanosensitive channels are subject to a calcium-dependent feedback inhibition, possibly involving calmodulin. Noteworthy, the model predictions are in accordance with literature results concerning some aspects of calcium homeostasis and signaling that were not specifically addressed within the model itself, suggesting that it actually depicts all the main cellular components and interactions that constitute the HTS calcium pathway, and thus can correctly reproduce the shaping of the calcium signature by calmodulin- and calcineurin-dependent complex regulations. The model predictions also allowed to provide an interpretation of different regulatory schemes involved in calcium handling in both wild-type and mutants yeast strains. The model could be easily extended to represent different calcium signals in other eukaryotic cells.
Oxford Nanopore Technologies (ONT) is a third generation sequencing approach that allows the analysis of individual, full-length nucleic acids. ONT records the alterations of an ionic current flowing ...across a nano-scaled pore while a DNA or RNA strand is threading through the pore. Basecalling methods are then leveraged to translate the recorded signal back to the nucleic acid sequence. However, basecall generally introduces errors that hinder the process of barcode demultiplexing, a pivotal task in single-cell RNA sequencing that allows for separating the sequenced transcripts on the basis of their cell of origin.
To solve this issue, we present a novel framework, called UNPLEX, designed to tackle the barcode demultiplexing problem by operating directly on the recorded signals. UNPLEX combines two unsupervised machine learning methods: autoencoders and self-organizing maps (SOM). The autoencoders extract compact, latent representations of the recorded signals that are then clustered by the SOM.
Our results, obtained on two datasets composed of
generated ONT-like signals, show that UNPLEX represents a promising starting point for the development of effective tools to cluster the signals corresponding to the same cell.
Ras oncoproteins play a crucial role in the onset, maintenance, and progression of the most common and deadly human cancers. Despite extensive research efforts, only a few mutant-specific Ras ...inhibitors have been reported. We show that cmp4-previously identified as a water-soluble Ras inhibitor- targets multiple steps in the activation and downstream signaling of different Ras mutants and isoforms. Binding of this pan-Ras inhibitor to an extended Switch II pocket on HRas and KRas proteins induces a conformational change that down-regulates intrinsic and GEF-mediated nucleotide dissociation and exchange and effector binding. A mathematical model of the Ras activation cycle predicts that the inhibitor severely reduces the proliferation of different Ras-driven cancer cells, effectively cooperating with Cetuximab to reduce proliferation even of Cetuximab-resistant cancer cell lines. Experimental data confirm the model prediction, indicating that the pan-Ras inhibitor is an appropriate candidate for medicinal chemistry efforts tailored at improving its currently unsatisfactory affinity.
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.
Advances in microscopy imaging technologies have enabled the visualization of live-cell dynamic processes using time-lapse microscopy imaging. However, modern methods exhibit several limitations ...related to the training phases and to time constraints, hindering their application in the laboratory practice. In this work, we present a novel method, named Automated Cell Detection and Counting (ACDC), designed for activity detection of fluorescent labeled cell nuclei in time-lapse microscopy. ACDC overcomes the limitations of the literature methods, by first applying bilateral filtering on the original image to smooth the input cell images while preserving edge sharpness, and then by exploiting the watershed transform and morphological filtering. Moreover, ACDC represents a feasible solution for the laboratory practice, as it can leverage multi-core architectures in computer clusters to efficiently handle large-scale imaging datasets. Indeed, our Parent-Workers implementation of ACDC allows to obtain up to a 3.7× speed-up compared to the sequential counterpart. ACDC was tested on two distinct cell imaging datasets to assess its accuracy and effectiveness on images with different characteristics. We achieved an accurate cell-count and nuclei segmentation without relying on large-scale annotated datasets, a result confirmed by the average Dice Similarity Coefficients of 76.84 and 88.64 and the Pearson coefficients of 0.99 and 0.96, calculated against the manual cell counting, on the two tested datasets.
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.
Among the existing global optimization algorithms, Particle Swarm Optimization (PSO) is one of the most effective methods for non-linear and complex high-dimensional problems. Since PSO performance ...strongly depends on the choice of its settings (i.e., inertia, cognitive and social factors, minimum and maximum velocity), Fuzzy Logic (FL) was previously exploited to select these values. So far, FL-based implementations of PSO aimed at the calculation of a unique settings for the whole swarm. In this work we propose a novel self-tuning algorithm—called Fuzzy Self-Tuning PSO (FST-PSO)—which exploits FL to calculate the inertia, cognitive and social factor, minimum and maximum velocity independently for each particle, thus realizing a complete settings-free version of PSO. The novelty and strength of FST-PSO lie in the fact that it does not require any expertise in PSO functioning, since the behavior of every particle is automatically and dynamically adjusted during the optimization. We compare the performance of FST-PSO with standard PSO, Proactive Particles in Swarm Optimization, Artificial Bee Colony, Covariance Matrix Adaptation Evolution Strategy, Differential Evolution and Genetic Algorithms. We empirically show that FST-PSO can basically outperform all tested algorithms with respect to the convergence speed and is competitive concerning the best solutions found, noticeably with a reduced computational effort.