Search

This work extends the work of Whitlock in examining the critical effective population sizes from the fixation of both deleterious and beneficial mutations under drift and selection to prevent mutation breakdown of the population. The validity of approximations for the probability of fixation depends on the nature of the assumed distribution for the fitness effect of both types of mutations. Using no approximation for the probability of fixation and assuming a heavy tailed fitness effect distribution, the current model indicates that the coefficients of variation for the fitness effect distributions of both types of mutations and the fitness effect distribution mean for the beneficial mutations are important predictors of the critical effective population size. The current model further predicts that very small populations can be sustained if the fitness effect variances for both types of mutations and the mean for beneficial mutations are large.

Earlier work rigorously derived a general probabilistic model for the PCR process that includes as a special case the Velikanov-Kapral model where all nucleotide reaction rates are the same. In this model, the probability of binding of deoxy-nucleoside triphosphate (dNTP) molecules with template strands is derived from the microscopic chemical kinetics. A recursive solution for the probability function of binding of dNTPs is developed for a single cycle and is used to calculate expected yield for a multicycle PCR. The model is able to reproduce important features of the PCR amplification process quantitatively. With a set of favorable reaction conditions, the amplification of the target sequence is fast enough to rapidly outnumber all side products. Furthermore, the final yield of the target sequence in a multicycle PCR run always approaches an asymptotic limit that is less than one. The amplification process itself is highly sensitive to initial concentrations and the reaction rates of addition to the template strand of each type of dNTP in the solution. This paper extends the earlier Saha model with a physics based model of the dependence of the reaction rates on temperature, and estimates parameters in this new model by nonlinear regression. The calibrated model is validated using RT-PCR data.

Parameter values for a kinetic model of the nuclear replication-division cycle in frog eggs are estimated by fitting solutions of the kinetic equations (nonlinear ordinary differential equations) to a suite of experimental observations. A set of optimal parameter values is found by minimizing an objective function defined as the orthogonal distance between the data and the model. The differential equations are solved by LSODAR and the objective function is minimized by ODRPACK. The optimal parameter values are close to the "guesstimates" of the modelers who first studied this problem. These tools are sufficiently general to attack more complicated problems, where guesstimation is impractical or unreliable.

The life of a cell is governed by the physicochemical properties of a complex network of interacting macromolecules (primarily genes and proteins). Hence, a full scientific understanding of and rational engineering approach to cell physiology require accurate mathematical models of the spatial and temporal dynamics of these macromolecular assemblies, especially the networks involved in integrating signals and regulating cellular responses. The Virginia Tech Consortium is involved in three specific goals of DARPA's computational biology program (Bio-COMP): to create effective software tools for modeling gene-protein-metabolite networks, to employ these tools in creating a new generation of realistic models, and to test and refine these models by well-conceived experimental studies. The special emphasis of this group is to understand the mechanisms of cell cycle control in eukaryotes (yeast cells and frog eggs). The software tools developed at Virginia Tech are designed to meet general requirements of modeling regulatory networks and are collected in a problem-solving environment called JigCell.