Biomedical subjects

Ursula Kummer

Publications and source records attributed to Ursula Kummer.

10 recordsLinked to original sources

PubMed2006

COPASI--a COmplex PAthway SImulator.

MOTIVATION: Simulation and modeling is becoming a standard approach to understand complex biochemical processes. Therefore, there is a big need for software tools that allow access to diverse simulation and modeling methods as well as support for the usage of these methods. RESULTS: Here, we present COPASI, a platform-independent and user-friendly biochemical simulator that offers several unique features. We discuss numerical issues with these features; in particular, the criteria to switch between stochastic and deterministic simulation methods, hybrid deterministic-stochastic methods, and the importance of random number generator numerical resolution in stochastic simulation. AVAILABILITY: The complete software is available in binary (executable) for MS Windows, OS X, Linux (Intel) and Sun Solaris (SPARC), as well as the full source code under an open source license from http://www.copasi.org.

Algorithms↗

PubMed2006

Quantification of cross correlations in complex spatiotemporal systems.

We propose a design of the equal time correlation matrix suitable for the analysis of multivariate time series with ill-defined phases. We present the cross-correlation analysis of model data sets taken from coupled stochastic oscillators and compare the concept with the results obtained from a conventional correlation matrix analysis. We show that the concept provides a higher sensitivity combined with a better statistical significance when quantifying weak cross correlations.

Journal Article↗

PubMed2005

A new dynamical layout algorithm for complex biochemical reaction networks.

BACKGROUND: To study complex biochemical reaction networks in living cells researchers more and more rely on databases and computational methods. In order to facilitate computational approaches, visualisation techniques are highly important. Biochemical reaction networks, e.g. metabolic pathways are often depicted as graphs and these graphs should be drawn dynamically to provide flexibility in the context of different data. Conventional layout algorithms are not sufficient for every kind of pathway in biochemical research. This is mainly due to certain conventions to which biochemists/biologists are used to and which are not in accordance to conventional layout algorithms. A number of approaches has been developed to improve this situation. Some of these are used in the context of biochemical databases and make more or less use of the information in these databases to aid the layout process. However, visualisation becomes also more and more important in modelling and simulation tools which mostly do not offer additional connections to databases. Therefore, layout algorithms used in these tools have to work independently of any databases. In addition, all of the existing algorithms face some limitations with respect to the number of edge crossings when it comes to larger biochemical systems due to the interconnectivity of these. Last but not least, in some cases, biochemical conventions are not met properly. RESULTS: Out of these reasons we have developed a new algorithm which tackles these problems by reducing the number of edge crossings in complex systems, taking further biological conventions into account to identify and visualise cycles. Furthermore the algorithm is independent from database information in order to be easily adopted in any application. It can also be tested as part of the SimWiz package (free to download for academic users at 1). CONCLUSION: The new algorithm reduces the complexity of pathways, as well as edge crossings and edge length in the resulting graphical representation. It also considers existing and further biological conventions to create a drawing most biochemists are familiar with. A lot of examples can be found on 2.

Algorithms↗

PubMed2005

Transition from stochastic to deterministic behavior in calcium oscillations.

Simulation and modeling is becoming more and more important when studying complex biochemical systems. Most often, ordinary differential equations are employed for this purpose. However, these are only applicable when the numbers of participating molecules in the biochemical systems are large enough to be treated as concentrations. For smaller systems, stochastic simulations on discrete particle basis are more accurate. Unfortunately, there are no general rules for determining which method should be employed for exactly which problem to get the most realistic result. Therefore, we study the transition from stochastic to deterministic behavior in a widely studied system, namely the signal transduction via calcium, especially calcium oscillations. We observe that the transition occurs within a range of particle numbers, which roughly corresponds to the number of receptors and channels in the cell, and depends heavily on the attractive properties of the phase space of the respective systems dynamics. We conclude that the attractive properties of a system, expressed, e.g., by the divergence of the system, are a good measure for determining which simulation algorithm is appropriate in terms of speed and realism.

Adenosine Triphosphate↗

PubMed2004

On the encoding and decoding of calcium signals in hepatocytes.

Many different agonists use calcium as a second messenger. Despite intensive research in intracellular calcium signalling it is an unsolved riddle how the different types of information represented by the different agonists, is encoded using the universal carrier calcium. It is also still not clear how the information encoded is decoded again into the intracellular specific information at the site of enzymes and genes. After the discovery of calcium oscillations, one likely mechanism is that information is encoded in the frequency, amplitude and waveform of the oscillations. This hypothesis has received some experimental support. However, the mechanism of decoding of oscillatory signals is still not known. Here, we study a mechanistic model of calcium oscillations, which is able to reproduce both spiking and bursting calcium oscillations. We use the model to study the decoding of calcium signals on the basis of co-operativity of calcium binding to various proteins. We show that this co-operativity offers a simple way to decode different calcium dynamics into different enzyme activities.

Animals↗

PubMed2004

CombAlign: a protein sequence comparison algorithm considering recombinations.

The basic linear treatment of sequence comparisons limits the ability of contemporary sequence alignment algorithms to detect non-order-conserving recombinations. Here, we introduce the algorithm combAlign which addresses the assessment of pairwise sequence similarity on non-order-conserving recombinations on a large scale. Emphasizing a two-level approach, combAlign first detects locally well conserved subsequences in a target and a source sequence. Subsequently, the relative placement of alignments is mapped to a graph. Concatenating local alignments to reassemble the target sequence to the fullest extent, the maximum scoring path through the graph denotes the best attainable combAlignment. Parameters influencing this process can be set to meet the user's specific demands. combAlign is applied to examples demonstrating the possibility to reflect evolutionary kinship of proteins even if their domains and motifs are strongly rearranged.

Algorithms↗

PubMed2003

Mechanism of melatonin-induced oscillations in the peroxidase-oxidase reaction.

Melatonin induces oscillations in the peroxidase-oxidase (PO) reaction catalyzed by horseradish peroxidase. We present here studies of the effect of pH, enzyme concentration, and concentration of melatonin on the oscillation frequency. We also present a mechanistic model to explain the experimentally observed changes in oscillation frequency. Using the data obtained here we are able to predict that oscillations will also occur in the PO reaction catalyzed by myeloperoxidase. Myeloperoxidase is an important protein in activated neutrophils and we provide evidence that the oscillations of NAD(P)H, superoxide and hydrogen peroxide in these cells may involve this enzyme. Thus, our experimental system can be considered a model system for the nonrespiratory oxygen metabolism in activated neutrophils and other similar cells participating in the defence against invading pathogens.

Dose-Response Relationship, Drug↗

PubMed2003

Concerted simulations reveal how peroxidase compound III formation results in cellular oscillations.

A major problem in mathematical modeling of the dynamics of complex biological systems is the frequent lack of knowledge of kinetic parameters. Here, we apply Brownian dynamics simulations, based on protein three-dimensional structures, to estimate a previously undetermined kinetic parameter, which is then used in biochemical network simulations. The peroxidase-oxidase reaction involves many elementary steps and displays oscillatory dynamics important for immune response. Brownian dynamics simulations were performed for three different peroxidases to estimate the rate constant for one of the elementary steps crucial for oscillations in the peroxidase-oxidase reaction, the association of superoxide with peroxidase. Computed second-order rate constants agree well with available experimental data and permit prediction of rate constants at physiological conditions. The simulations show that electrostatic interactions depress the rate of superoxide association with myeloperoxidase, bringing it into the range necessary for oscillatory behavior in activated neutrophils. Such negative electrostatic steering of enzyme-substrate association presents a novel control mechanism and lies in sharp contrast to the electrostatically-steered fast association of superoxide and Cu/Zn superoxide dismutase, which is also simulated here. The results demonstrate the potential of an integrated and concerted application of structure-based simulations and biochemical network simulations in cellular systems biology.

Animals↗

PubMed2003

A model of the oscillatory metabolism of activated neutrophils.

We present a two-compartment model to explain the oscillatory behavior observed experimentally in activated neutrophils. Our model is based mainly on the peroxidase-oxidase reaction catalyzed by myeloperoxidase with melatonin as a cofactor and NADPH oxidase, a major protein in the phagosome membrane of the leukocyte. The model predicts that after activation of a neutrophil, an increase in the activity of the hexose monophosphate shunt and the delivery of myeloperoxidase into the phagosome results in oscillations in oxygen and NAD(P)H concentration. The period of oscillation changes from >200 s to 10-30 s. The model is consistent with previously reported oscillations in cell metabolism and oxidant production. Key features and predictions of the model were confirmed experimentally. The requirement of the hexose monophosphate pathway for 10 s oscillations was verified using 6-aminonicotinamide and dexamethasone, which are inhibitors of glucose-6-phosphate dehydrogenase. The role of the NADPH oxidase in promoting oscillations was confirmed by dose-response studies of the effect of diphenylene iodonium, an inhibitor of the NADPH oxidase. Moreover, the model predicted an increase in the amplitude of NADPH oscillations in the presence of melatonin, which was confirmed experimentally. Successful computer modeling of complex chemical dynamics within cells and their chemical perturbation will enhance our ability to identify new antiinflammatory compounds.

Catalysis↗

PubMed2003

Mechanism of protection of peroxidase activity by oscillatory dynamics.

The peroxidase-oxidase reaction is known to involve reactive oxygen species as intermediates. These intermediates inactivate many types of biomolecules, including peroxidase itself. Previously, we have shown that oscillatory dynamics in the peroxidase-oxidase reaction seem to protect the enzyme from inactivation. It was suggested that this is due to a lower average concentration of reactive oxygen species in the oscillatory state compared to the steady state. Here, we studied the peroxidase-oxidase reaction with either 4-hydroxybenzoic acid or melatonin as cofactors. We show that the protective effect of oscillatory dynamics is present in both cases. We also found that the enzyme degradation depends on the concentration of the cofactor and on the pH of the reaction mixture. We simulated the oscillatory behaviour, including the oscillation/steady state bistability observed experimentally, using a detailed reaction scheme. The computational results confirm the hypothesis that protection is due to lower average concentrations of superoxide radical during oscillations. They also show that the shape of the oscillations changes with increasing cofactor concentration resulting in a further decrease in the average concentration of radicals. We therefore hypothesize that the protective effect of oscillatory dynamics is a general effect in this system.

Animals↗