Nucleosome sliding can influence the spreading of histone modifications.

Nucleosomes are the fundamental building blocks of chromatin that not only help in the folding of chromatin, but also in carrying epigenetic information. It is known that nucleosome sliding is responsible for dynamically organizing chromatin structure and the resulting gene regulation. Since sliding can move two neighboring nucleosomes physically close or away, can it play a role in the spreading of histone modifications? We investigate this by simulating a stochastic model that couples nucleosome dynamics with the kinetics of histone modifications. We show that the sliding of nucleosomes can affect the modification pattern as well as the time it takes to modify a given region of chromatin. Exploring different nucleosome densities and modification kinetic parameters, we show that nucleosome sliding can be important for creating histone modification domains. Our model predicts that nucleosome density coupled with sliding dynamics can create an asymmetric histone modification profile around regulatory regions. We also compute the probability distribution of modified nucleosomes and relaxation kinetics of modifications. Our predictions are comparable with known experimental results.

Continuous variable responses and signal gating form kinetic bases for pulsatile insulin signaling and emergence of resistance.

Understanding kinetic control of biological processes is as important as identifying components that constitute pathways. Insulin signaling is central for almost all metazoans, and its perturbations are associated with various developmental disorders, metabolic diseases, and aging. While temporal phosphorylation changes and kinetic constants have provided some insights, constant or variable parameters that establish and maintain signal topology are poorly understood. Here, we report kinetic parameters that encode insulin concentration and nutrient-dependent flow of information using iterative experimental and mathematical simulation-based approaches. Our results illustrate how dynamics of distinct phosphorylation events collectively contribute to selective kinetic gating of signals and maximum connectivity of the signaling cascade under normo-insulinemic but not hyper-insulinemic states. In addition to identifying parameters that provide predictive value for maintaining the balance between metabolic and growth-factor arms, we posit a kinetic basis for the emergence of insulin resistance. Given that pulsatile insulin secretion during a fasted state precedes a fed response, our findings reveal rewiring of insulin signaling akin to memory and anticipation, which was hitherto unknown. Striking disparate temporal behavior of key phosphorylation events that destroy the topology under hyper-insulinemic states underscores the importance of unraveling regulatory components that act as bandwidth filters. In conclusion, besides providing fundamental insights, our study will help in identifying therapeutic strategies that conserve coupling between metabolic and growth-factor arms, which is lost in diseases and conditions of hyper-insulinemia.

Solving the problem of negative populations in approximate accelerated stochastic simulations using the representative reaction approach.

Methods based on the stochastic formulation of chemical kinetics have the potential to accurately reproduce the dynamical behavior of various biochemical systems of interest. However, the computational expense makes them impractical for the study of real systems. Attempts to render these methods practical have led to the development of accelerated methods, where the reaction numbers are modeled by Poisson random numbers. However, for certain systems, such methods give rise to physically unrealistic negative numbers for species populations. The methods which make use of binomial variables, in place of Poisson random numbers, have since become popular, and have been partially successful in addressing this problem. In this manuscript, the development of two new computational methods, based on the representative reaction approach (RRA), has been discussed. The new methods endeavor to solve the problem of negative numbers, by making use of tools like the stochastic simulation algorithm and the binomial method, in conjunction with the RRA. It is found that these newly developed methods perform better than other binomial methods used for stochastic simulations, in resolving the problem of negative populations.

A new approximate method for the stochastic simulation of chemical systems: the representative reaction approach.

We have developed two new approximate methods for stochastically simulating chemical systems. The methods are based on the idea of representing all the reactions in the chemical system by a single reaction, i.e., by the "representative reaction approach" (RRA). Discussed in the article are the concepts underlying the new methods along with flowchart with all the steps required for their implementation. It is shown that the two RRA methods {with the reaction 2A −> B as the representative reaction (RR)} perform creditably with regard to accuracy and computational efficiency, in comparison to the exact stochastic simulation algorithm (SSA) developed by Gillespie and are able to successfully reproduce at least the first two moments of the probability distribution of each species in the systems studied. As such, the RRA methods represent a promising new approach for stochastically simulating chemical systems.

A theoretical study of metal-metal cooperativity in the homogeneous water gas shift reaction.

The possibility of metal-metal cooperativity in improving the yield of the homogeneous water gas shift reaction (WGSR) has been investigated through full quantum mechanical density functional theory calculations. The calculations indicate that bimetallic catalysts would be likely to be more highly active than mononuclear metal-based catalysts for the WGSR. The results have implications for the design of improved WGSR catalysts in the future.