Local Uniqueness of Alignments with a Fixed Proportion of Gaps

We consider two independent random strings with i.i.d. characters and examine their optimal alignments containing a fixed proportion of gaps. We prove that when the proportion of gaps is small then with high probability optimal alignments differ only in a small number of places and are locally unique everywhere else. The result is somewhat surprising, as one might expect unrelated sequences to admit many near-optimal alignments, whereas for related sequences one might expect an optimal alignment that is unique in many places

An Upper Bound on the Convergence Rate of a Second Functional in Optimal Sequence Alignment

Consider finite sequences $X_{[1,n]}=X_1\dots X_n$ and $Y_{[1,n]}=Y_1\dots Y_n$ of length $n$, consisting of i.i.d.\ samples of random letters from a finite alphabet, and let $S$ and $T$ be chosen i.i.d.\ randomly from the unit ball in the space of symmetric scoring functions over this alphabet augmented by a gap symbol. We prove a probabilistic upper bound of linear order in $n^{0.75}$ for the deviation of the score relative to $T$ of optimal alignments with gaps of $X_{[1,n]}$ and $Y_{[1,n]}$ relative to $S$. It remains an open problem to prove a lower bound. Our result contributes to the understanding of the microstructure of optimal alignments relative to one given scoring function, extending a theory begun by the first two authors

A Monte Carlo Approach to the Fluctuation Problem in Optimal Alignments of Random Strings

The problem of determining the correct order of fluctuation of the optimal alignment score of two random strings of length $n$ has been open for several decades. It is known [12] that the biased expected effect of a random letter-change on the optimal score implies an order of fluctuation linear in √$n$. However, in many situations where such a biased effect is observed empirically, it has been impossible to prove analytically. The main result of this paper shows that when the rescaled-limit of the optimal alignment score increases in a certain direction, then the biased effect exists. On the basis of this result one can quantify a confidence level for the existence of such a biased effect and hence of an order √$n$ fluctuation based on simulation of optimal alignments scores. This is an important step forward, as the correct order of fluctuation was previously known only for certain special distributions [12],[13],[5],[10]. To illustrate the usefulness of our new methodology, we apply it to optimal alignments of strings written in the DNA-alphabet. As scoring function, we use the BLASTZ default-substitution matrix together with a realistic gap penalty. BLASTZ is one of the most widely used sequence alignment methodologies in bioinformatics. For this DNA-setting, we show that with a high level of confidence, the fluctuation of the optimal alignment score is of order Θ(√$n$). An important special case of optimal alignment score is the Longest Common Subsequence (LCS) of random strings. For binary sequences with equiprobably symbols the question of the fluctuation of the LCS remains open. The symmetry in that case does not allow for our method. On the other hand, in real-life DNA sequences, it is not the case that all letters occur with the same frequency. So, for many real life situations, our method allows to determine the order of the fluctuation up to a high confidence level

Large deviation based upper bounds for the LCS-problem

We analyse and apply a large deviation and Montecarlo simulation based method for the computation of improved upper bounds on the Chvatal-Sankoff constant for i.i.d. random sequences over a finite alphabet. Our theoretical results show that this method converges to the exact value of when a control parameter converges to infinity. We also give upper bounds on the complexity for numerically computing the Chvatal-Sankoff constant to any given precision via this method. Our numerical experiments confirm the theory and allow us to give new upper bounds that are correct to two digits

Model evaluation for glycolytic oscillations in yeast biotransformations of xenobiotics

Anaerobic glycolysis in yeast perturbed by the reduction of xenobiotic ketones is studied numerically in two models which possess the same topology but different levels of complexity. By comparing both models' predictions for concentrations and fluxes as well as steady or oscillatory temporal behavior we answer the question what phenomena require what kind of minimum model abstraction. While mean concentrations and fluxes are predicted in agreement by both models we observe different domains of oscillatory behavior in parameter space. Generic properties of the glycolytic response to ketones are discussed

Measurement of the Associated γ+μ±\gamma + \mu^\pm Production Cross Section in ppˉp \bar p Collisions at s=1.8\sqrt{s} = 1.8 TeV

We present the first measurement of associated direct photon + muon production in hadronic collisions, from a sample of 1.8 TeV $p \bar p$ collisions recorded with the Collider Detector at Fermilab. Quantum chromodynamics (QCD) predicts that these events are primarily from the Compton scattering process $cg \to c\gamma$, with the final state charm quark producing a muon. Hence this measurement is sensitive to the charm quark content of the proton. The measured cross section of $29\pm 9 pb^{-1}$ is compared to a leading-order QCD parton shower model as well as a next-to-leading-order QCD calculation.Comment: 12 pages, 4 figures Added more detailed description of muon background estimat

Search for the Supersymmetric Partner of the Top-Quark in ppˉp \bar{p} Collisions at s=1.8TeV\sqrt{s} = 1.8 {\rm TeV}

We report on a search for the supersymmetric partner of the top quark (stop) produced in $t \bar{t}$ events using $110 {\rm pb}^{-1}$ of $p \bar{p}$ collisions at $\sqrt{s} = 1.8 {\rm TeV}$ recorded with the Collider Detector at Fermilab. In the case of a light stop squark, the decay of the top quark into stop plus the lightest supersymmetric particle (LSP) could have a significant branching ratio. The observed events are consistent with Standard Model $t \bar{t}$ production and decay. Hence, we set limits on the branching ratio of the top quark decaying into stop plus LSP, excluding branching ratios above 45% for a LSP mass up to 40 {\rm GeV/c}$^{2}$.Comment: 11 pages, 4 figure

Search for Second-Generation Scalar Leptoquarks in ppˉ\bm{p \bar{p}} Collisions at s\sqrt{s}=1.96 TeV

Results on a search for pair production of second generation scalar leptoquark in $p \bar{p}$ collisions at $\sqrt{s}$=1.96 TeV are reported. The data analyzed were collected by the CDF detector during the 2002-2003 Tevatron Run II and correspond to an integrated luminosity of 198 pb$^{-1}$. Leptoquarks (LQ) are sought through their decay into (charged) leptons and quarks, with final state signatures represented by two muons and jets and one muon, large transverse missing energy and jets. We observe no evidence for $LQ$ production and derive 95% C.L. upper limits on the $LQ$ production cross sections as well as lower limits on their mass as a function of $\beta$, where $\beta$ is the branching fraction for $LQ \to \mu q$.Comment: 9 pages (3 author list) 5 figure

Search for squarks and gluinos with the ATLAS detector in final states with jets and missing transverse momentum using √s=8 TeV proton-proton collision data

A search for squarks and gluinos in final states containing high-p T jets, missing transverse momentum and no electrons or muons is presented. The data were recorded in 2012 by the ATLAS experiment in s√=8 TeV proton-proton collisions at the Large Hadron Collider, with a total integrated luminosity of 20.3 fb−1. Results are interpreted in a variety of simplified and specific supersymmetry-breaking models assuming that R-parity is conserved and that the lightest neutralino is the lightest supersymmetric particle. An exclusion limit at the 95% confidence level on the mass of the gluino is set at 1330 GeV for a simplified model incorporating only a gluino and the lightest neutralino. For a simplified model involving the strong production of first- and second-generation squarks, squark masses below 850 GeV (440 GeV) are excluded for a massless lightest neutralino, assuming mass degenerate (single light-flavour) squarks. In mSUGRA/CMSSM models with tan β = 30, A 0 = −2m 0 and μ > 0, squarks and gluinos of equal mass are excluded for masses below 1700 GeV. Additional limits are set for non-universal Higgs mass models with gaugino mediation and for simplified models involving the pair production of gluinos, each decaying to a top squark and a top quark, with the top squark decaying to a charm quark and a neutralino. These limits extend the region of supersymmetric parameter space excluded by previous searches with the ATLAS detector