Intervention in gene regulatory networks via greedy control policies based on long-run behavior

<p>Abstract</p> <p>Background</p> <p>A salient purpose for studying gene regulatory networks is to derive intervention strategies, the goals being to identify potential drug targets and design gene-based therapeutic intervention. Optimal stochastic control based on the transition probability matrix of the underlying Markov chain has been studied extensively for probabilistic Boolean networks. Optimization is based on minimization of a cost function and a key goal of control is to reduce the steady-state probability mass of undesirable network states. Owing to computational complexity, it is difficult to apply optimal control for large networks.</p> <p>Results</p> <p>In this paper, we propose three new greedy stationary control policies by directly investigating the effects on the network long-run behavior. Similar to the recently proposed mean-first-passage-time (MFPT) control policy, these policies do not depend on minimization of a cost function and avoid the computational burden of dynamic programming. They can be used to design stationary control policies that avoid the need for a user-defined cost function because they are based directly on long-run network behavior; they can be used as an alternative to dynamic programming algorithms when the latter are computationally prohibitive; and they can be used to predict the best control gene with reduced computational complexity, even when one is employing dynamic programming to derive the final control policy. We compare the performance of these three greedy control policies and the MFPT policy using randomly generated probabilistic Boolean networks and give a preliminary example for intervening in a mammalian cell cycle network.</p> <p>Conclusion</p> <p>The newly proposed control policies have better performance in general than the MFPT policy and, as indicated by the results on the mammalian cell cycle network, they can potentially serve as future gene therapeutic intervention strategies.</p

Filling Kinetic Gaps: Dynamic Modeling of Metabolism Where Detailed Kinetic Information Is Lacking

Integrative analysis between dynamical modeling of metabolic networks and data obtained from high throughput technology represents a worthy effort toward a holistic understanding of the link among phenotype and dynamical response. Even though the theoretical foundation for modeling metabolic network has been extensively treated elsewhere, the lack of kinetic information has limited the analysis in most of the cases. To overcome this constraint, we present and illustrate a new statistical approach that has two purposes: integrate high throughput data and survey the general dynamical mechanisms emerging for a slightly perturbed metabolic network.This paper presents a statistic framework capable to study how and how fast the metabolites participating in a perturbed metabolic network reach a steady-state. Instead of requiring accurate kinetic information, this approach uses high throughput metabolome technology to define a feasible kinetic library, which constitutes the base for identifying, statistical and dynamical properties during the relaxation. For the sake of illustration we have applied this approach to the human Red blood cell metabolism (hRBC) and its capacity to predict temporal phenomena was evaluated. Remarkable, the main dynamical properties obtained from a detailed kinetic model in hRBC were recovered by our statistical approach. Furthermore, robust properties in time scale and metabolite organization were identify and one concluded that they are a consequence of the combined performance of redundancies and variability in metabolite participation.In this work we present an approach that integrates high throughput metabolome data to define the dynamic behavior of a slightly perturbed metabolic network where kinetic information is lacking. Having information of metabolite concentrations at steady-state, this method has significant relevance due its potential scope to analyze others genome scale metabolic reconstructions. Thus, I expect this approach will significantly contribute to explore the relationship between dynamic and physiology in other metabolic reconstructions, particularly those whose kinetic information is practically nulls. For instances, I envisage that this approach can be useful in genomic medicine or pharmacogenomics, where the estimation of time scales and the identification of metabolite organization may be crucial to characterize and identify (dis)functional stages

A Viral Dynamic Model for Treatment Regimens with Direct-acting Antivirals for Chronic Hepatitis C Infection

We propose an integrative, mechanistic model that integrates in vitro virology data, pharmacokinetics, and viral response to a combination regimen of a direct-acting antiviral (telaprevir, an HCV NS3-4A protease inhibitor) and peginterferon alfa-2a/ribavirin (PR) in patients with genotype 1 chronic hepatitis C (CHC). This model, which was parameterized with on-treatment data from early phase clinical studies in treatment-naïve patients, prospectively predicted sustained virologic response (SVR) rates that were comparable to observed rates in subsequent clinical trials of regimens with different treatment durations in treatment-naïve and treatment-experienced populations. The model explains the clinically-observed responses, taking into account the IC50, fitness, and prevalence prior to treatment of viral resistant variants and patient diversity in treatment responses, which result in different eradication times of each variant. The proposed model provides a framework to optimize treatment strategies and to integrate multifaceted mechanistic information and give insight into novel CHC treatments that include direct-acting antiviral agents

A cellular defense memory imprinted by early life toxic stress

Stress exposure early in life is implicated in various behavioural and somatic diseases. Experiences during the critical perinatal period form permanent, imprinted memories promoting adult survival. Although imprinting is widely recognized to dictate behaviour, whether it actuates specific transcriptional responses at the cellular level is unknown. Here we report that in response to early life stresses, Caenorhabditis elegans nematodes form an imprinted cellular defense memory. We show that exposing newly-born worms to toxic antimycin A and paraquat, respectively, stimulates the expression of toxin-specific cytoprotective reporters. Toxin exposure also induces avoidance of the toxin-containing bacterial lawn. In contrast, adult worms do not exhibit aversive behaviour towards stress-associated bacterial sensory cues. However, the mere re-encounter with the same cues reactivates the previously induced cytoprotective reporters. Learned adult defenses require memory formation during the L1 larval stage and do not appear to confer increased protection against the toxin. Thus, exposure of C. elegans to toxic stresses in the critical period elicits adaptive behavioural and cytoprotective responses, which do not form imprinted aversive behaviour, but imprint a cytoprotective memory. Our findings identify a novel form of imprinting and suggest that imprinted molecular defenses might underlie various pathophysiological alterations related to early life stress. © 2019, The Author(s)

Integrating Quantitative Knowledge into a Qualitative Gene Regulatory Network

Despite recent improvements in molecular techniques, biological knowledge remains incomplete. Any theorizing about living systems is therefore necessarily based on the use of heterogeneous and partial information. Much current research has focused successfully on the qualitative behaviors of macromolecular networks. Nonetheless, it is not capable of taking into account available quantitative information such as time-series protein concentration variations. The present work proposes a probabilistic modeling framework that integrates both kinds of information. Average case analysis methods are used in combination with Markov chains to link qualitative information about transcriptional regulations to quantitative information about protein concentrations. The approach is illustrated by modeling the carbon starvation response in Escherichia coli. It accurately predicts the quantitative time-series evolution of several protein concentrations using only knowledge of discrete gene interactions and a small number of quantitative observations on a single protein concentration. From this, the modeling technique also derives a ranking of interactions with respect to their importance during the experiment considered. Such a classification is confirmed by the literature. Therefore, our method is principally novel in that it allows (i) a hybrid model that integrates both qualitative discrete model and quantities to be built, even using a small amount of quantitative information, (ii) new quantitative predictions to be derived, (iii) the robustness and relevance of interactions with respect to phenotypic criteria to be precisely quantified, and (iv) the key features of the model to be extracted that can be used as a guidance to design future experiments