Statistical modeling of large microarray data sets to identify stimulus-response profiles

A statistical modeling approach is proposed for use in searching large microarray data sets for genes that have a transcriptional response to a stimulus. The approach is unrestricted with respect to the timing, magnitude or duration of the response, or the overall abundance of the transcript. The statistical model makes an accommodation for systematic heterogeneity in expression levels. Corresponding data analyses provide gene-specific information, and the approach provides a means for evaluating the statistical significance of such information. To illustrate this strategy we have derived a model to depict the profile expected for a periodically transcribed gene and used it to look for budding yeast transcripts that adhere to this profile. Using objective criteria, this method identifies 81% of the known periodic transcripts and 1,088 genes, which show significant periodicity in at least one of the three data sets analyzed. However, only one-quarter of these genes show significant oscillations in at least two data sets and can be classified as periodic with high confidence. The method provides estimates of the mean activation and deactivation times, induced and basal expression levels, and statistical measures of the precision of these estimates for each periodic transcript.

Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets

An integrated understanding of molecular and developmental biology must consider the large number of molecular species involved and the low concentrations of many species in vivo. Quantitative stochastic models of molecular interaction networks can be expressed as stochastic Petri nets (SPNs), a mathematical formalism developed in computer science. Existing software can be used to define molecular interaction networks as SPNs and solve such models for the probability distributions of molecular species. This approach allows biologists to focus on the content of models and their interpretation, rather than their implementation. The standardized format of SPNs also facilitates the replication, extension, and transfer of models between researchers. A simple chemical system is presented to demonstrate the link between stochastic models of molecular interactions and SPNs. The approach is illustrated with examples of models of genetic and biochemical phenomena where the UltraSAN package is used to present results from numerical analysis and the outcome of simulations.

A numerical strategy for efficient modeling of the equatorial wave guide

Convection in the tropics is observed to involve a wide-ranging hierarchy of scales from a few kilometers to the planetary scales and also has a profound impact on short-term climate. The mechanisms responsible for this behavior present a major unsolved problem. A promising emerging approach to address these issues is cloud-resolving modeling. Here a family of numerical models is introduced specifically to model the feedback of small-scale deep convection on tropical planetary waves and tropical circulation in a highly efficient manner compatible with the approach through cloud-resolving modeling. Such a procedure is also useful for theoretical purposes. The basic idea in the approach is to use low-order truncation in the meriodonal direction through Gauss–Hermite quadrature projected onto a simple discrete radiation condition. In this fashion, the cloud-resolving modeling of equatorially trapped planetary waves reduces to the solution of a small number of purely zonal two-dimensional wave systems along a few judiciously chosen meriodonal layers that are coupled only by some additional source terms. The approach is analyzed in detail with full mathematical rigor for linearized equatorial primitive equations with source terms.

A combined kinetic and modeling study of the catalytic center subsites of human angiogenin.

Kinetic analysis and molecular modeling have been used to map the ribonucleolytic center of angiogenin (Ang). Pyrimidine nucleotides were found to interact very weakly with Ang, consistent with the inaccessible B1 pyrimidine binding site revealed by x-ray crystallography. Ang also lacks an effective phosphate binding site on the 5' side of B1. Although the B2 site that preferentially binds purines on the 3' side of B1 is also weak, its associated phosphate subsites make substantial contributions: both 3',5'-ADP and 5'-ADP have Ki values 6-fold lower than for 5'-AMP, and adding a 3'-phosphate to the substrate CpA increases Kcat/Km by 9-fold. Thus Ang has a functional P2 site on the 3' side of B2 and a site for a second phosphate on the 5' side of B2. Modeling of an Ang-d(ApTpApA) complex suggested that Arg-5 forms part of the P2 site and that a 2'-phosphate might bind more tightly than a 3'-phosphate. Both predictions were confirmed kinetically. The subsite map obtained by this combined approach indicated that 5'-diphosphoadenosine 2'-phosphate might be a more potent inhibitor than any of the nucleotides tested thus far. Indeed, its Ki value of 150 microM is 50-fold lower than that for the best nucleotide previously reported and 400-fold lower than the Km for the best dinucleotide substrate. This compound may serve as a suitable starting point for the eventual design of tight-binding inhibitors of Ang as antiangiogenic agents for human therapy.