Transient responses of budding yeast populations

The unequal-division model for budding yeast is used to formulate a population-balance model for the transient behavior of populations of these organisms. The model consists of linear partial differential equations coupled through algebraic equations. It is shown how the solution of this system of equations can be obtained in a systematic stepwise fashion. The special case of a population subjected to a step change in growth rate is described in some detail, and solutions for two special cases are determined for transients following an age-distribution perturbation. It is shown how experimental data on transient behavior of a cell population can yield information on single-cell mass-synthesis kinetics and on the manner in which individual cells control certain critical parameters in the cell cycle.     

Sedimentation Patterns of Rapidly Reversible Protein Interactions

The transport behavior of macromolecular mixtures with rapidly reversible complex formation is of great interest in the study of protein interactions by many different methods. Complicated transport patterns arise even for simple bimolecular reactions, when all species exhibit different migration velocities. Although partial differential equations are available to describe the spatial and temporal evolution of the interacting system given particular initial conditions, a general overview of the phase behavior of the systems in parameter space has not yet been reported. In the case of sedimentation of two-component mixtures, this study presents simple analytical solutions that solve the underlying equations in the diffusion-free limit previously subject to Gilbert-Jenkins theory. The new expressions describe, with high precision, the average sedimentation coefficients and composition of each boundary, which allow the examination of features of the whole parameter space at once, and may be used for experimental design and robust analysis of experimental boundary patterns to derive the stoichiometry and affinity of the complex. This study finds previously unrecognized features, including a phase transition between boundary patterns. The model reveals that the time-average velocities of all components in the reaction mixture must match—a condition that suggests an intuitive physical picture of an effective particle of the coupled cosedimentation of an interacting system. Adding to the existing numerical solutions of the relevant partial differential equations, the effective particle model provides physical insights into the relationships of the parameters that govern sedimentation patterns.     

Capturing Intracellular pH Dynamics by Coupling Its Molecular Mechanisms within a Fully Tractable Mathematical Model

We describe the construction of a fully tractable mathematical model for intracellular pH. This work is based on coupling the kinetic equations depicting the molecular mechanisms for pumps, transporters and chemical reactions, which determine this parameter in eukaryotic cells. Thus, our system also calculates the membrane potential and the cytosolic ionic composition. Such a model required the development of a novel algebraic method that couples differential equations for slow relaxation processes to steady-state equations for fast chemical reactions. Compared to classical heuristic approaches based on fitted curves and ad hoc constants, this yields significant improvements. This model is mathematically self-consistent and allows for the first time to establish analytical solutions for steady-state pH and a reduced differential equation for pH regulation. Because of its modular structure, it can integrate any additional mechanism that will directly or indirectly affect pH. In addition, it provides mathematical clarifications for widely observed biological phenomena such as overshooting in regulatory loops. Finally, instead of including a limited set of experimental results to fit our model, we show examples of numerical calculations that are extremely consistent with the wide body of intracellular pH experimental measurements gathered by different groups in many different cellular systems.     

Stability analysis for a peri-implant osseointegration model

We investigate stability of the solution of a set of partial differential equations, which is used to model a peri-implant osseointegration process. For certain parameter values, the solution has a ‘wave-like’ profile, which appears in the distribution of osteogenic cells, osteoblasts, growth factor and bone matrix. This ‘wave-like’ profile contradicts experimental observations. In our study we investigate the conditions, under which such profile appears in the solution. Those conditions are determined in terms of model parameters, by means of linear stability analysis, carried out at one of the constant solutions of the simplified system. The stability analysis was carried out for the reduced system of PDE’s, of which we prove, that it is equivalent to the original system of equations, with respect to the stability properties of constant solutions. The conclusions, derived from the linear stability analysis, are extended for the case of large perturbations. If the constant solution is unstable, then the solution of the system never converges to this constant solution. The analytical results are validated with finite element simulations. The simulations show, that stability of the constant solution could determine the behavior of the solution of the whole system, if certain initial conditions are considered.     

The geometrical analysis of a predator-prey model with two state impulses

Using successor functions and Poincaré-Bendixson theorem of impulsive differential equations, the existence of periodical solutions to a predator-prey model with two state impulses is investigated. By stability theorem of periodic solution to impulsive differential equations, the stability conditions of periodic solutions to the system are given. Some simulations are exerted to prove the results.     

An analytical model of the hydraulic aspects of stomatal dynamics

An analytical model of the hydraulic aspects of stomatal dynamics is formulated in this paper. The model consists of a coupled system of non-linear, ordinary differential equations, written in terms of water potentials, hydrostatic pressures, osmotic potentials, water vapor resistances and water fluxes. The model is validated by comparisons with the experimental literature. Numerical solutions of the model show qualitative agreement with most known stomatal responses.Stomatal opening in the model is dependent on the interaction of the guard and subsidiary cells in the following manner. Pore opening is initiated by a rise in the guard cell hydrostatic pressure. As the stomate opens, transpiration increases, causing the cell wall water potential to drop. The drop in cell wall water potential then causes the subsidiary cell pressure to drop, opening is accelerated, and the stomate literally “pops” open. Simulated opening proceeds in two distinct phases: a stress phase and a motor phase. During the stress phase, guard cell pressure rises but the pore remains closed. The motor phase commences when the guard cell pressure has risen sufficiently to initiate pore opening, beyond which point opening progresses rapidly.Hydropassive stomatal movements are found to be insufficient to regulate water loss at low leaf water potentials. Stable, hydraulically-based oscillations in stomatal aperture are shown in the model by the existence of a stable limit cycle. The period of these oscillations is strongly influenced by the cell membrane hydraulic conductivity. An increased conductivity results in a shorter period oscillation. Environmental conditions promoting oscillatory behavior are in qualitative agreement with the experimental literature.     

Electrical properties of frog skeletal muscle fibers interpreted with a mesh model of the tubular system.

This paper presents the construction, derivation, and test of a mesh model for the electrical properties of the transverse tubular system (T-system) in skeletal muscle. We model the irregular system of tubules as a random network of miniature transmission lines, using differential equations to describe the potential between the nodes and difference equations to describe the potential at the nodes. The solution to the equations can be accurately represented in several approximate forms with simple physical and graphical interpretations. All the parameters of the solution are specified by impedance and morphometric measurements. The effect of wide circumferential spacing between T-system openings is analyzed and the resulting restricted mesh model is shown to be approximated by a mesh with an access resistance. The continuous limit of the mesh model is shown to have the same form as the disk model of the T-system, but with a different expression for the tortuosity factor. The physical meaning of the tortuosity factor is examined, and a short derivation of the disk model is presented that gives results identical to the continuous limit of the mesh model. Both the mesh and restricted mesh models are compared with experimental data on the impedance of muscle fibers of the frog sartorius. The derived value for the resistivity of the lumen of the tubules is not too different from that of the bathing solution, the difference probably arising from the sensitivity of this value to errors in the morphometric measurements.     

A mathematical model for fibroblast and collagen orientation

Due to the increasing importance of the extracellular matrix in many biological problems, in this paper we develop a model for fibroblast and collagen orientation with the ultimate objective of understanding how fibroblasts form and remodel the extracellular matrix, in particular its collagen component. The model uses integrodifferential equations to describe the interaction between the cells and fibers at a point in space with various orientations. The equations are studied both analytically and numerically to discover different types of solutions and their behavior. In particular we examine solutions where all the fibroblasts and collagen have discrete orientations, a localized continuum of orientations and a continuous distribution of orientations with several maxima. The effect of altering the parameters in the system is explored, including the angular diffusion coefficient for the fibroblasts, as well as the strength and range of the interaction between fibroblasts and collagen. We find the initial conditions and the range of influence between the collagen and the fibroblasts are the two factors which determine the behavior of the solutions. The implications of this for wound healing and cancer are discussed including the conclusion that the major factor in determining the degree of scarring is the initial deposition of collagen.     

The effect of the yield expression on the existence of oscillatory behavior in a three-variable model of a continuous fermentation system subject to product inhibition

A three-variable model of a continuous fermentation process characterised by product inhibition is studied. It is shown that if the cell to substrate yield is constant, the system cannot have periodic solutions. If, on the other hand, the yield term is a variable function of substrate concentration, the model will exhibit oscillations in the cells, substrate and product concentrations in the form of Hopf bifurcation in the underlying system of three nonlinear, ordinary differential equations which comprise the model.     

Current status and challenges in connecting models of erythrocyte metabolism to experimental reality

Detailed kinetic models of human erythrocyte metabolism have served to summarize the vast literature and to predict outcomes from laboratory and “Nature's” experiments on this simple cell. Mathematical methods for handling the large array of nonlinear ordinary differential equations that describe the time dependence of this system are well developed, but experimental methods that can guide the evolution of the models are in short supply. NMR spectroscopy is one method that is non-selective with respect to analyte detection but is highly specific with respect to their identification and quantification. Thus time courses of metabolism are readily recorded for easily changed experimental conditions. While the data can be simulated, the systems of equations are too complex to allow solutions of the inverse problem, namely parameter-value estimation for the large number of enzyme and membrane-transport reactions operating in situ as opposed to in vitro. Other complications with the modelling include the dependence of cell volume on time, and the rates of membrane transport processes are often dependent on the membrane potential. These matters are discussed in the light of new modelling strategies.     

A new method of identifying a systems based on the minimal quadratic discrepancy criterion for biophysics problems

A wide range of biophysical systems are described by nonlinear dynamic models mathematically presented as a set of ordinary differential equations in the Cauchy explicit form: [formula: see text] Fij(X1(t),..,XN(t),t), (i = 1,...,N, j = 1,...,M), where Fij (X1(t), ..., XN(t), t) is a set of basis functions satisfying the Lipschitz condition. We investigate the problem of evaluation of model constants aij (the system identification) using experimental data about the time dependence of the dynamic parameters of the system Xi(t). A new method of system identification for the class of similar nonlinear dynamic models is proposed. It is shown that the problem of identifying an initial nonlinear model can be reduced to the solution of a system of linear equations for the matrix of the dynamic model constants [aj]i. It is proposed to determine the set of dynamic model constants aij using the criterion of minimal quadratic discrepancy for the time dependence of the set of dynamic parameters Xi(t). An important special case of the nonlinear model, the quadratic model, is considered. Test problems of identification using this method are presented for two nonlinear systems: the Van der Pol type multiparametric nonlinear oscillator and the strange attractor of Ressler, a widely known example of dynamic systems showing the stochastic behavior.     

Positive solutions of an elliptic system arising from a model in evolutionary ecology

We give sufficient and almost necessary conditions for the existence of positive solutions to an elliptic system satisfying various Dirichlet boundary conditions. The elliptic system consists of the steady-state equations of a parabolic system used to model the growth and spread of a particular gene and population living in a bounded region. The model takes into account the fact that the fitness of the individuals in the population may depend on the population size. Some non-existence results are also included.Research partially supported by NSF grant no. DMS-8801968     

A model of cell cycle behavior dominated by kinetics of a pathway stimulated by growth factors

A modified version of a previously developed mathematical model [Obeyesekere et al., Cell Prolif. (1997)] of the G1-phase of the cell cycle is presented. This model describes the regulation of the G1-phase that includes the interactions of the nuclear proteins, RB, cyclin E, cyclin D, cdk2, cdk4 and E2F. The effects of the growth factors on cyclin D synthesis under saturated or unsaturated growth factor conditions are investigated based on this model. The solutions to this model (a system of nonlinear ordinary differential equations) are discussed with respect to existing experiments. Predictions based on mathematical analysis of this model are presented. In particular, results are presented on the existence of two stablesolutions, i. e., bistability within the G1-phase. It is shown that this bistability exists under unsaturated growth factor concentration levels. This phenomenon is very noticeable if the efficiency of the signal transduction, initiated by the growth factors leading to cyclin D synthesis, is low. The biological significance of this result as well as possible experimental designs to test these predictions are presented.     

Qualitative Analysis of an ODE Model of a Class of Enzymatic Reactions

The present paper analyzes an ODE model of a certain class of (open) enzymatic reactions. This type of model is used, for instance, to describe the interactions between messenger RNAs and microRNAs. It is shown that solutions defined by positive initial conditions are well defined and bounded on \([0, \infty )\) and that the positive octant of \({\mathbb {R}}^3\) is a positively invariant set. We prove further that in this positive octant there exists a unique equilibrium point, which is asymptotically stable and a global attractor for any initial state with positive components; a controllability property is emphasized. We also investigate the qualitative behavior of the QSSA system in the phase plane \({\mathbb {R}}^2\). For this planar system we obtain similar results regarding global stability by using Lyapunov theory, invariant regions and controllability.