PubMed Health⌕ Search

Biomedical subjects

Jeffrey J Morgan

Publications and source records attributed to Jeffrey J Morgan.

4 recordsLinked to original sources

Whole-cell modeling framework in which biochemical dynamics impact aspects of cellular geometry.

A mathematical framework for modeling biological cells from a physicochemical perspective is described. Cells modeled within this framework consist of at least two regions, including a cytosolic volume encapsulated by a membrane surface. The cytosol is viewed as a well-stirred chemical reactor capable of changing volume while the membrane is assumed to be an oriented 2-D surface capable of changing surface area. Two physical properties of the cell, namely volume and surface area, are determined by (and determine) the reaction dynamics generated from a set of chemical reactions designed to be occurring in the cell. This framework allows the modeling of complex cellular behaviors, including self-replication. This capability is illustrated by constructing two self-replicating prototypical whole-cell models. One protocell was designed to be of minimal complexity; the other to incorporate a previously reported well-known mechanism of the eukaryotic cell cycle. In both cases, self-replicative behavior was achieved by seeking stable physically possible oscillations in concentrations and surface-to-volume ratio, and by synchronizing the period of such oscillations to the doubling of cytosolic volume and membrane surface area. Rather than being enforced externally or artificially, growth and division occur naturally as a consequence of the assumed chemical mechanism operating within the framework.

Animals↗

A framework for whole-cell mathematical modeling.

The default framework for modeling biochemical processes is that of a constant-volume reactor operating under steady-state conditions. This is satisfactory for many applications, but not for modeling growth and division of cells. In this study, a whole-cell modeling framework is developed that assumes expanding volumes and a cell-division cycle. A spherical newborn cell is designed to grow in volume during the growth phase of the cycle. After 80% of the cycle period, the cell begins to divide by constricting about its equator, ultimately affording two spherical cells with total volume equal to twice that of the original. The cell is partitioned into two regions or volumes, namely the cytoplasm (Vcyt) and membrane (Vmem), with molecular components present in each. Both volumes change during the cell cycle; Vcyt changes in response to osmotic pressure changes as nutrients enter the cell from the environment, while Vmem changes in response to this osmotic pressure effect such that membrane thickness remains invariant. The two volumes change at different rates; in most cases, this imposes periodic or oscillatory behavior on all components within the cell. Since the framework itself rather than a particular set of reactions and components is responsible for this behavior, it should be possible to model various biochemical processes within it, affording stable periodic solutions without requiring that the biochemical process itself generates oscillations as an inherent feature. Given that these processes naturally occur in growing and dividing cells, it is reasonable to conclude that the dynamics of component concentrations will be more realistic than when modeled within constant-volume and/or steady-state frameworks. This approach is illustrated using a symbolic whole cell model.

Animals↗

Dynamic responses of protein homeostatic regulatory mechanisms to perturbations from steady state.

Nineteen hypothetical protein homeostatic regulatory mechanisms were constructed and analysed in terms of the rate at which they recovered from a perturbation in the steady-state concentration of any component. Systems were constructed to symbolize transcription/translation processes of the average protein from Escherichia coli (1000 copies of protein P along with 1 gene G per cell). In some model systems, G catalysed the synthesis of P directly, while in others G catalysed the synthesis of mRNA (called M), and M catalysed the synthesis of P in a subsequent step. Recovery rates for each regulatory mechanism were obtained by generating the corresponding system of differential equations, linearizing the system about the steady state, and determining eigenvalues of the associated coefficient matrix. The optimal rate of recovery for a given mechanism, R(D), was determined by combining random and gradient search approaches to find rate constants for which the system recovered fastest. Regulatory elements that improved dynamic regulation were identified. These consisted of negative feedback relationships that involved P binding to either G (to shut off the synthesis of P) or M (to stimulate its degradation). Regulation improved as increasing numbers of P's bound to either G or M; however, the binding to M was more effective. In other mechanisms PP dimers bound G. Dimer-binding mechanisms were roughly twice as effective in terms of regulation as those that bound P monomers. The effect of linking two regulatory "modules" was also investigated. Linking had no effect on R(D), but optimal rate constants for the linked system were similar to those of the unlinked modules, suggesting that it may be feasible to construct regulatory networks by linking individual modules of this type.

Animals↗

Analysis of protein homeostatic regulatory mechanisms in perturbed environments at steady state.

Nine different protein homeostatic regulatory mechanisms were analysed for their ability to maintain a generic protein P within a specified range of a set-point steady-state concentration while perturbed by external processes that altered the rates at which P was produced and/or consumed. Steady state regulatory effectiveness was defined by the area within a rectangular region of "perturbation space", where axes correspond to rates of positive and negative perturbations. The size of this region differed in accordance with the regulatory elements composing the homeostatic mechanism. Such elements included basic negative feedback control of transcription (in which P, at some high concentration relative to its set-point value, binds to the gene G that encodes it, thereby inhibiting transcription), multiple sequential binding of a feedback effector (two P's bind sequentially to G), and dimerization of a feedback effector (a P(2) dimer binds to G). Two homeostatic mechanisms included a cascade structure, one with and one without translational feedback control. Another mechanism included feedback control of P degradation. Finally, two mechanisms illustrated the limits of regulatory systems. One lacked all regulatory elements (and included only an invariant rate of P synthesis and degradation) while the other assumed perfect (Boolean) regulation, in which transcription is completely inhibited at [P]>[P](sp) and is fully active at [P]<[P](sp). All of the systems evaluated are known, but the analytical expressions developed here allow quantitative comparisons between them. These expressions were evaluated at values typical of the average protein in Escherichia coli. A method for building regulatory networks by linking semi-independent regulatory modules is discussed.

Animals↗