PeTTSy: a computational tool for perturbation analysis of complex systems biology models.

BACKGROUND: Over the last decade sensitivity analysis techniques have been shown to be very useful to analyse complex and high dimensional Systems Biology models. However, many of the currently available toolboxes have either used parameter sampling, been focused on a restricted set of model observables of interest, studied optimisation of a objective function, or have not dealt with multiple simultaneous model parameter changes where the changes can be permanent or temporary. RESULTS: Here we introduce our new, freely downloadable toolbox, PeTTSy (Perturbation Theory Toolbox for Systems). PeTTSy is a package for MATLAB which implements a wide array of techniques for the perturbation theory and sensitivity analysis of large and complex ordinary differential equation (ODE) based models. PeTTSy is a comprehensive modelling framework that introduces a number of new approaches and that fully addresses analysis of oscillatory systems. It examines sensitivity analysis of the models to perturbations of parameters, where the perturbation timing, strength, length and overall shape can be controlled by the user. This can be done in a system-global setting, namely, the user can determine how many parameters to perturb, by how much and for how long. PeTTSy also offers the user the ability to explore the effect of the parameter perturbations on many different types of outputs: period, phase (timing of peak) and model solutions. PeTTSy can be employed on a wide range of mathematical models including free-running and forced oscillators and signalling systems. To enable experimental optimisation using the Fisher Information Matrix it efficiently allows one to combine multiple variants of a model (i.e. a model with multiple experimental conditions) in order to determine the value of new experiments. It is especially useful in the analysis of large and complex models involving many variables and parameters. CONCLUSIONS: PeTTSy is a comprehensive tool for analysing large and complex models of regulatory and signalling systems. It allows for simulation and analysis of models under a variety of environmental conditions and for experimental optimisation of complex combined experiments. With its unique set of tools it makes a valuable addition to the current library of sensitivity analysis toolboxes. We believe that this software will be of great use to the wider biological, systems biology and modelling communities.

Impact of Serotonin Transporter Absence on Brain Insulin Receptor Expression, Plasma Metabolome Changes, and ADHD-like Behavior in Mice fed a Western Diet.

The impaired function of the serotonin transporter (SERT) in humans has been linked to a higher risk of obesity and type 2 diabetes, especially as people age. Consuming a "Western diet" (WD), which is high in saturated fats, cholesterol, and sugars, can induce metabolic syndrome. Previous research indicated that mice carrying a targeted inactivation of the Sert gene (knockout, KO) and fed a WD display significant metabolic disturbances and behaviors reminiscent of ADHD. These abnormalities might be mediated via a dysfunction in insulin receptor (IR) signaling, which is also associated with adult ADHD. However, the impact of Sert deficiency on IR signaling and systemic metabolic changes has not been thoroughly explored. In this study, we conducted a detailed analysis of locomotor behavior in wild-type (WT) and KO mice fed a WD or control diet. We investigated changes in the blood metabolome and examined, via PCR, the expression of insulin receptor A and B isoforms and key regulators of their function in the brain. Twelve-month-old KO mice and their WT littermates were fed a WD for three weeks. Nuclear magnetic resonance spectroscopy analysis of plasma samples showed that KO mice on a WD had higher levels of lipids and lipoproteins and lower levels of glucose, lactate, alanine, valine, and isoleucine compared to other groups. SERT-KO mice on the control diet exhibited increased brain levels of both IR A and B isoforms, accompanied by a modest increase in the negative regulator ENPP. The KO mice also displayed anxiety-like behavior and reduced exploratory activity in an open field test. However, when the KO animals were fed a WD, the aberrant expression levels of IR isoforms in the KO mice and locomotor behavior were ameliorated indicating a complex interaction between genetic and dietary factors that might contribute to ADHD-like symptoms. Overall, our findings suggest that the lack of Sert leads to a unique metabolic phenotype in aged mice, characterized by dysregulated IR-related pathways. These changes are exacerbated by WD in the blood metabolome and are associated with behavioral abnormalities.

Genotypes of the UCP1 gene polymorphisms and cardiometabolic diseases: A multifactorial study of association with disease probability.

Cardiometabolic diseases (CMDs) are complex disorders with a heterogenous phenotype, which are caused by multiple factors including genetic factors. Single nucleotide polymorphisms (SNPs) rs45539933 (p.Ala64Thr), rs10011540 (c.-112A>C), rs3811791 (c.-1766A>G), and rs1800592 (c.-3826A>G) in the UCP1 gene have been analyzed for association with CMDs in many studies providing controversial results. However, previous studies only considered individual UCP1 SNPs and did not evaluate them in an integrated manner, which is a more powerful approach to uncover genetic component of complex diseases. This study aimed to investigate associations between UCP1 genotype combinations and CMDs or CMD risk factors in the context of non-genetic factors. We performed multiple logistic regression analysis and proposed new methodology of testing different combinations of SNP genotypes. We found that probability of CMDs increased in presence of the three-SNP combination of genotypes with minor alleles of c.-3826A>G and p.Ala64Thr and wild allele of c.-112A>C, with increasing age, body mass index (BMI), body fat percentage (BF%) and may differ between sexes and between countries. The combination of genotypes with c.-3826A>G minor allele and wild homozygotes of c.-112A>C and p.Ala64Thr was associated with increased probability of diabetes. While combination of genotypes with minor alleles of all three SNPs reduced the CMD probability. The present results suggest that age, BMI, sex, and UCP1 three-SNP combinations of genotypes significantly contribute to CMD probability. Varying of c.-112A>C alleles in the genotype combination with minor alleles of c.-3826A>G and p.Ala64Thr markedly changes CMD probability.

Optimization of the MACE endpoint composition to increase power in studies of lipid-lowering therapies-a model-based meta-analysis.

Aims: To develop a model-informed methodology for the optimization of the Major Adverse Cardiac Events (MACE) composite endpoint, based on a model-based meta-analysis across anti-hypercholesterolemia trials of statin and anti-PCSK9 drugs. Methods and results: Mixed-effects meta-regression modeling of stand-alone MACE outcomes was performed, with therapy type, population demographics, baseline and change over time in lipid biomarkers as predictors. Randomized clinical trials up to June 28, 2022, of either statins or anti-PCSK9 therapies were identified through a systematic review process in PubMed and ClinicalTrials.gov databases. In total, 54 studies (270,471 patients) were collected, reporting 15 different single cardiovascular events. Treatment-mediated decrease in low density lipoprotein cholesterol, baseline levels of remnant and high-density lipoprotein cholesterol as well as non-lipid population characteristics and type of therapy were identified as significant covariates for 10 of the 15 outcomes. The required sample size per composite 3- and 4-point MACE endpoint was calculated based on the estimated treatment effects in a population and frequencies of the incorporated events in the control group, trial duration, and uncertainty in model parameters. Conclusion: A quantitative tool was developed and used to benchmark different compositions of 3- and 4-point MACE for statins and anti-PCSK9 therapies, based on the minimum population size required to achieve statistical significance in relative risk reduction, following meta-regression modeling of the single MACE components. The approach we developed may be applied towards the optimization of the design of future trials in dyslipidemia disorders as well as in other therapeutic areas.

Dose dependence of treatment-related adverse events for immune checkpoint inhibitor therapies: a model-based meta-analysis.

Programmed cell death-1 (PD-1) and/or cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) immune checkpoint inhibitor (ICI) treatments are associated with adverse events (AEs), which may be dependent on ICI dose. Applying a model-based meta-analysis to evaluate safety data from published clinical trials from 2005 to 2018, we analyzed the dose/exposure dependence of ICI treatment-related AE (trAE) and immune-mediated AE (imAE) rates. Unlike with PD-1 inhibitor monotherapy, CTLA-4 inhibitor monotherapy exhibited a dose/exposure dependence on most AE types evaluated. Furthermore, combination therapy with PD-1 inhibitor significantly strengthened the dependence of trAE and imAE rates on CTLA-4 inhibitor dose/exposure.

A Generic Mechanism for Enhanced Cytokine Signaling via Cytokine-Neutralizing Antibodies.

Enhancement or inhibition of cytokine signaling and corresponding immune cells responses are critical factors in various disease treatments. Cytokine signaling may be inhibited by cytokine-neutralizing antibodies (CNAs), which prevents further activation of cytokine receptors. However, CNAs may result in enhanced-instead of inhibitory-cytokine signaling (an "agonistic effect") in various in vitro and in vivo experiments. This may lead to lack of efficacy or adverse events for cytokine-inhibiting based medicines. Alternatively, cytokine-antibody complexes may produce stronger signaling vs. cytokine alone, thereby increasing the efficacy of stimulating cytokine-based drugs, at equal or lower cytokine doses. In this paper, the effect of cytokine signaling enhancement by a CNA was studied in a generic mathematical model of interleukin-4 (IL-4) driven T-cell proliferation. The occurrence of the agonistic effect depends upon the antibody-to-cytokine binding affinity and initial concentrations of antibody and cytokine. Model predictions were in agreement with experimental studies. When the cytokine receptor consists of multiple subunits with substantially differing affinities (e.g., IL-4 case), the choice of the receptor chain to be blocked by the antibody is critical, for the agonistic effect to appear. We propose a generic mechanism for the effect: initially, binding of the CNA to the cytokine reduces free cytokine concentration; yet, cytokine molecules bound within the cytokine-CNA complex-and released later and over time-are "rescued" from earlier clearance via cellular internalization. Hence, although free cytokine-dependent signalling may be less potent initially, it will also be more sustained over time; and given non-linear dynamics, it will lead ultimately to larger cellular effector responses, vs. the same amount of free cytokine in the absence of CNA. We suggest that the proposed mechanism is a generic property of {cytokine, CNA, receptor} triads, both in vitro and in vivo, and can occur in a predictable fashion for a variety of cytokines of the immune system.