Tailored parameter optimization methods for ordinary differential equation models with steady-state constraints

Background: Ordinary differential equation (ODE) models are widely used to describe (bio-)chemical and biological processes. To enhance the predictive power of these models, their unknown parameters are estimated from experimental data. These experimental data are mostly collected in perturbation experiments, in which the processes are pushed out of steady state by applying a stimulus. The information that the initial condition is a steady state of the unperturbed process provides valuable information, as it restricts the dynamics of the process and thereby the parameters. However, implementing steady-state constraints in the optimization often results in convergence problems. Results: In this manuscript, we propose two new methods for solving optimization problems with steady-state constraints. The first method exploits ideas from optimization algorithms on manifolds and introduces a retraction operator, essentially reducing the dimension of the optimization problem. The second method is based on the continuous analogue of the optimization problem. This continuous analogue is an ODE whose equilibrium points are the optima of the constrained optimization problem. This equivalence enables the use of adaptive numerical methods for solving optimization problems with steady-state constraints. Both methods are tailored to the problem structure and exploit the local geometry of the steady-state manifold and its stability properties. A parameterization of the steady-state manifold is not required. The efficiency and reliability of the proposed methods is evaluated using one toy example and two applications. The first application example uses published data while the second uses a novel dataset for Raf/MEK/ERK signaling. The proposed methods demonstrated better convergence properties than state-of-the-art methods employed in systems and computational biology. Furthermore, the average computation time per converged start is significantly lower. In addition to the theoretical results, the analysis of the dataset for Raf/MEK/ERK signaling provides novel biological insights regarding the existence of feedback regulation. Conclusion: Many optimization problems considered in systems and computational biology are subject to steady-state constraints. While most optimization methods have convergence problems if these steady-state constraints are highly nonlinear, the methods presented recover the convergence properties of optimizers which can exploit an analytical expression for the parameter-dependent steady state. This renders them an excellent alternative to methods which are currently employed in systems and computational biology

Toward system innovation for more sustainable chemistry: insights into consumers’ perceptions, knowledge, and behavior related to traceability and product design strategies along leather supply chains

The leather industry is a complex system with multiple actors that faces a fundamental transition toward more sustainable chemistry. To support this process, this article analyzes challenges of the industry and consumers’ roles as a nexus of transition-relevant developments. We present findings of an empirical study (N = 439) among consumers on their perception of leather, related knowledge, and purchasing behavior. We found that participants perceived leather as natural, robust, and of high quality. Knowledge about the manufacturing of leather products was overall limited but varied. Applying a psychological behavior theory, we found that being aware of environmental and health consequences from conventional manufacturing of leather products was positively associated with a personal norm to purchase leather products that are less harmful to environment and health. The perceived ease of buying such products was positively associated with their purchase. Our findings shed light on consumers’ roles in the current leather system and their support of niche innovations toward more sustainable chemistry. Against this backdrop, we discuss implications for product design, consumer information, and needs for traceability along supply chains

EpiScanpy: integrated single-cell epigenomic analysis

EpiScanpy is a toolkit for the analysis of single-cell epigenomic data, namely single-cell DNA methylation and single-cell ATAC-seq data. To address the modality specific challenges from epigenomics data, epiScanpy quantifies the epigenome using multiple feature space constructions and builds a nearest neighbour graph using epigenomic distance between cells. EpiScanpy makes the many existing scRNA-seq workflows from scanpy available to large-scale single-cell data from other -omics modalities, including methods for common clustering, dimension reduction, cell type identification and trajectory learning techniques, as well as an atlas integration tool for scATAC-seq datasets. The toolkit also features numerous useful downstream functions, such as differential methylation and differential openness calling, mapping epigenomic features of interest to their nearest gene, or constructing gene activity matrices using chromatin openness. We successfully benchmark epiScanpy against other scATAC-seq analysis tools and show its outperformance at discriminating cell types. The authors present epiScanpy: a computational framework for the analysis of single-cell epigenomic data, both ATAC-seq and DNA methylation data, with examples for clustering, cell type identification, trajectory learning and atlas integration - and show its performance in distinguishing cell types

Inhibition of Aromatase Induces Partial Sex Change in a Cichlid Fish: Distinct Functions for Sex Steroids in Brains and Gonads

Sex steroids are major drivers of sexual development and also responsible for the maintenance of the established gender. Especially fishes exhibit great plasticity and less conservation in sex determination and sexual development compared to other vertebrate groups. In addition, fishes have a constant sex steroid production throughout their entire lifespan, which makes them particularly susceptible to interferences with the endogenous sex steroid system. This susceptibility has recently been used to show that inhibition of the key enzyme of estrogen synthesis, aromatase Cyp19a1, can induce functional sex reversal even in adult fish. Here, we investigated the impact of the aromatase inhibitor (AI) fadrozole in adult females of the East African cichlid fish Astatotilapia burtoni. Using gene expression, phenotypic measurements, behavioral experiments, and hormone measurements, we assessed if females treated with fadrozole develop a male-like phenotype. We found that AI treatment has a different effect on gene expression in the gonad compared to the brain, the 2 tissues mostly implicated in sexual development. In contrast to observations in other gonochoristic species, A. burtoni ovaries cannot be transformed into functional testis by AI. However, rapid changes towards a male-like phenotype can be induced with AI in coloration, hormone levels, and behavior

Inhibitory framing in hypersexual patients with Parkinson's disease. An fMRI pilot study

Hypersexuality in medicated patients with PD is caused by an increased influence of motivational drive areas and a decreased influence of inhibitory control areas due to dopaminergic medication. In this pilot study, we test a newly developed paradigm investigating the influence of dopaminergic medication on brain activation elicited by sexual pictures with and without inhibitory contextual framing. Twenty PD patients with and without hypersexuality were examined with fMRI either OFF or ON standardized dopaminergic medication. The paradigm consisted of a priming phase where either a neutral context or an inhibitory context was presented. This priming phase was either followed by a sexual or a neutral target. Sexual, compared to neutral pictures resulted in a BOLD activation of various brain regions implicated in sexual processing. Hypersexual PD patients showed increased activity compared to PD controls in these regions. There was no relevant effect of medication between the two groups. The inhibitory context elicited less activation in inhibition-related areas in hypersexual PD, but had no influence on the perception of sexual cues. The paradigm partially worked: reactivity of motivational brain areas to sexual cues was increased in hypersexual PD and inhibitory contextual framing lead to decreased activation of inhibitory control areas in PD. We could not find a medication effect and the length of the inhibitory stimulus was not optimal to suppress reactivity to sexual cues. Our data provide new insights into the mechanisms of hypersexuality and warrant a replication with a greater cohort and an optimized stimulus length in the future

Intra- and inter-rater reliability of joint range of motion tests using tape measure, digital inclinometer and inertial motion capturing

Background In clinical practice range of motion (RoM) is usually assessed with low-cost devices such as a tape measure (TM) or a digital inclinometer (DI). However, the intra- and inter-rater reliability of typical RoM tests differ, which impairs the evaluation of therapy progress. More objective and reliable kinematic data can be obtained with the inertial motion capture system (IMC) by Xsens. The aim of this study was to obtain the intra- and inter-rater reliability of the TM, DI and IMC methods in five RoM tests: modified Thomas test (DI), shoulder test modified after Janda (DI), retroflexion of the trunk modified after Janda (DI), lateral inclination (TM) and fingertip-to-floor test (TM). Methods Two raters executed the RoM tests (TM or DI) in a randomized order on 22 healthy individuals while, simultaneously, the IMC data (Xsens MVN) was collected. After 15 warm-up repetitions, each rater recorded five measurements. Findings Intra-rater reliabilities were (almost) perfect for tests in all three devices (ICCs 0.886–0.996). Inter-rater reliability was substantial to (almost) perfect in the DI (ICCs 0.71–0.87) and the IMC methods (ICCs 0.61–0.993) and (almost) perfect in the TM methods (ICCs 0.923–0.961). The measurement error (ME) for the tests measured in degree (°) was 0.9–3.3° for the DI methods and 0.5–1.2° for the IMC approaches. In the tests measured in centimeters the ME was 0.5–1.3cm for the TM methods and 0.6–2.7cm for the IMC methods. Pearson correlations between the results of the DI or the TM respectively with the IMC results were significant in all tests except for the shoulder test on the right body side (r = 0.41–0.81). Interpretation Measurement repetitions of either one or multiple trained raters can be considered reliable in all three devices