Ronny Feuer

Properties of Boolean networks and methods for their tests

Transcriptional regulation networks are often modeled as Boolean networks. We discuss certain properties of Boolean functions (BFs), which are considered as important in such networks, namely, membership to the classes of unate or canalizing functions. Of further interest is the average sensitivity (AS) of functions. In this article, we discuss several algorithms to test the properties of interest. To test canalizing properties of functions, we apply spectral techniques, which can also be used to characterize the AS of functions as well as the influences of variables in unate BFs. Further, we provide and review upper and lower bounds on the AS of unate BFs based on the spectral representation. Finally, we apply these methods to a transcriptional regulation network of Escherichia coli, which controls central parts of the E. coli metabolism. We find that all functions are unate. Also the analysis of the AS of the network reveals an exceptional robustness against transient fluctuations of the binary variables.a

LEMming: A Linear Error Model to Normalize Parallel Quantitative Real-Time PCR (qPCR) Data as an Alternative to Reference Gene Based Methods

Background Gene expression analysis is an essential part of biological and medical investigations. Quantitative real-time PCR (qPCR) is characterized with excellent sensitivity, dynamic range, reproducibility and is still regarded to be the gold standard for quantifying transcripts abundance. Parallelization of qPCR such as by microfluidic Taqman Fluidigm Biomark Platform enables evaluation of multiple transcripts in samples treated under various conditions. Despite advanced technologies, correct evaluation of the measurements remains challenging. Most widely used methods for evaluating or calculating gene expression data include geNorm and ΔΔC t, respectively. They rely on one or several stable reference genes (RGs) for normalization, thus potentially causing biased results. We therefore applied multivariable regression with a tailored error model to overcome the necessity of stable RGs. Results We developed a RG independent data normalization approach based on a tailored linear error model for parallel qPCR data, called LEMming. It uses the assumption that the mean C t values within samples of similarly treated groups are equal. Performance of LEMming was evaluated in three data sets with different stability patterns of RGs and compared to the results of geNorm normalization. Data set 1 showed that both methods gave similar results if stable RGs are available. Data set 2 included RGs which are stable according to geNorm criteria, but became differentially expressed in normalized data evaluated by a t-test. geNorm-normalized data showed an effect of a shifted mean per gene per condition whereas LEMming-normalized data did not. Comparing the decrease of standard deviation from raw data to geNorm and to LEMming, the latter was superior. In data set 3 according to geNorm calculated average expression stability and pairwise variation, stable RGs were available, but t-tests of raw data contradicted this. Normalization with RGs resulted in distorted data contradicting literature, while LEMming normalized data did not. Conclusions If RGs are coexpressed but are not independent of the experimental conditions the stability criteria based on inter- and intragroup variation fail. The linear error model developed, LEMming, overcomes the dependency of using RGs for parallel qPCR measurements, besides resolving biases of both technical and biological nature in qPCR. However, to distinguish systematic errors per treated group from a global treatment effect an additional measurement is needed. Quantification of total cDNA content per sample helps to identify systematic errors.

High definition urethral pressure profilometry: Evaluating a novel microtip catheter

Urethral pressure profilometry (UPP) is used in the diagnosis of stress urinary incontinence (SUI). SUI is a significant medical, social, and economic problem, affecting about 12.5% of the population. A novel microtip catheter was developed for UPP featuring an inclination sensor and higher angular resolution compared to systems in clinical use today. Therewith, the location of each measured pressure sample can be determined and the spatial pressure distribution inside the urethra reconstructed. In order to assess the performance and plausibility of data from the microtip catheter, we compare it to data from a double balloon air charged system.

Development of an arm support system to improve ergonomics in laparoscopic surgery: study design and provisional results.

BackgroundLaparoscopic surgery (LS) induces physical stress to the surgeon that is associated with an increased prevalence of musculoskeletal pain and injury in the shoulder–neck region. The aim of this research project is to develop an arm support system (ASsyst) that reduces physical stress and is applicable to various laparoscopic interventions and operation room settings.MethodsA systematic approach to develop an ASsyst started in October 2012 consisting of five consecutive steps. In step 1, 14 laparoscopic interventions were observed using subjective and objective measures to determine key indicators for the conception of an ASsyst in LS. In step 2, an expert workshop was held to find and evaluate solutions to generate concepts for a support system based on the results of step 1 and general methods. During the third step, prototypes of ASsyst were tested in an experimental setting. Steps 4 and 5 are currently in process and include the final development of the ASsyst using the most promising concept for the evaluation during simulated LS.ResultsIncreased levels of physical stress were found in LS. Asymmetric strains were common. Three prototypes of ASsyst emerged from step 1 and 2. These prototypes were a cable construction with a noose for the lower arm, a support from below the elbow and a pneumatic vest supporting the upper arm. The experimental testing of these prototypes demonstrated reduced physical stress when compared to the unsupported environment. The support from below the elbow seemed to be the most practical in terms of implementation in various operation room settings and acceptance by surgeons. Step 4 and 5 are still in process.ConclusionsErgonomic problems have been identified in LS that could be addressed by an ASsyst. The concept of supporting the elbow from below has been found to be the most promising approach.

Modeling and Simulation of Optimal Resource Management during the Diurnal Cycle in Emiliania huxleyi by Genome-Scale Reconstruction and an Extended Flux Balance Analysis Approach

The coccolithophorid unicellular alga Emiliania huxleyi is known to form large blooms, which have a strong effect on the marine carbon cycle. As a photosynthetic organism, it is subjected to a circadian rhythm due to the changing light conditions throughout the day. For a better understanding of the metabolic processes under these periodically-changing environmental conditions, a genome-scale model based on a genome reconstruction of the E. huxleyi strain CCMP 1516 was created. It comprises 410 reactions and 363 metabolites. Biomass composition is variable based on the differentiation into functional biomass components and storage metabolites. The model is analyzed with a flux balance analysis approach called diurnal flux balance analysis (diuFBA) that was designed for organisms with a circadian rhythm. It allows storage metabolites to accumulate or be consumed over the diurnal cycle, while keeping the structure of a classical FBA problem. A feature of this approach is that the production and consumption of storage metabolites is not defined externally via the biomass composition, but the result of optimal resource management adapted to the diurnally-changing environmental conditions. The model in combination with this approach is able to simulate the variable biomass composition during the diurnal cycle in proximity to literature data.

Interleukin-1β enhances FasL-induced caspase-3/-7 activity without increasing apoptosis in primary mouse hepatocytes.

Sustained inflammation may increase the susceptibility of hepatocytes to apoptotic cell death and therefore exacerbate liver damage. Here we report that the pro-inflammatory cytokine IL-1β sensitizes primary murine hepatocytes to Fas ligand (FasL)-induced caspase-3/-7 activity. This process was dependent on JNK1/2 and the BH3-only proteins Bim and Bid. Mathematical modeling revealed that incubation of hepatocytes with IL-1β depleted the anti-apoptotic Bcl-2 protein pool and thus shifted hepatocytes to mitochondrial type II apoptosis following Fas activation. As a consequence, IL-1β and FasL treatment enhanced cytochrome c release. Surprisingly, despite increased caspase-3/-7 activation, FasL-induced cell death was reduced by IL-1β pre-treatment. This protective effect was independent of JNK1/2, Bim or Bid. Furthermore, elevated caspase-3/-7 activity upon IL-1β and FasL treatment did not result in enhanced PARP cleavage. The protective effect of IL-1β was seen after 3 h of pre-incubation, indicating an anti-apoptotic transcriptional response. Indeed, NF-κB DNA binding was increased in response to IL-1β plus FasL and gene-expression profiling of NF-κB regulated genes revealed a transcriptional and translational upregulation of the caspase-8 inhibitor A20. A mathematical model was developed to explain the contradictious occurrence of both increased caspase-3/-7 activity and elevated cell viability by including a heterogeneous distribution of Bcl-2 proteins and variations in Fas signaling resulting in different subpopulations of hepatocytes.

Transition of an Anaerobic Escherichia coli Culture to Aerobiosis: Balancing mRNA and Protein Levels in a Demand-Directed Dynamic Flux Balance Analysis

The facultative anaerobic bacterium Escherichia coli is frequently forced to adapt to changing environmental conditions. One important determinant for metabolism is the availability of oxygen allowing a more efficient metabolism. Especially in large scale bioreactors, the distribution of oxygen is inhomogeneous and individual cells encounter frequent changes. This might contribute to observed yield losses during process upscaling. Short-term gene expression data exist of an anaerobic E. coli batch culture shifting to aerobic conditions. The data reveal temporary upregulation of genes that are less efficient in terms of energy conservation than the genes predicted by conventional flux balance analyses. In this study, we provide evidence for a positive correlation between metabolic fluxes and gene expression. We then hypothesize that the more efficient enzymes are limited by their low expression, restricting flux through their reactions. We define a demand that triggers expression of the demanded enzymes that we explicitly include in our model. With these features we propose a method, demand-directed dynamic flux balance analysis, dddFBA, bringing together elements of several previously published methods. The introduction of additional flux constraints proportional to gene expression provoke a temporary demand for less efficient enzymes, which is in agreement with the transient upregulation of these genes observed in the data. In the proposed approach, the applied objective function of growth rate maximization together with the introduced constraints triggers expression of metabolically less efficient genes. This finding is one possible explanation for the yield losses observed in large scale bacterial cultivations where steady oxygen supply cannot be warranted.

Model-based analysis of an adaptive evolution experiment with Escherichia coli in a pyruvate limited continuous culture with glycerol

Bacterial strains that were genetically blocked in important metabolic pathways and grown under selective conditions underwent a process of adaptive evolution: certain pathways may have been deregulated and therefore allowed for the circumvention of the given block. A block of endogenous pyruvate synthesis from glycerol was realized by a knockout of pyruvate kinase and phosphoenolpyruvate carboxylase in E. coli. The resulting mutant strain was able to grow on a medium containing glycerol and lactate, which served as an exogenous pyruvate source. Heterologous expression of a pyruvate carboxylase gene from Corynebacterium glutamicum was used for anaplerosis of the TCA cycle. Selective conditions were controlled in a continuous culture with limited lactate feed and an excess of glycerol feed. After 200–300 generations pyruvate-prototrophic mutants were isolated. The genomic analysis of an evolved strain revealed that the genotypic basis for the regained pyruvate-prototrophy was not obvious. A constraint-based model of the metabolism was employed to compute all possible detours around the given metabolic block by solving a hierarchy of linear programming problems. The regulatory network was expected to be responsible for the adaptation process. Hence, a Boolean model of the transcription factor network was connected to the metabolic model. Our model analysis only showed a marginal impact of transcriptional control on the biomass yield on substrate which is a key variable in the selection process. In our experiment, microarray analysis confirmed that transcriptional control probably played a minor role in the deregulation of the alternative pathways for the circumvention of the block.

Signal processing in urodynamics: towards high definition urethral pressure profilometry

BackgroundUrethral pressure profilometry (UPP) is used in the diagnosis of stress urinary incontinence (SUI) which is a significant medical, social, and economic problem. Low spatial pressure resolution, common occurrence of artifacts, and uncertainties in data location limit the diagnostic value of UPP. To overcome these limitations, high definition urethral pressure profilometry (HD-UPP) combining enhanced UPP hardware and signal processing algorithms has been developed. In this work, we present the different signal processing steps in HD-UPP and show experimental results from female minipigs.MethodsWe use a special microtip catheter with high angular pressure resolution and an integrated inclination sensor. Signals from the catheter are filtered and time-correlated artifacts removed. A signal reconstruction algorithm processes pressure data into a detailed pressure image on the urethra’s inside. Finally, the pressure distribution on the urethra’s outside is calculated through deconvolution. A mathematical model of the urethra is contained in a point-spread-function (PSF) which is identified depending on geometric and material properties of the urethra. We additionally investigate the PSF’s frequency response to determine the relevant frequency band for pressure information on the urinary sphincter.ResultsExperimental pressure data are spatially located and processed into high resolution pressure images. Artifacts are successfully removed from data without blurring other details. The pressure distribution on the urethra’s outside is reconstructed and compared to the one on the inside. Finally, the pressure images are mapped onto the urethral geometry calculated from inclination and position data to provide an integrated image of pressure distribution, anatomical shape, and location.ConclusionsWith its advanced sensing capabilities, the novel microtip catheter collects an unprecedented amount of urethral pressure data. Through sequential signal processing steps, physicians are provided with detailed information on the pressure distribution in and around the urethra. Therefore, HD-UPP overcomes many current limitations of conventional UPP and offers the opportunity to evaluate urethral structures, especially the sphincter, in context of the correct anatomical location. This could enable the development of focal therapy approaches in the treatment of SUI.

Sampling lattice and signal reconstruction in urodynamics

Urethral pressure profilometry is a tool in the diagnosis of urinary incontinence. The pressure profile along the urethra is measured by a special catheter in order to asses the contraction strength of the sphincter muscle. However, the diagnostic value of pressure profilometry is limited. We seek to increase the diagnostic value by providing a detailed spatial reconstruction of the pressure profile on the surface of the urethra. Therefore, we investigate the sampling lattice (i.e., spatial sample distribution) and its properties with respect to signal reconstruction depending on the catheters sensor configuration and retraction motion. A sampling lattice with beneficial properties for reconstruction is generated through downsampling and a coordinate transformation of the original samples. We investigate its stability properties and present a reconstruction algorithm for the two-dimensional pressure distribution. Simulation results for signal reconstruction and stability properties conclude this work.