Modularization of biochemical networks based on classification of Petri net t-invariants

<p>Abstract</p> <p>Background</p> <p>Structural analysis of biochemical networks is a growing field in bioinformatics and systems biology. The availability of an increasing amount of biological data from molecular biological networks promises a deeper understanding but confronts researchers with the problem of combinatorial explosion. The amount of qualitative network data is growing much faster than the amount of quantitative data, such as enzyme kinetics. In many cases it is even impossible to measure quantitative data because of limitations of experimental methods, or for ethical reasons. Thus, a huge amount of qualitative data, such as interaction data, is available, but it was not sufficiently used for modeling purposes, until now. New approaches have been developed, but the complexity of data often limits the application of many of the methods. Biochemical Petri nets make it possible to explore static and dynamic qualitative system properties. One Petri net approach is model validation based on the computation of the system's invariant properties, focusing on t-invariants. T-invariants correspond to subnetworks, which describe the basic system behavior.</p> <p>With increasing system complexity, the basic behavior can only be expressed by a huge number of t-invariants. According to our validation criteria for biochemical Petri nets, the necessary verification of the biological meaning, by interpreting each subnetwork (t-invariant) manually, is not possible anymore. Thus, an automated, biologically meaningful classification would be helpful in analyzing t-invariants, and supporting the understanding of the basic behavior of the considered biological system.</p> <p>Methods</p> <p>Here, we introduce a new approach to automatically classify t-invariants to cope with network complexity. We apply clustering techniques such as UPGMA, Complete Linkage, Single Linkage, and Neighbor Joining in combination with different distance measures to get biologically meaningful clusters (t-clusters), which can be interpreted as modules. To find the optimal number of t-clusters to consider for interpretation, the cluster validity measure, Silhouette Width, is applied.</p> <p>Results</p> <p>We considered two different case studies as examples: a small signal transduction pathway (pheromone response pathway in <it>Saccharomyces cerevisiae</it>) and a medium-sized gene regulatory network (gene regulation of Duchenne muscular dystrophy). We automatically classified the t-invariants into functionally distinct t-clusters, which could be interpreted biologically as functional modules in the network. We found differences in the suitability of the various distance measures as well as the clustering methods. In terms of a biologically meaningful classification of t-invariants, the best results are obtained using the Tanimoto distance measure. Considering clustering methods, the obtained results suggest that UPGMA and Complete Linkage are suitable for clustering t-invariants with respect to the biological interpretability.</p> <p>Conclusion</p> <p>We propose a new approach for the biological classification of Petri net t-invariants based on cluster analysis. Due to the biologically meaningful data reduction and structuring of network processes, large sets of t-invariants can be evaluated, allowing for model validation of qualitative biochemical Petri nets. This approach can also be applied to elementary mode analysis.</p

Search for a Technicolor omega_T Particle in Events with a Photon and a b-quark Jet at CDF

If the Technicolor omega_T particle exists, a likely decay mode is omega_T -> gamma pi_T, followed by pi_T -> bb-bar, yielding the signature gamma bb-bar. We have searched 85 pb^-1 of data collected by the CDF experiment at the Fermilab Tevatron for events with a photon and two jets, where one of the jets must contain a secondary vertex implying the presence of a b quark. We find no excess of events above standard model expectations. We express the result of an exclusion region in the M_omega_T - M_pi_T mass plane.Comment: 14 pages, 2 figures. Available from the CDF server (PS with figs): http://www-cdf.fnal.gov/physics/pub98/cdf4674_omega_t_prl_4.ps FERMILAB-PUB-98/321-

Measurement of the B0 anti-B0 oscillation frequency using l- D*+ pairs and lepton flavor tags

The oscillation frequency Delta-md of B0 anti-B0 mixing is measured using the partially reconstructed semileptonic decay anti-B0 -> l- nubar D*+ X. The data sample was collected with the CDF detector at the Fermilab Tevatron collider during 1992 - 1995 by triggering on the existence of two lepton candidates in an event, and corresponds to about 110 pb-1 of pbar p collisions at sqrt(s) = 1.8 TeV. We estimate the proper decay time of the anti-B0 meson from the measured decay length and reconstructed momentum of the l- D*+ system. The charge of the lepton in the final state identifies the flavor of the anti-B0 meson at its decay. The second lepton in the event is used to infer the flavor of the anti-B0 meson at production. We measure the oscillation frequency to be Delta-md = 0.516 +/- 0.099 +0.029 -0.035 ps-1, where the first uncertainty is statistical and the second is systematic.Comment: 30 pages, 7 figures. Submitted to Physical Review

Search for New Particles Decaying to top-antitop in proton-antiproton collisions at squareroot(s)=1.8 TeV

We use 106 \ipb of data collected with the Collider Detector at Fermilab to search for narrow-width, vector particles decaying to a top and an anti-top quark. Model independent upper limits on the cross section for narrow, vector resonances decaying to \ttbar are presented. At the 95% confidence level, we exclude the existence of a leptophobic \zpr boson in a model of topcolor-assisted technicolor with mass M_{\zpr} $<$ 480 \gev for natural width $\Gamma$ = 0.012 M_{\zpr}, and M_{\zpr} $<$ 780 \gev for $\Gamma$ = 0.04 M_{\zpr}.Comment: The CDF Collaboration, submitted to PRL 25-Feb-200

A Measurement of the Differential Dijet Mass Cross Section in p-pbar Collisions at sqrt{s}=1.8 TeV

We present a measurement of the cross section for production of two or more jets as a function of dijet mass, based on an integrated luminosity of 86 pb^-1 collected with the Collider Detector at Fermilab. Our dijet mass spectrum is described within errors by next-to-leading order QCD predictions using CTEQ4HJ parton distributions, and is in good agreement with a similar measurement from the D0 experiment.Comment: 18 pages including 2 figures and 3 tables. Submitted to Phys. Rev. D Rapid Communication

Diffractive Dijet Production at sqrt(s)=630 and 1800 GeV at the Fermilab Tevatron

We report a measurement of the diffractive structure function $F_{jj}^D$ of the antiproton obtained from a study of dijet events produced in association with a leading antiproton in $\bar pp$ collisions at $\sqrt s=630$ GeV at the Fermilab Tevatron. The ratio of $F_{jj}^D$ at $\sqrt s=630$ GeV to $F_{jj}^D$ obtained from a similar measurement at $\sqrt s=1800$ GeV is compared with expectations from QCD factorization and with theoretical predictions. We also report a measurement of the $\xi$ ($x$-Pomeron) and $\beta$ ($x$ of parton in Pomeron) dependence of $F_{jj}^D$ at $\sqrt s=1800$ GeV. In the region $0.035<\xi<0.095$, $|t|<1$ GeV$^2$ and $\beta<0.5$, $F_{jj}^D(\beta,\xi)$ is found to be of the form $\beta^{-1.0\pm 0.1} \xi^{-0.9\pm 0.1}$, which obeys $\beta$-$\xi$ factorization.Comment: LaTeX, 9 pages, Submitted to Phys. Rev. Letter

Search for Kaluza-Klein Graviton Emission in ppˉp\bar{p} Collisions at s=1.8\sqrt{s}=1.8 TeV using the Missing Energy Signature

We report on a search for direct Kaluza-Klein graviton production in a data sample of 84 ${pb}^{-1}$ of \ppb collisions at $\sqrt{s}$ = 1.8 TeV, recorded by the Collider Detector at Fermilab. We investigate the final state of large missing transverse energy and one or two high energy jets. We compare the data with the predictions from a $3+1+n$-dimensional Kaluza-Klein scenario in which gravity becomes strong at the TeV scale. At 95% confidence level (C.L.) for $n$=2, 4, and 6 we exclude an effective Planck scale below 1.0, 0.77, and 0.71 TeV, respectively.Comment: Submitted to PRL, 7 pages 4 figures/Revision includes 5 figure

Measurement of the average time-integrated mixing probability of b-flavored hadrons produced at the Tevatron

We have measured the number of like-sign (LS) and opposite-sign (OS) lepton pairs arising from double semileptonic decays of $b$ and $\bar{b}$-hadrons, pair-produced at the Fermilab Tevatron collider. The data samples were collected with the Collider Detector at Fermilab (CDF) during the 1992-1995 collider run by triggering on the existence of $\mu \mu$ and $e \mu$ candidates in an event. The observed ratio of LS to OS dileptons leads to a measurement of the average time-integrated mixing probability of all produced $b$-flavored hadrons which decay weakly, $\bar{\chi} = 0.152 \pm 0.007$ (stat.) $\pm 0.011$ (syst.), that is significantly larger than the world average $\bar{\chi} = 0.118 \pm 0.005$.Comment: 47 pages, 10 figures, 15 tables Submitted to Phys. Rev.

A Study of B0 -> J/psi K(*)0 pi+ pi- Decays with the Collider Detector at Fermilab

We report a study of the decays B0 -> J/psi K(*)0 pi+ pi-, which involve the creation of a u u-bar or d d-bar quark pair in addition to a b-bar -> c-bar(c s-bar) decay. The data sample consists of 110 1/pb of p p-bar collisions at sqrt{s} = 1.8 TeV collected by the CDF detector at the Fermilab Tevatron collider during 1992-1995. We measure the branching ratios to be BR(B0 -> J/psi K*0 pi+ pi-) = (8.0 +- 2.2 +- 1.5) * 10^{-4} and BR(B0 -> J/psi K0 pi+ pi-) = (1.1 +- 0.4 +- 0.2) * 10^{-3}. Contributions to these decays are seen from psi(2S) K(*)0, J/psi K0 rho0, J/psi K*+ pi-, and J/psi K1(1270)