Search

Abstract Background Mechanistic models are becoming more and more popular in Systems Biology; identification and control of models underlying biochemical pathways of interest in oncology is a primary goal in this field. Unfortunately the scarce availability of data still limits our understanding of the intrinsic characteristics of complex pathologies like cancer: acquiring information for a system understanding of complex reaction networks is time consuming and expensive. Stimulus response experiments (SRE) have been used to gain a deeper insight into the details of biochemical mechanisms underlying cell life and functioning. Optimisation of the input time-profile, however, still remains a major area of research due to the complexity of the problem and its relevance for the task of information retrieval in systems biology-related experiments. Results We have addressed the problem of quantifying the information associated to an experiment using the Fisher Information Matrix and we have proposed an optimal experimental design strategy based on evolutionary algorithm to cope with the problem of information gathering in Systems Biology. On the basis of the theoretical results obtained in the field of control systems theory, we have studied the dynamical properties of the signals to be used in cell stimulation. The results of this study have been used to develop a microfluidic device for the automation of the process of cell stimulation for system identification. Conclusion We have applied the proposed approach to the Epidermal Growth Factor Receptor pathway and we observed that it minimises the amount of parametric uncertainty associated to the identified model. A statistical framework based on Monte-Carlo estimations of the uncertainty ellipsoid confirmed the superiority of optimally designed experiments over canonical inputs. The proposed approach can be easily extended to multiobjective formulations that can also take advantage of identifiability analysis. Moreover, the availability of fully automated microfluidic platforms explicitly developed for the task of biochemical model identification will hopefully reduce the effects of the 'data rich-data poor' paradox in Systems Biology.

Pompe disease, a rare genetic neuromuscular disorder, is caused by a deficiency of acid alpha-glucosidase (GAA), leading to an accumulation of glycogen in lysosomes, and resulting in the progressive development of muscle weakness. The current standard treatment, enzyme replacement therapy (ERT), is not curative and has limitations such as poor penetration into skeletal muscle and both the central and peripheral nervous systems, a risk of immune responses against the recombinant enzyme, and the requirement for high doses and frequent infusions. To overcome these limitations, lentiviral vector-mediated hematopoietic stem and progenitor cell (HSPC) gene therapy has been proposed as a next-generation approach for treating Pompe disease. This study demonstrates the potential of lentiviral HSPC gene therapy to reverse the pathological effects of Pompe disease in a preclinical mouse model. It includes a comprehensive safety assessment via integration site analysis, along with single-cell RNA sequencing analysis of central nervous tissue samples to gain insights into the underlying mechanisms of phenotype correction., Competing Interests: Declaration of interests All authors were former employees of AVROBIO, Inc., Cambridge, MA, USA during the conception and writing of the manuscript, except V.M., M.R., and A.S. N.P.v.T. and C.M. are inventors on patents in the field of HSC gene therapy. AVROBIO, Inc., has a preclinical gene therapy program for Pompe disease (AVR-RD-03) based on a genetically modified HSPC platform using lentiviral vectors. Collection of data and analysis was performed as part of the program. This research received no external funding and was sponsored by AVROBIO, Inc., (Copyright © 2024 The American Society of Gene and Cell Therapy. Published by Elsevier Inc. All rights reserved.)

Pompe disease is a rare genetic neuromuscular disorder caused by acid α-glucosidase (GAA) deficiency resulting in lysosomal glycogen accumulation and progressive myopathy. Enzyme replacement therapy, the current standard of care, penetrates poorly into the skeletal muscles and the peripheral and central nervous system (CNS), risks recombinant enzyme immunogenicity, and requires high doses and frequent infusions. Lentiviral vector-mediated hematopoietic stem and progenitor cell (HSPC) gene therapy was investigated in a Pompe mouse model using a clinically relevant promoter driving nine engineered GAA coding sequences incorporating distinct peptide tags and codon optimizations. Vectors solely including glycosylation-independent lysosomal targeting tags enhanced secretion and improved reduction of glycogen, myofiber, and CNS vacuolation in key tissues, although GAA enzyme activity and protein was consistently lower compared with native GAA. Genetically modified microglial cells in brains were detected at low levels but provided robust phenotypic correction. Furthermore, an amino acid substitution introduced in the tag reduced insulin receptor-mediated signaling with no evidence of an effect on blood glucose levels in Pompe mice. This study demonstrated the therapeutic potential of lentiviral HSPC gene therapy exploiting optimized GAA tagged coding sequences to reverse Pompe disease pathology in a preclinical mouse model, providing promising vector candidates for further investigation., Competing Interests: All authors are current or former employees of AVROBIO, Inc., Cambridge, MA, USA during the conception and writing of the manuscript, except A.S. AVROBIO, Inc. has a preclinical gene therapy program for Pompe disease (AVR-RD-03) based on a genetically modified HSPC platform using lentiviral vectors. Collection of data and analysis was performed as part of the program. This research received no external funding and was sponsored by AVROBIO, Inc., (© 2022 The Authors.)

Background: Even though clinical trials are indispensable for medical research, they are frequently impaired by delayed or incomplete patient recruitment, resulting in cost overruns or aborted studies. Study protocols based on real-world data with precisely expressed eligibility criteria and realistic cohort estimations are crucial for successful study execution. The increasing availability of routine clinical data in electronic health records (EHRs) provides the opportunity to also support patient recruitment during the prescreening phase. While solutions for electronic recruitment support have been published, to our knowledge, no method for the prioritization of eligibility criteria in this context has been explored., Methods: In the context of the Electronic Health Records for Clinical Research (EHR4CR) project, we examined the eligibility criteria of the KATHERINE trial. Criteria were extracted from the study protocol, deduplicated, and decomposed. A paper chart review and data warehouse query were executed to retrieve clinical data for the resulting set of simplified criteria separately from both sources. Criteria were scored according to disease specificity, data availability, and discriminatory power based on their content and the clinical dataset., Results: The study protocol contained 35 eligibility criteria, which after simplification yielded 70 atomic criteria. For a cohort of 106 patients with breast cancer and neoadjuvant treatment, 47.9% of data elements were captured through paper chart review, with the data warehouse query yielding 26.9% of data elements. Score application resulted in a prioritized subset of 17 criteria, which yielded a sensitivity of 1.00 and specificity 0.57 on EHR data (paper charts, 1.00 and 0.80) compared with actual recruitment in the trial., Conclusion: It is possible to prioritize clinical trial eligibility criteria based on real-world data to optimize prescreening of patients on a selected subset of relevant and available criteria and reduce implementation efforts for recruitment support. The performance could be further improved by increasing EHR data coverage., Competing Interests: T.G. reports grants from European Commission during the conduct of the study. G.M. reports grants from European Commission, during the conduct of the study. H.U.P. reports grants from European Commission, during the conduct of the study., (Thieme. All rights reserved.)

The induction of an interferon-mediated response is the first line of defense against pathogens such as viruses. Yet, the dynamics and extent of interferon alpha (IFNα)-induced antiviral genes vary remarkably and comprise three expression clusters: early, intermediate and late. By mathematical modeling based on time-resolved quantitative data, we identified mRNA stability as well as a negative regulatory loop as key mechanisms endogenously controlling the expression dynamics of IFNα-induced antiviral genes in hepatocytes. Guided by the mathematical model, we uncovered that this regulatory loop is mediated by the transcription factor IRF2 and showed that knock-down of IRF2 results in enhanced expression of early, intermediate and late IFNα-induced antiviral genes. Co-stimulation experiments with different pro-inflammatory cytokines revealed that this amplified expression dynamics of the early, intermediate and late IFNα-induced antiviral genes can also be achieved by co-application of IFNα and interleukin1 beta (IL1β). Consistently, we found that IL1β enhances IFNα-mediated repression of viral replication. Conversely, we observed that in IL1β receptor knock-out mice replication of viruses sensitive to IFNα is increased. Thus, IL1β is capable to potentiate IFNα-induced antiviral responses and could be exploited to improve antiviral therapies., Competing Interests: The authors have declared that no competing interests exist.

Tightly interlinked feedback regulators control the dynamics of intracellular responses elicited by the activation of signal transduction pathways. Interferon alpha (IFNα) orchestrates antiviral responses in hepatocytes, yet mechanisms that define pathway sensitization in response to prestimulation with different IFNα doses remained unresolved. We establish, based on quantitative measurements obtained for the hepatoma cell line Huh7.5, an ordinary differential equation model for IFNα signal transduction that comprises the feedback regulators STAT1, STAT2, IRF9, USP18, SOCS1, SOCS3, and IRF2. The model-based analysis shows that, mediated by the signaling proteins STAT2 and IRF9, prestimulation with a low IFNα dose hypersensitizes the pathway. In contrast, prestimulation with a high dose of IFNα leads to a dose-dependent desensitization, mediated by the negative regulators USP18 and SOCS1 that act at the receptor. The analysis of basal protein abundance in primary human hepatocytes reveals high heterogeneity in patient-specific amounts of STAT1, STAT2, IRF9, and USP18. The mathematical modeling approach shows that the basal amount of USP18 determines patient-specific pathway desensitization, while the abundance of STAT2 predicts the patient-specific IFNα signal response., (© 2020 The Authors. Published under the terms of the CC BY 4.0 license.)

Blindness caused by the eye diseases Retinitis-Pigmentosa and Age-Related-Macular-Degeneration leads to a degeneration of the photoreceptor layer while postsynaptic cells mostly stay intact. In this Paper a new concept for retinal implants is proposed. Instead of converting the incident light to a gray-scale picture with corresponding continuous-value stimulation levels, we here suggest to produce a binary image picture that only highlight edges in order to stimulate the retina solely at points which belong to an edge. An integrated test circuit is designed with a 130 nm BiCMOS process by using cellular neural networks for binary image generation. The circuit yields a simulated maximum rated power consumption of 2.61 mW for a 1000 information processing cells.

Activation of the Met receptor tyrosine kinase, either by its ligand, hepatocyte growth factor (HGF), or via ligand-independent mechanisms, such as MET amplification or receptor overexpression, has been implicated in driving tumor proliferation, metastasis, and resistance to therapy. Clinical development of Met-targeted antibodies has been challenging, however, as bivalent antibodies exhibit agonistic properties, whereas monovalent antibodies lack potency and the capacity to down-regulate Met. Through computational modeling, we found that the potency of a monovalent antibody targeting Met could be dramatically improved by introducing a second binding site that recognizes an unrelated, highly expressed antigen on the tumor cell surface. Guided by this prediction, we engineered MM-131, a bispecific antibody that is monovalent for both Met and epithelial cell adhesion molecule (EpCAM). MM-131 is a purely antagonistic antibody that blocks ligand-dependent and ligand-independent Met signaling by inhibiting HGF binding to Met and inducing receptor down-regulation. Together, these mechanisms lead to inhibition of proliferation in Met-driven cancer cells, inhibition of HGF-mediated cancer cell migration, and inhibition of tumor growth in HGF-dependent and -independent mouse xenograft models. Consistent with its design, MM-131 is more potent in EpCAM-high cells than in EpCAM-low cells, and its potency decreases when EpCAM levels are reduced by RNAi. Evaluation of Met, EpCAM, and HGF levels in human tumor samples reveals that EpCAM is expressed at high levels in a wide range of Met-positive tumor types, suggesting a broad opportunity for clinical development of MM-131., Competing Interests: Conflict of interest statement: Several authors are or were employees of Merrimack Pharmaceuticals, Inc. at the time of contributing to the manuscript.

In systems biology, one of the major tasks is to tailor model complexity to information content of the data. A useful model should describe the data and produce well-determined parameter estimates and predictions. Too small of a model will not be able to describe the data whereas a model which is too large tends to overfit measurement errors and does not provide precise predictions. Typically, the model is modified and tuned to fit the data, which often results in an oversized model. To restore the balance between model complexity and available measurements, either new data has to be gathered or the model has to be reduced. In this manuscript, we present a data-based method for reducing non-linear models. The profile likelihood is utilised to assess parameter identifiability and designate likely candidates for reduction. Parameter dependencies are analysed along profiles, providing context-dependent suggestions for the type of reduction. We discriminate four distinct scenarios, each associated with a specific model reduction strategy. Iterating the presented procedure eventually results in an identifiable model, which is capable of generating precise and testable predictions. Source code for all toy examples is provided within the freely available, open-source modelling environment Data2Dynamics based on MATLAB available at http://www.data2dynamics.org/, as well as the R packages dMod/cOde available at https://github.com/dkaschek/. Moreover, the concept is generally applicable and can readily be used with any software capable of calculating the profile likelihood., Competing Interests: The authors have declared that no competing interests exist. Author AR received support in the form of salaries from Merrimack Pharmaceuticals. This does not alter the authors’ adherence to PLOS ONE policies on sharing data and materials.

Signaling pathways are characterized by crosstalk, feedback and feedforward mechanisms giving rise to highly complex and cell-context specific signaling networks. Dissecting the underlying relations is crucial to predict the impact of targeted perturbations. However, a major challenge in identifying cell-context specific signaling networks is the enormous number of potentially possible interactions. Here, we report a novel hybrid mathematical modeling strategy to systematically unravel hepatocyte growth factor (HGF) stimulated phosphoinositide-3-kinase (PI3K) and mitogen activated protein kinase (MAPK) signaling, which critically contribute to liver regeneration. By combining time-resolved quantitative experimental data generated in primary mouse hepatocytes with interaction graph and ordinary differential equation modeling, we identify and experimentally validate a network structure that represents the experimental data best and indicates specific crosstalk mechanisms. Whereas the identified network is robust against single perturbations, combinatorial inhibition strategies are predicted that result in strong reduction of Akt and ERK activation. Thus, by capitalizing on the advantages of the two modeling approaches, we reduce the high combinatorial complexity and identify cell-context specific signaling networks.

Background: Mathematical models of dynamical systems facilitate the computation of characteristic properties that are not accessible experimentally. In cell biology, two main properties of interest are (1) the time-period a protein is accessible to other molecules in a certain state - its half-life - and (2) the time it spends when passing through a subsystem - its transit-time. We discuss two approaches to quantify the half-life, present the novel method of in silico labeling, and introduce the label half-life and label transit-time. The developed method has been motivated by laboratory tracer experiments. To investigate the kinetic properties and behavior of a substance of interest, we computationally label this species in order to track it throughout its life cycle. The corresponding mathematical model is extended by an additional set of reactions for the labeled species, avoiding any double-counting within closed circuits, correcting for the influences of upstream fluxes, and taking into account combinatorial multiplicity for complexes or reactions with several reactants or products. A profile likelihood approach is used to estimate confidence intervals on the label half-life and transit-time., Results: Application to the JAK-STAT signaling pathway in Epo-stimulated BaF3-EpoR cells enabled the calculation of the time-dependent label half-life and transit-time of STAT species. The results were robust against parameter uncertainties., Conclusions: Our approach renders possible the estimation of species and label half-lives and transit-times. It is applicable to large non-linear systems and an implementation is provided within the PottersWheel modeling framework (http://www.potterswheel.de).

The program PottersWheel has been developed to provide an intuitive and yet powerful framework for data-based modeling of dynamical systems like biochemical reaction networks. Its key functionality is multi-experiment fitting, where several experimental data sets from different laboratory conditions are fitted simultaneously in order to improve the estimation of unknown model parameters, to check the validity of a given model, and to discriminate competing model hypotheses. New experiments can be designed interactively. Models are either created text-based or using a visual model designer. Dynamically generated and compiled C files provide fast simulation and fitting procedures. Each function can either be accessed using a graphical user interface or via command line, allowing for batch processing within custom Matlab scripts. PottersWheel is designed as a Matlab toolbox, comprises 250,000 lines of Matlab and C code, and is freely available for academic usage at www.potterswheel.de .

Epidermal growth factor receptors (ErbB1-4) are oncogenic receptor tyrosine kinases (RTKs) that regulate diverse cellular processes. In this study, we combine measurement and mathematical modeling to quantify phospho-turnover at ErbB receptors in human cells and to determine the consequences for signaling and drug binding. We find that phosphotyrosine residues on ErbB1 have half-lives of a few seconds and therefore turn over 100-1000 times in the course of a typical immediate-early response to ligand. Rapid phospho-turnover is also observed for EGF-activated ErbB2 and ErbB3, unrelated RTKs, and multiple intracellular adaptor proteins and signaling kinases. Thus, the complexes formed on the cytoplasmic tail of active receptors and the downstream signaling kinases they control are highly dynamic and antagonized by potent phosphatases. We develop a kinetic scheme for binding of anti-ErbB1 drugs to receptors and show that rapid phospho-turnover significantly impacts their mechanisms of action., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Up to 29% of patients with epilepsy report "prodromal" sensations more than 30 minutes prior to seizures. We developed and implemented an objective methodology to prospectively assess the sensitivity and specificity of these subjective experiences using personal digital assistants (PDAs). The key property, in contrast to paper-based diaries, is the internal recording of the patient's entering time of prodromes and seizures. Of 500 patients with epilepsy interviewed, 31 claimed to sense prodromal symptoms at least 30 minutes before seizure onset. Eleven of them agreed to participate in a 4-week study to objectively measure their prospective prediction performance. In 9 patients returning data, the majority of prodrome entries were not followed by seizures or were identified only retrospectively. Statistical analysis revealed that no patient could outperform a nonspecific random predictor when predicting seizures based on the occurrence of prodromes, and that the group performance matched precisely the expected result for a by-chance prediction. These results question the predictive value of "prodromes" and the specificity of their occurrence in the preictal period., (Copyright © 2011. Published by Elsevier Inc.)