Quantitative Biology

[abridged] Background: The distribution of chemical species in an open system at metastable equilibrium can be expressed as a function of environmental variables which can include temperature, oxidation-reduction potential and others. Calculations of metastable equilibrium for various model systems were used to characterize chemical transformations among proteins and groups of proteins found in different compartments of yeast cells. Results: With increasing oxygen fugacity, the relative metastability fields of model proteins for major subcellular compartments go as mitochondrion, endoplasmic reticulum, cytoplasm, nucleus. In a metastable equilibrium setting at relatively high oxygen fugacity, proteins making up actin are predominant, but those constituting the microtubule occur with a low chemical activity. A reaction sequence involving the microtubule and spindle pole proteins was predicted by combining the known intercompartmental interactions with a hypothetical program of oxygen fugacity changes in the local environment. In further calculations, the most-abundant proteins within compartments generally occur in relative abundances that only weakly correspond to a metastable equilibrium distribution. However, physiological populations of proteins that form complexes often show an overall positive or negative correlation with the relative abundances of proteins in metastable assemblages. Conclusions: This study explored the outlines of a thermodynamic description of chemical transformations among interacting proteins in yeast cells. The results suggest that these methods can be used to measure the degree of departure of a natural biochemical process or population from a local minimum in Gibbs energy.

The plane of bacterial cell division must be precisely positioned. In the bacterium Myxococcus xanthus, the proteins PomX and PomY form a large cluster, which is tethered to the nucleoid by the ATPase PomZ and moves in a stochastic, but biased manner towards midcell, where it initiates cell division. Previously, a positioning mechanism based on the fluxes of PomZ on the nucleoid was proposed. However, the cluster dynamics was analyzed in a reduced, one-dimensional geometry. Here we introduce a mathematical model that accounts for the three-dimensional shape of the nucleoid, such that nucleoid-bound PomZ dimers can diffuse past the cluster without interacting with it. Using stochastic simulations, we find that the cluster still moves to and localizes at midcell. Redistribution of PomZ by diffusion in the cytosol is essential for this cluster dynamics. Our mechanism also positions two clusters equidistantly on the nucleoid. We conclude that a flux-based mechanism allows for cluster positioning in a biologically realistic three-dimensional cell geometry.

This article presents a physical biology approach to understanding organization and segregation of bacterial chromosomes. The author uses a "piston" analogy for bacterial chromosomes in a cell, which leads to a phase diagram for the organization of two athermal chains confined in a closed geometry characterized by two length scales (length and width). When applied to rod-shaped bacteria such as Escherichia coli, this phase diagram predicts that, despite strong confinement, duplicated chromosomes will demix, i.e., there exists a primordial physical driving force for chromosome segregation. The author discusses segregation of duplicating chromosomes using the concentric-shell model, which predicts that newly synthesized DNA will be found in the periphery of the chromosome during replication. In contrast to chromosomes, these results suggest that most plasmids will be randomly distributed inside the cell because of their small sizes. An active partitioning system is therefore required for accurate segregation of low-copy number plasmids. Implications of these results are also sketched, e.g., on the role of proteins, segregation mechanisms for bacteria of diverse shapes, cell cycle of an artificial cell, and evolution.

The prospects of application of graphene and related structures as the membrane mimetic materials, capable of reproducing several biomembrane functions up to the certain limit, are analyzed in the series of our papers. This paper considers the possibility of the ion channel function modeling using graphene and its derivatives. The physical mechanisms providing selective permeability for different membrane mimetic materials, as well as the limits of the adequate simulation of the transport, catalytic, sensing and electrogenic properties of the cell membrane ion channels using bilayered graphene-based structures are discussed.

An endogenous molecular-cellular network for both normal and abnormal functions is assumed to exist. This endogenous network forms a nonlinear stochastic dynamical system, with many stable attractors in its functional landscape. Normal or abnormal robust states can be decided by this network in a manner similar to the neural network. In this context cancer is hypothesized as one of its robust intrinsic states. This hypothesis implies that a nonlinear stochastic mathematical cancer model is constructible based on available experimental data and its quantitative prediction is directly testable. Within such model the genesis and progression of cancer may be viewed as stochastic transitions between different attractors. Thus it further suggests that progressions are not arbitrary. Other important issues on cancer, such as genetic vs epigenetics, double-edge effect, dormancy, are discussed in the light of present hypothesis. A different set of strategies for cancer prevention, cure, and care, is therefore suggested.

Transmembrane receptor clustering is a ubiquitous phenomenon in pro- and eukaryotic cells to physically sense receptor/ligand interactions and subsequently translate an exogenous signal into a cellular response. Despite that receptor cluster formation has been described for a wide variety of receptors, ranging from chemotactic receptors in bacteria to growth factor and neurotransmitter receptors in mammalian cells, a mechanistic understanding of the underlying molecular processes is still puzzling. In an attempt to fill this gap we followed a combined experimental and theoretical approach by dissecting and modulating cargo binding, internalization and cellular response mediated by KDEL receptors (KDELRs) at the mammalian cell surface after interaction with a model cargo/ligand. Using a fluorescent variant of ricin toxin A chain as KDELR-ligand, eGFP-RTA (H/KDEL), we demonstrate that cargo binding induces dose-dependent receptor cluster formation at and subsequent internalization from the membrane which is associated and counteracted by anterograde and microtubule-assisted receptor transport to preferred docking sites at the plasma membrane. By means of analytical arguments and extensive numerical simulations we show that cargo-synchronized receptor transport from and to the membrane is causative for KDELR/cargo cluster formation at the mammalian cell surface.

Induced effects by direct exposure to ionizing radiation (IR) are a central issue in many fields like radiation protection, clinic diagnosis and oncological therapies. Direct irradiation at certain doses induce cell death, but similar effects can also occur in cells no directly exposed to IR, a mechanism known as bystander effect. Non-IR (radiofrequency waves) can induce the death of cells loaded with MNPs in a focused oncological therapy known as magnetic hyperthermia. Indirect mechanisms are also able to induce the death of unloaded MNPs cells. Using in vitro cell models, we found that colocalization of the MNPs at the lysosomes and the non-increase of the temperature induces bystander effect under non-IR. Our results provide a landscape in which bystander effects are a more general mechanism, up to now only observed and clinically used in the field of radiotherapy.

Eukaryotic cell development has been optimized by natural selection to obey maximal intracellular flux of messenger proteins. This, in turn, implies maximum Fisher information on angular position about a target nuclear pore complex (NPR). The cell is simply modeled as spherical, with cell membrane (CM) diameter 10 micron and concentric nuclear membrane (NM) diameter 6 micron. The NM contains about 3000 nuclear pore complexes (NPCs). Development requires messenger ligands to travel from the CM-NPC-DNA target binding sites. Ligands acquire negative charge by phosphorylation, passing through the cytoplasm over Newtonian trajectories toward positively charged NPCs (utilizing positive nuclear localization sequences). The CM-NPC channel obeys maximized mean protein flux F and Fisher information I at the NPC, with first-order delta I = 0 and approximate 2nd-order delta I = 0 stability to environmental perturbations. Many of its predictions are confirmed, including the dominance of protein pathways of from 1-4 proteins, a 4nm size for the EGFR protein and the approximate flux value F =10^16 proteins/m2-s. After entering the nucleus, each protein ultimately delivers its ligand information to a DNA target site with maximum probability, i.e. maximum Kullback-Liebler entropy HKL. In a smoothness limit HKL approaches IDNA/2, so that the total CM-NPC-DNA channel obeys maximum Fisher I. Thus maximum information approaches non-equilibrium, one condition for life.

As a cheap and safe antimalarial agent, chloroquine (CQ) has been used in the battle against malaria for more than half century. However, the mechanism of CQ action and resistance in Plasmodium falciparum remains elusive. Based on further analysis of our published experimental results, we propose that the mechanism of CQ action and resistance might be closely linked with cell-cycle-associated amplified genomic-DNA fragments (CAGFs, singular form = CAGF) as CQ induces CAGF production in P. falciparum, which could affect multiple biological processes of the parasite, and thus might contribute to parasite death and CQ resistance. Recently, we found that CQ induced one of CAGFs, UB1- CAGF, might downregulate a probable P. falciparum cystine transporter (Pfct) gene expression, which could be used to understand the mechanism of CQ action and resistance in P. falciparum.

How cells regulate the number of organelles is a fundamental question in cell biology. While decades of experimental work have uncovered four fundamental processes that regulate organelle biogenesis, namely, de novo synthesis, fission, fusion and decay, a comprehensive understanding of how these processes together control organelle abundance remains elusive. Recent fluorescence microscopy experiments allow for the counting of organelles at the single-cell level. These measurements provide information about the cell-to-cell variability in organelle abundance in addition to the mean level. Motivated by such measurements, we build upon a recent study and analyze a general stochastic model of organelle biogenesis. We compute the exact analytical expressions for the probability distribution of organelle numbers, their mean, and variance across a population of single cells. It is shown that different mechanisms of organelle biogenesis lead to distinct signatures in the distribution of organelle numbers which allows us to discriminate between these various mechanisms. By comparing our theory against published data for peroxisome abundance measurements in yeast, we show that a widely believed model of peroxisome biogenesis that involves de novo synthesis, fission, and decay is inadequate in explaining the data. Also, our theory predicts bimodality in certain limits of the model. Overall, the framework developed here can be harnessed to gain mechanistic insights into the process of organelle biogenesis.