Stochastic Simulations of the Repressilator Circuit

The genetic repressilator circuit consists of three transcription factors, or repressors, which negatively regulate each other in a cyclic manner. This circuit was synthetically constructed on plasmids in {\it Escherichia coli} and was found to exhibit oscillations in the concentrations of the three repressors. Since the repressors and their binding sites often appear in low copy numbers, the oscillations are noisy and irregular. Therefore, the repressilator circuit cannot be fully analyzed using deterministic methods such as rate-equations. Here we perform stochastic analysis of the repressilator circuit using the master equation and Monte Carlo simulations. It is found that fluctuations modify the range of conditions in which oscillations appear as well as their amplitude and period, compared to the deterministic equations. The deterministic and stochastic approaches coincide only in the limit in which all the relevant components, including free proteins, plasmids and bound proteins, appear in high copy numbers. We also find that subtle features such as cooperative binding and bound-repressor degradation strongly affect the existence and properties of the oscillations.Comment: Accepted to PR

Endothelin receptor B antagonists decrease glioma cell viability independently of their cognate receptor

Background: Endothelin receptor antagonists inhibit the progression of many cancers, but research into their influence on glioma has been limited. Methods: We treated glioma cell lines, LN-229 and SW1088, and melanoma cell lines, A375 and WM35, with two endothelin receptor type B (ETRB)-specific antagonists, A-192621 and BQ788, and quantified viable cells by the capacity of their intracellular esterases to convert non-fluorescent calcein AM into green-fluorescent calcein. We assessed cell proliferation by labeling cells with carboxyfluorescein diacetate succinimidyl ester and quantifying the fluorescence by FACS analysis. We also examined the cell cycle status using BrdU/propidium iodide double staining and FACS analysis. We evaluated changes in gene expression by microarray analysis following treatment with A-192621 in glioma cells. We examined the role of ETRB by reducing its expression level using small interfering RNA (siRNA). Results: We report that two ETRB-specific antagonists, A-192621 and BQ788, reduce the number of viable cells in two glioma cell lines in a dose- and time-dependent manner. We describe similar results for two melanoma cell lines. The more potent of the two antagonists, A-192621, decreases the mean number of cell divisions at least in part by inducing a G2/M arrest and apoptosis. Microarray analysis of the effects of A-192621 treatment reveals up-regulation of several DNA damage-inducible genes. These results were confirmed by real-time RT-PCR. Importantly, reducing expression of ETRB with siRNAs does not abrogate the effects of either A-192621 or BQ788 in glioma or melanoma cells. Furthermore, BQ123, an endothelin receptor type A (ETRA)-specific antagonist, has no effect on cell viability in any of these cell lines, indicating that the ETRB-independent effects on cell viability exhibited by A-192621 and BQ788 are not a result of ETRA inhibition. Conclusion: While ETRB antagonists reduce the viability of glioma cells in vitro, it appears unlikely that this effect is mediated by ETRB inhibition or cross-reaction with ETRA. Instead, we present evidence that A-192621 affects glioma and melanoma viability by activating stress/DNA damage response pathways, which leads to cell cycle arrest and apoptosis. This is the first evidence linking ETRB antagonist treatment to enhanced expression of DNA damage-inducible genes

Deterministic and stochastic descriptions of gene expression dynamics

A key goal of systems biology is the predictive mathematical description of gene regulatory circuits. Different approaches are used such as deterministic and stochastic models, models that describe cell growth and division explicitly or implicitly etc. Here we consider simple systems of unregulated (constitutive) gene expression and compare different mathematical descriptions systematically to obtain insight into the errors that are introduced by various common approximations such as describing cell growth and division by an effective protein degradation term. In particular, we show that the population average of protein content of a cell exhibits a subtle dependence on the dynamics of growth and division, the specific model for volume growth and the age structure of the population. Nevertheless, the error made by models with implicit cell growth and division is quite small. Furthermore, we compare various models that are partially stochastic to investigate the impact of different sources of (intrinsic) noise. This comparison indicates that different sources of noise (protein synthesis, partitioning in cell division) contribute comparable amounts of noise if protein synthesis is not or only weakly bursty. If protein synthesis is very bursty, the burstiness is the dominant noise source, independent of other details of the model. Finally, we discuss two sources of extrinsic noise: cell-to-cell variations in protein content due to cells being at different stages in the division cycles, which we show to be small (for the protein concentration and, surprisingly, also for the protein copy number per cell) and fluctuations in the growth rate, which can have a significant impact.Comment: 23 pages, 5 figures; Journal of Statistical physics (2012

The Persistence Length of a Strongly Charged, Rod-like, Polyelectrolyte in the Presence of Salt

The persistence length of a single, intrinsically rigid polyelectrolyte chain, above the Manning condensation threshold is investigated theoretically in presence of added salt. Using a loop expansion method, the partition function is consistently calculated, taking into account corrections to mean-field theory. Within a mean-field approximation, the well-known results of Odijk, Skolnick and Fixman are reproduced. Beyond mean-field, it is found that density correlations between counterions and thermal fluctuations reduce the stiffness of the chain, indicating an effective attraction between monomers for highly charged chains and multivalent counterions. This attraction results in a possible mechanical instability (collapse), alluding to the phenomenon of DNA condensation. In addition, we find that more counterions condense on slightly bent conformations of the chain than predicted by the Manning model for the case of an infinite cylinder. Finally, our results are compared with previous models and experiments.Comment: 13 pages, 2 ps figure

Measurement of the cross-section ratio sigma_{psi(2S)}/sigma_{J/psi(1S)} in deep inelastic exclusive ep scattering at HERA

The exclusive deep inelastic electroproduction of $\psi(2S)$ and $J/\psi(1S)$ at an $ep$ centre-of-mass energy of 317 GeV has been studied with the ZEUS detector at HERA in the kinematic range $2 < Q^2 < 80$ GeV$^2$, $30 < W < 210$ GeV and $|t| < 1$ GeV$^2$, where $Q^2$ is the photon virtuality, $W$ is the photon-proton centre-of-mass energy and $t$ is the squared four-momentum transfer at the proton vertex. The data for $2 < Q^2 < 5$ GeV$^2$ were taken in the HERA I running period and correspond to an integrated luminosity of 114 pb$^{-1}$. The data for $5 < Q^2 < 80$ GeV$^2$ are from both HERA I and HERA II periods and correspond to an integrated luminosity of 468 pb$^{-1}$. The decay modes analysed were $\mu^+\mu^-$ and $J/\psi(1S) \,\pi^+\pi^-$ for the $\psi(2S)$ and $\mu^+\mu^-$ for the $J/\psi(1S)$. The cross-section ratio $\sigma_{\psi(2S)}/\sigma_{J/\psi(1S)}$ has been measured as a function of $Q^2, W$ and $t$. The results are compared to predictions of QCD-inspired models of exclusive vector-meson production.Comment: 24 pages, 8 figure

Measurement of neutral current e+/-p cross sections at high Bjorken x with the ZEUS detector

The neutral current e+/-p cross section has been measured up to values of Bjorken x of approximately 1 with the ZEUS detector at HERA using an integrated luminosity of 187 inv. pb of e-p and 142 inv. pb of e+p collisions at sqrt(s) = 318GeV. Differential cross sections in x and Q2, the exchanged boson virtuality, are presented for Q2 geq 725GeV2. An improved reconstruction method and greatly increased amount of data allows a finer binning in the high-x region of the neutral current cross section and leads to a measurement with much improved precision compared to a similar earlier analysis. The measurements are compared to Standard Model expectations based on a variety of recent parton distribution functions.Comment: 39 pages, 9 figure

Rare region effects at classical, quantum, and non-equilibrium phase transitions

Rare regions, i.e., rare large spatial disorder fluctuations, can dramatically change the properties of a phase transition in a quenched disordered system. In generic classical equilibrium systems, they lead to an essential singularity, the so-called Griffiths singularity, of the free energy in the vicinity of the phase transition. Stronger effects can be observed at zero-temperature quantum phase transitions, at nonequilibrium phase transitions, and in systems with correlated disorder. In some cases, rare regions can actually completely destroy the sharp phase transition by smearing. This topical review presents a unifying framework for rare region effects at weakly disordered classical, quantum, and nonequilibrium phase transitions based on the effective dimensionality of the rare regions. Explicit examples include disordered classical Ising and Heisenberg models, insulating and metallic random quantum magnets, and the disordered contact process.Comment: Topical review, 68 pages, 14 figures, final version as publishe

Search for a narrow charmed baryonic state decaying to D^*+/- p^-/+ in ep collisions at HERA

A resonance search has been made in the D^*+/- p^-/+ invariant-mass spectrum with the ZEUS detector at HERA using an integrated luminosity of 126 pb^-1. The decay channels D^*+ -> D^0 pi^+_s -> (K^- pi^+) pi^+_s and D^*+ -> D^0 pi^+_s -> (K^- pi^+ pi^+ pi^-) pi^+_s (and the corresponding antiparticle decays) were used to identify D^*+/- mesons. No resonance structure was observed in the D^*+/- p^-/+ mass spectrum from more than 60000 reconstructed D^*+/- mesons. The results are not compatible with a report of the H1 Collaboration of a charmed pentaquark, Theta^0_c.Comment: 22 pages, 7 figures, 1 table; minor text revisions; 2 references adde

Search for lepton-flavor violation at HERA

A search for lepton-flavor-violating interactions $e p \to \mu X$ and $e p\to \tau X$ has been performed with the ZEUS detector using the entire HERA I data sample, corresponding to an integrated luminosity of 130 pb^{-1}. The data were taken at center-of-mass energies, $\sqrt{s}$, of 300 and 318 GeV. No evidence of lepton-flavor violation was found, and constraints were derived on leptoquarks (LQs) that could mediate such interactions. For LQ masses below $\sqrt{s}$, limits were set on $\lambda_{eq_1} \sqrt{\beta_{\ell q}}$, where $\lambda_{eq_1}$ is the coupling of the LQ to an electron and a first-generation quark $q_1$, and $\beta_{\ell q}$ is the branching ratio of the LQ to the final-state lepton $\ell$ ($\mu$ or $\tau$) and a quark $q$. For LQ masses much larger than $\sqrt{s}$, limits were set on the four-fermion interaction term $\lambda_{e q_\alpha} \lambda_{\ell q_\beta} / M_{\mathrm{LQ}}^2$ for LQs that couple to an electron and a quark $q_\alpha$ and to a lepton $\ell$ and a quark $q_\beta$, where $\alpha$ and $\beta$ are quark generation indices. Some of the limits are also applicable to lepton-flavor-violating processes mediated by squarks in $R$-Parity-violating supersymmetric models. In some cases, especially when a higher-generation quark is involved and for the process $e p\to \tau X$, the ZEUS limits are the most stringent to date.Comment: 37 pages, 10 figures, Accepted by EPJC. References and 1 figure (Fig. 6) adde