Bistability versus Bimodal Distributions in Gene Regulatory Processes from Population Balance

In recent times, stochastic treatments of gene regulatory processes have appeared in the literature in which a cell exposed to a signaling molecule in its environment triggers the synthesis of a specific protein through a network of intracellular reactions. The stochastic nature of this process leads to a distribution of protein levels in a population of cells as determined by a Fokker-Planck equation. Often instability occurs as a consequence of two (stable) steady state protein levels, one at the low end representing the “off” state, and the other at the high end representing the “on” state for a given concentration of the signaling molecule within a suitable range. A consequence of such bistability has been the appearance of bimodal distributions indicating two different populations, one in the “off” state and the other in the “on” state. The bimodal distribution can come about from stochastic analysis of a single cell. However, the concerted action of the population altering the extracellular concentration in the environment of individual cells and hence their behavior can only be accomplished by an appropriate population balance model which accounts for the reciprocal effects of interaction between the population and its environment. In this study, we show how to formulate a population balance model in which stochastic gene expression in individual cells is incorporated. Interestingly, the simulation of the model shows that bistability is neither sufficient nor necessary for bimodal distributions in a population. The original notion of linking bistability with bimodal distribution from single cell stochastic model is therefore only a special consequence of a population balance model

Multi-Scale Modeling of HIV Infection in vitro and APOBEC3G-Based Anti-Retroviral Therapy

The human APOBEC3G is an innate restriction factor that, in the absence of Vif, restricts HIV-1 replication by inducing excessive deamination of cytidine residues in nascent reverse transcripts and inhibiting reverse transcription and integration. To shed light on impact of A3G-Vif interactions on HIV replication, we developed a multi-scale computational system consisting of intracellular (single-cell), cellular and extracellular (multicellular) events by using ordinary differential equations. The single-cell model describes molecular-level events within individual cells (such as production and degradation of host and viral proteins, and assembly and release of new virions), whereas the multicellular model describes the viral dynamics and multiple cycles of infection within a population of cells. We estimated the model parameters either directly from previously published experimental data or by running simulations to find the optimum values. We validated our integrated model by reproducing the results of in vitro T cell culture experiments. Crucially, both downstream effects of A3G (hypermutation and reduction of viral burst size) were necessary to replicate the experimental results in silico. We also used the model to study anti-HIV capability of several possible therapeutic strategies including: an antibody to Vif; upregulation of A3G; and mutated forms of A3G. According to our simulations, A3G with a mutated Vif binding site is predicted to be significantly more effective than other molecules at the same dose. Ultimately, we performed sensitivity analysis to identify important model parameters. The results showed that the timing of particle formation and virus release had the highest impacts on HIV replication. The model also predicted that the degradation of A3G by Vif is not a crucial step in HIV pathogenesis