The Goldbeter-Koshland switch in the first-order region and its response to dynamic disorder

In their classical work (Proc. Natl. Acad. Sci. USA, 1981, 78:6840-6844), Goldbeter and Koshland mathematically analyzed a reversible covalent modification system which is highly sensitive to the concentration of effectors. Its signal-response curve appears sigmoidal, constituting a biochemical switch. However, the switch behavior only emerges in the "zero-order region", i.e. when the signal molecule concentration is much lower than that of the substrate it modifies. In this work we showed that the switching behavior can also occur under comparable concentrations of signals and substrates, provided that the signal molecules catalyze the modification reaction in cooperation. We also studied the effect of dynamic disorders on the proposed biochemical switch, in which the enzymatic reaction rates, instead of constant, appear as stochastic functions of time. We showed that the system is robust to dynamic disorder at bulk concentration. But if the dynamic disorder is quasi-static, large fluctuations of the switch response behavior may be observed at low concentrations. Such fluctuation is relevant to many biological functions. It can be reduced by either increasing the conformation interconversion rate of the protein, or correlating the enzymatic reaction rates in the network.Comment: 23 pages, 4 figures, accepted by PLOS ON

Real-time PCR complements immunohistochemistry in the determination of HER-2/neu status in breast cancer

BACKGROUND: The clinical benefit of determining the status of HER-2/neu amplification in breast cancer patients is well accepted. Although immunohistochemistry (IHC) is the most frequently used method to assess the over-expression of HER-2 protein, fluorescent in-situ hybridization (FISH) is recognized as the "gold standard" for the determining of HER-2/neu status. The greatest discordance between the two methods occurs among breast tumors that receive an indeterminate IHC score of 2+. More recently, a real-time polymerase chain reaction (PCR) assay using the LightCycler(® )has been developed for quantifying HER-2/neu gene amplification. In this study, we evaluated the sensitivity and specificity of a commercially available LightCycler assay as it compares to FISH. To determine whether this assay provides an accurate alternative for the determination of HER-2/neu status, we focused primarily on tumors that were deemed indeterminate or borderline status by IHC. METHODS: Thirty-nine breast tumors receiving an IHC score of 2+ were evaluated by both FISH and LightCycler(® )technologies in order to determine whether quantitative real-time PCR provides an accurate alternative for the determination of HER-2/neu status. RESULTS: We found a high concordance (92%) between FISH and real-time PCR results. We also observed that 10% of these tumors were positive for gene amplification by both FISH and real-time PCR. CONCLUSION: The data show that the results obtained for the gene amplification of HER-2/neu by real-time PCR on the LightCycler(® )instrument is comparable to results obtained by FISH. These results therefore suggest that real-time PCR analysis, using the LightCycler(®), is a viable alternative to FISH for reassessing breast tumors which receive an IHC score of 2+, and that a combined IHC and real-time PCR approach for the determination of HER-2 status in breast cancer patients may be an effective and efficient strategy

On Simulated Annealing and the Construction of Linear Spline Approximations for Scattered Data

Abstract. We describe a method to create optimal linear spline approximations to arbitrary functions of one or two variables, given as scattered data without known connectivity. We start with an initial approximation consisting of a fixed number of vertices and improve this approximation by choosing different ver-tices, governed by a simulated annealing algorithm. In the case of one variable, the approximation is defined by line segments; in the case of two variables, the vertices are connected to define a Delaunay triangulation of the selected subset of sites in the plane. In a second version of this algorithm, specifically designed for the bivariate case, we choose vertex sets and also change the triangulation to achieve both optimal vertex placement and optimal triangulation. We then cre-ate a hierarchy of linear spline approximations, each one being a superset of all lower-resolution ones.

A unifying probabilistic framework for analyzing residual dipolar couplings

Residual dipolar couplings provide complementary information to the nuclear Overhauser effect measurements that are traditionally used in biomolecular structure determination by NMR. In a de novo structure determination, however, lack of knowledge about the degree and orientation of molecular alignment complicates the analysis of dipolar coupling data. We present a probabilistic framework for analyzing residual dipolar couplings and demonstrate that it is possible to estimate the atomic coordinates, the complete molecular alignment tensor, and the error of the couplings simultaneously. As a by-product, we also obtain estimates of the uncertainty in the coordinates and the alignment tensor. We show that our approach encompasses existing methods for determining the alignment tensor as special cases, including least squares estimation, histogram fitting, and elimination of an explicit alignment tensor in the restraint energy

Identifying a High Fraction of the Human Genome to be under Selective Constraint Using GERP++

Computational efforts to identify functional elements within genomes leverage comparative sequence information by looking for regions that exhibit evidence of selective constraint. One way of detecting constrained elements is to follow a bottom-up approach by computing constraint scores for individual positions of a multiple alignment and then defining constrained elements as segments of contiguous, highly scoring nucleotide positions. Here we present GERP++, a new tool that uses maximum likelihood evolutionary rate estimation for position-specific scoring and, in contrast to previous bottom-up methods, a novel dynamic programming approach to subsequently define constrained elements. GERP++ evaluates a richer set of candidate element breakpoints and ranks them based on statistical significance, eliminating the need for biased heuristic extension techniques. Using GERP++ we identify over 1.3 million constrained elements spanning over 7% of the human genome. We predict a higher fraction than earlier estimates largely due to the annotation of longer constrained elements, which improves one to one correspondence between predicted elements with known functional sequences. GERP++ is an efficient and effective tool to provide both nucleotide- and element-level constraint scores within deep multiple sequence alignments

Antisense DNA parameters derived from next-nearest-neighbor analysis of experimental data

<p>Abstract</p> <p>Background</p> <p>The enumeration of tetrameric and other sequence motifs that are positively or negatively correlated with <it>in vivo </it>antisense DNA effects has been a useful addition to the arsenal of information needed to predict effective targets for antisense DNA control of gene expression. Such retrospective information derived from <it>in vivo </it>cellular experiments characterizes aspects of the sequence dependence of antisense inhibition that are not predicted by nearest-neighbor (NN) thermodynamic parameters derived from <it>in vitro </it>experiments. However, quantitation of the antisense contributions of motifs is problematic, since individual motifs are not isolated from the effects of neighboring nucleotides, and motifs may be overlapping. These problems are circumvented by a next-nearest-neighbor (NNN) analysis of antisense DNA effects in which the overlapping nature of nearest-neighbors is taken into account.</p> <p>Results</p> <p>Next-nearest-neighbor triplet combinations of nucleotides are the simplest that include overlapping sequence effects and therefore can encompass interactions beyond those of nearest neighbors. We used singular value decomposition (SVD) to fit experimental data from our laboratory in which phosphorothioate-modified antisense DNAs (S-DNAs) 20 nucleotides long were used to inhibit cellular protein expression in 112 experiments involving four gene targets and two cell lines. Data were fitted using a NNN model, neglecting end effects, to derive NNN inhibition parameters that could be combined to give parameters for a set of 49 sequences that represents the inhibitory effects of all possible overlapping triplet interactions in the cellular targets of these antisense S-DNAs. We also show that parameters to describe subsets of the data, such as the mRNAs being targeted and the cell lines used, can be included in such a derivation. While NNN triplet parameters provided an adequate model to fit our data, NN doublet parameters did not.</p> <p>Conclusions</p> <p>The methodology presented illustrates how NNN antisense inhibitory information can be derived from <it>in vivo </it>cellular experiments. Subsequent calculations of the antisense inhibitory parameters for any mRNA target sequence automatically take into account the effects of all possible overlapping combinations of nearest-neighbors in the sequence. This procedure is more robust than the tallying of tetrameric motifs that have positive or negative antisense effects. The specific parameters derived in this work are limited in their applicability by the relatively small database of experiments that was used in their derivation.</p

131I-metaiodobenzylguanidine (131I-MIBG) therapy for residual neuroblastoma: a mono-institutional experience with 43 patients

Incomplete response to therapy may compromise the outcome of children with advanced neuroblastoma. In an attempt to improve tumour response we incorporated 131I-metaiodobenzylguanidine (131I-MIBG) in the treatment regimens of selected stage 3 and stage 4 patients. Between 1986 and 1997, 43 neuroblastoma patients older than 1 year at diagnosis, 13 with stage 3 (group A) and 30 with stage 4 disease (group B) who had completed the first-line protocol without achieving complete response entered in this study. 131I-MIBG dose/course ranged from 2.5 to 5.5 Gbq (median, 3.7). The number of courses ranged from 1 to 5 (median 3) depending on the tumour response and toxicity. The most common acute side-effect was thrombocytopenia. Later side-effects included severe interstitial pneumonia in one patient, acute myeloid leukaemia in two, reduced thyroid reserve in 21. Complete response was documented in one stage 4 patient, partial response in 12 (two stage 3, 10 stage 4), mixed or no response in 25 (ten stage 3, 15 stage 4) and disease progression in five (one stage 3, four stage 4) Twenty-four patients (12/13 stage 3, 12/30 stage 4) are alive at 22–153 months (median, 59) from diagnosis. 131I-MIBG therapy may increase the cure rate of stage 3 and improve the response of stage 4 neuroblastoma patients with residual disease after first-line therapy. A larger number of patients should be treated to confirm these results but logistic problems hamper prospective and coordinated studies. Long-term toxicity can be severe. © 1999 Cancer Research Campaig

The Formation and Evolution of the First Massive Black Holes

The first massive astrophysical black holes likely formed at high redshifts (z>10) at the centers of low mass (~10^6 Msun) dark matter concentrations. These black holes grow by mergers and gas accretion, evolve into the population of bright quasars observed at lower redshifts, and eventually leave the supermassive black hole remnants that are ubiquitous at the centers of galaxies in the nearby universe. The astrophysical processes responsible for the formation of the earliest seed black holes are poorly understood. The purpose of this review is threefold: (1) to describe theoretical expectations for the formation and growth of the earliest black holes within the general paradigm of hierarchical cold dark matter cosmologies, (2) to summarize several relevant recent observations that have implications for the formation of the earliest black holes, and (3) to look into the future and assess the power of forthcoming observations to probe the physics of the first active galactic nuclei.Comment: 39 pages, review for "Supermassive Black Holes in the Distant Universe", Ed. A. J. Barger, Kluwer Academic Publisher

The Pioneer Anomaly

Radio-metric Doppler tracking data received from the Pioneer 10 and 11 spacecraft from heliocentric distances of 20-70 AU has consistently indicated the presence of a small, anomalous, blue-shifted frequency drift uniformly changing with a rate of ~6 x 10^{-9} Hz/s. Ultimately, the drift was interpreted as a constant sunward deceleration of each particular spacecraft at the level of a_P = (8.74 +/- 1.33) x 10^{-10} m/s^2. This apparent violation of the Newton's gravitational inverse-square law has become known as the Pioneer anomaly; the nature of this anomaly remains unexplained. In this review, we summarize the current knowledge of the physical properties of the anomaly and the conditions that led to its detection and characterization. We review various mechanisms proposed to explain the anomaly and discuss the current state of efforts to determine its nature. A comprehensive new investigation of the anomalous behavior of the two Pioneers has begun recently. The new efforts rely on the much-extended set of radio-metric Doppler data for both spacecraft in conjunction with the newly available complete record of their telemetry files and a large archive of original project documentation. As the new study is yet to report its findings, this review provides the necessary background for the new results to appear in the near future. In particular, we provide a significant amount of information on the design, operations and behavior of the two Pioneers during their entire missions, including descriptions of various data formats and techniques used for their navigation and radio-science data analysis. As most of this information was recovered relatively recently, it was not used in the previous studies of the Pioneer anomaly, but it is critical for the new investigation.Comment: 165 pages, 40 figures, 16 tables; accepted for publication in Living Reviews in Relativit

Statistical Analyses Support Power Law Distributions Found in Neuronal Avalanches

The size distribution of neuronal avalanches in cortical networks has been reported to follow a power law distribution with exponent close to −1.5, which is a reflection of long-range spatial correlations in spontaneous neuronal activity. However, identifying power law scaling in empirical data can be difficult and sometimes controversial. In the present study, we tested the power law hypothesis for neuronal avalanches by using more stringent statistical analyses. In particular, we performed the following steps: (i) analysis of finite-size scaling to identify scale-free dynamics in neuronal avalanches, (ii) model parameter estimation to determine the specific exponent of the power law, and (iii) comparison of the power law to alternative model distributions. Consistent with critical state dynamics, avalanche size distributions exhibited robust scaling behavior in which the maximum avalanche size was limited only by the spatial extent of sampling (“finite size” effect). This scale-free dynamics suggests the power law as a model for the distribution of avalanche sizes. Using both the Kolmogorov-Smirnov statistic and a maximum likelihood approach, we found the slope to be close to −1.5, which is in line with previous reports. Finally, the power law model for neuronal avalanches was compared to the exponential and to various heavy-tail distributions based on the Kolmogorov-Smirnov distance and by using a log-likelihood ratio test. Both the power law distribution without and with exponential cut-off provided significantly better fits to the cluster size distributions in neuronal avalanches than the exponential, the lognormal and the gamma distribution. In summary, our findings strongly support the power law scaling in neuronal avalanches, providing further evidence for critical state dynamics in superficial layers of cortex