Gene expression is the result of the balance between transcription and degradation. Recent experimental findings have shown fine and specific regulation of RNA degradation and the presence of various ...molecular machinery purposely devoted to this task, such as RNA binding proteins, non-coding RNAs, etc. A biological process can be studied by measuring time-courses of RNA abundance in response of internal and/or external stimuli, using recent technologies, such as the microarrays or the Next Generation Sequencing devices. Unfortunately, the picture provided by looking only at the transcriptome abundance may not gain insight into its dynamic regulation. By contrast, independent simultaneous measurement of RNA expression and half-lives could provide such valuable additional insight. A computational approach to the estimation of RNAs half-lives from RNA expression time profiles data, can be a low-cost alternative to its experimental measurement which may be also affected by various artifacts.
Here we present a computational methodology, called StaRTrEK (STAbility Rates ThRough Expression Kinetics), able to estimate half-life values basing only on genome-wide gene expression time series without transcriptional inhibition. The StaRTrEK algorithm makes use of a simple first order kinetic model and of a Formula: see text-norm regularized least square optimization approach to find its parameter values. Estimates provided by StaRTrEK are validated using simulated data and three independent experimental datasets of two short (6 samples) and one long (48 samples) time-courses.
We believe that our algorithm can be used as a fast valuable computational complement to time-course experimental gene expression studies by adding a relevant kinetic property, i.e. the RNA half-life, with a strong biological interpretation, thus providing a dynamic picture of what is going in a cell during the biological process under study.
An epidemic multi-group model formed by interconnected SEIR-like structures is formulated and used for data fitting to gain insight into the COVID-19 dynamics and into the role of non-pharmaceutical ...control actions implemented to limit the infection spread since its outbreak in Italy. The single submodels provide a rather accurate description of the COVID-19 evolution in each subpopulation by an extended SEIR model including the class of asymptomatic infectives, which is recognized as a determinant for disease diffusion. The multi-group structure is specifically designed to investigate the effects of the inter-regional mobility restored at the end of the first strong lockdown in Italy (June 3, 2020). In its time-invariant version, the model is shown to enjoy some analytical stability properties which provide significant insights on the efficacy of the implemented control measurements. In order to highlight the impact of human mobility on the disease evolution in Italy between the first and second wave onset, the model is applied to fit real epidemiological data of three geographical macro-areas in the period March–October 2020, including the mass departure for summer holidays. The simulation results are in good agreement with the data, so that the model can represent a useful tool for predicting the effects of the combination of containment measures in triggering future pandemic scenarios. Particularly, the simulation shows that, although the unrestricted mobility alone appears to be insufficient to trigger the second wave, the human transfers were crucial to make uniform the spatial distribution of the infection throughout the country and, combined with the restart of the production, trade, and education activities, determined a time advance of the contagion increase since September 2020.
The recent outbreak of COVID-19 underlined the need for a fast and trustworthy methodology to identify the features of a pandemic, whose early identification is of help for designing ...non-pharmaceutical interventions (including lockdown and social distancing) to limit the progression of the disease. A common approach in this context is the parameter identification from deterministic epidemic models, which, unfortunately, cannot take into account the inherent randomness of the epidemic phenomenon, especially in the initial stage; on the other hand, the use of raw data within the framework of a stochastic model is not straightforward. This note investigates the stochastic approach applied to a basic SIR (Susceptible, Infected, Recovered) epidemic model to enhance information from raw data generated in silico. The stochastic model consists of a Continuous-Time Markov Model, describing the epidemic outbreak in terms of stochastic discrete infection and recovery events in a given region, and where independent random paths are associated to different provinces of the same region, which are assumed to share the same set of model parameters. The estimation procedure is based on the building of a loss function that symmetrically weighs first-order and second-order moments, differently from the standard approach that considers a highly asymmetrical choice, exploiting only first-order moments. Instead, we opt for an innovative symmetrical identification approach which exploits both moments. The new approach is specifically proposed to enhance the statistical information content of the raw epidemiological data.
Insulin resistance is the common denominator of several diseases including type 2 diabetes and cancer, and investigating the mechanisms responsible for insulin signaling impairment is of primary ...importance. A mathematical model of the insulin signaling network (ISN) is proposed and used to investigate the dose-response curves of components of this network. Experimental data of C2C12 myoblasts with phosphatase and tensin homologue (PTEN) suppressed and data of L6 myotubes with induced insulin resistance have been analyzed by the model. We focused particularly on single and double Akt phosphorylation and pointed out insulin signaling changes related to insulin resistance. Moreover, a new characterization of the upstream signaling of the mammalian target of rapamycin complex 2 (mTORC2) is presented. As it is widely recognized that ISN proteins have a crucial role also in cell proliferation and death, the ISN model was linked to a cell population model and applied to data of a cell line of acute myeloid leukemia treated with a mammalian target of rapamycin inhibitor with antitumor activity. The analysis revealed simple relationships among the concentrations of ISN proteins and the parameters of the cell population model that characterize cell cycle progression and cell death.
Current train delay (TD) prediction systems do not take advantage of state-of-the-art tools and techniques for handling and extracting useful and actionable information from the large amount of ...endogenous (i.e., generated by the railway system itself) and exogenous (i.e., related to railway operation but generated by external phenomena) data available. Additionally, they are not designed in order to deal with the intrinsic time varying nature of the problem (e.g., regular changes in the nominal timetable, etc.). The purpose of this paper is to build a dynamic data-driven TD prediction system that exploits the most recent tools and techniques in the field of time varying big data analysis. In particular, we map the TD prediction problem into a time varying multivariate regression problem that allows exploiting both historical data about the train movements and exogenous data about the weather provided by the national weather services. The performance of these methods have been tuned through the state-of-the-art thresholdout technique, a very powerful procedure which relies on the differential privacy theory. Finally, the performance of two efficient implementations of shallow and deep extreme learning machines that fully exploit the recent in-memory large-scale data processing technologies have been compared with the current state-of-the-art TD prediction systems. Results on real-world data coming from the Italian railway network show that the proposal of this paper is able to remarkably improve the state-of-the-art systems.
In this work, we investigate some theoretical aspects related to the estimation approach proposed by Liebermeister and Klipp, 2006, in which general rate laws, derived from standardized enzymatic ...mechanisms, are exploited to kinetically describe the fluxes of a metabolic reaction network, and multiple metabolic steady-state measurements are exploited to estimate the unknown kinetic parameters. Further mathematical details are deeply investigated, and necessary conditions on the amount of information required to solve the identification problem are given. Moreover, theoretical results for the parameter identifiability are provided, and symmetrical and modular properties of the proposed approach are highlighted when the global identification problem is decoupled into smaller and simpler identification problems related to the single reactions of the network. Among the advantages of the proposed innovative approach are (i) non-restrictive conditions to guarantee the solvability of the parameter estimation problem, (ii) the unburden of the usual computational complexity for such identification problems, and (iii) the ease of obtaining the required number of measurements, which are actually steady-state data, experimentally easier to obtain with respect to the time-dependent ones. A simple example concludes the paper, highlighting the mentioned advantages of the method and the implementation of the related theoretical result.
The present work of review collects and evidences the main results of our previous papers on the optimization of fractionated radiotherapy protocols. The problem under investigation is presented here ...in a unitary framework as a nonlinear programming application that aims to determine the optimal schemes of dose fractionation commonly used in external beam radiotherapy. The radiation responses of tumor and normal tissues are described by means of the linear quadratic model. We formulate a nonlinear, non-convex optimization problem including two quadratic constraints to limit the collateral normal tissue damages and linear box constraints on the fractional dose sizes. The general problem is decomposed into two subproblems: (1) analytical determination of the optimal fraction dose sizes as a function of the model parameters for arbitrarily fixed treatment lengths; and (2) numerical determination of the optimal fraction number, and of the optimal treatment time, in different parameter settings. After establishing the boundedness of the optimal number of fractions, we investigate by numerical simulation the optimal solution behavior for experimentally meaningful parameter ranges, recognizing the crucial role of some parameters, such as the radiosensitivity ratio, in determining the optimality of hypo- or equi-fractionated treatments. Our results agree with findings of the theoretical and clinical literature.
(1) Background: Plantar fasciitis (PF) is the most common cause of heel pain in adults. Extracorporeal shockwave therapy (ESWT) is a minimally invasive treatments commonly used for treating PF. Our ...aim is to provide a complete overview of which treatments have been compared to ESWT, with a focus on the modalities that have been used. (2) Methods: A thorough search of the literature was performed on Medline via Pubmed, Cochrane Database of Systematic Reviews (CDSR) of the Cochrane Library and Physiotherapy Evidence Databases (PEDro) up to 18 November 2021. In the study were included only systematic reviews and meta-analysis in English language, published from 2010 to date. (3) Results: A total of 14 systematic reviews and meta-analysis were included in the umbrella review. A total of eight studies compared the efficacy of ESWT treatment with placebo, three studies compared ESWT with another therapy (two studies compared ESWT and corticosteroids, one study ESWT and ultrasound therapy), and three studies had more than one comparison. (4) Conclusions: When compared to placebo, ESWT demonstrated to be effective. More randomized trials with specific comparisons between different types and intensity of SW are needed to obtain more precise information on SW effectiveness.
Mathematical modeling is a powerful approach supporting the investigation of open problems in natural sciences, in particular physics, biology and medicine. Applied mathematics allows to translate ...the available information about real-world phenomena into mathematical objects and concepts. Mathematical models are useful descriptive tools that allow to gather the salient aspects of complex biological systems along with their fundamental governing laws, by elucidating the system behavior in time and space, also evidencing symmetry, or symmetry breaking, in geometry and morphology. Additionally, mathematical models are useful predictive tools able to reliably forecast the future system evolution or its response to specific inputs. More importantly, concerning biomedical systems, such models can even become prescriptive tools, allowing effective, sometimes optimal, intervention strategies for the treatment and control of pathological states to be planned. The application of mathematical physics, nonlinear analysis, systems and control theory to the study of biological and medical systems results in the formulation of new challenging problems for the scientific community. This Special Issue includes innovative contributions of experienced researchers in the field of mathematical modelling applied to biology and medicine.
In this work we address the problem of tumour growth control by properly exploiting a low-dimensional model that grounds on the Chemical Reaction Network (CRN) formalism. Originally conceived to work ...both in deterministic and stochastic frameworks, it is shown that, except for the case of very low number of tumour cells, the deterministic approach is appropriate to characterize the system behaviour, especially for control planning purposes. Two alternative control approaches are here investigated. One trivially assumes a constant infusion of external drug administration, the other is designed according to a state-feedback control scheme, with complete or partial knowledge of the state. Pros and cons of both control laws are investigated, showing that the tumour size at the beginning of the therapy plays a role of paramount importance for fixed infusion therapies, whilst only state-feedback laws can eradicate arbitrarily large tumours.
•We perform a qualitative analysis of a compact ODE tumour growth model.•Two alternative control approaches are investigated, with full/partial information.•Too large tumours have no chance to be tackled by constant infusion therapies.•The proposed state-feedback control is able to eradicate arbitrarily large tumours.