Preservation of dynamic properties in qualitative modeling frameworks for gene regulatory networks

Mathematical modeling often helps to provide a systems perspective on gene regulatory networks. In particular, qualitative approaches are useful when detailed kinetic information is lacking. Multiple methods have been developed that implement qualitative information in different ways, e.g., in purely discrete or hybrid discrete/continuous models. In this paper, we compare the discrete asynchronous logical modeling formalism for gene regulatory networks due to R. Thomas with piecewise affine differential equation models. We provide a local characterization of the qualitative dynamics of a piecewise affine differential equation model using the discrete dynamics of a corresponding Thomas model. Based on this result, we investigate the consistency of higher-level dynamical properties such as attractor characteristics and reachability. We show that although the two approaches are based on equivalent information, the resulting qualitative dynamics are different. In particular, the dynamics of the piecewise affine differential equation model is not a simple refinement of the dynamics of the Thomas model     

An integrative and practical evolutionary optimization for a complex, dynamic model of biological networks

Computer simulation is an important technique to capture the dynamics of biochemical networks. Numerical optimization is the key to estimate the values of kinetic parameters so that the dynamic model reproduces the behaviors of the existing experimental data. It is required to develop general strategies for the optimization of complex biochemical networks with a huge space of search parameters, under the condition that kinetic and quantitative data are hardly available. We propose an integrative and practical strategy for optimizing a complex dynamic model by using qualitative and incomplete experimental data. The key technologies are the divide and conquer method for reducing the search space, handling of multiple objective functions representing different types of biological behaviors, and design of rule-based objective functions that are suitable for qualitative and error-prone experimental data. This strategy is applied to optimizing a dynamic model of the yeast cell cycle to demonstrate the feasibility of it.     

A kinetic model describes metabolic response to perturbations and distribution of flux control in the benzenoid network of Petunia hybrida

In recent years there has been much interest in the genetic enhancement of plant metabolism; however, attempts at genetic modification are often unsuccessful due to an incomplete understanding of network dynamics and their regulatory properties. Kinetic modeling of plant metabolic networks can provide predictive information on network control and response to genetic perturbations, which allow estimation of flux at any concentration of intermediate or enzyme in the system. In this research, a kinetic model of the benzenoid network was developed to simulate whole network responses to different concentrations of supplied phenylalanine (Phe) in petunia flowers and capture flux redistributions caused by genetic manipulations. Kinetic parameters were obtained by network decomposition and non‐linear least squares optimization of data from petunia flowers supplied with either 75 or 150 mm ²H₅‐Phe. A single set of kinetic parameters simultaneously accommodated labeling and pool size data obtained for all endogenous and emitted volatiles at the two concentrations of supplied ²H₅‐Phe. The generated kinetic model was validated using flowers from transgenic petunia plants in which benzyl CoA:benzyl alcohol/phenylethanol benzoyltransferase (BPBT) was down‐regulated via RNAi. The determined in vivo kinetic parameters were used for metabolic control analysis, in which flux control coefficients were calculated for fluxes around the key branch point at Phe and revealed that phenylacetaldehyde synthase activity is the primary controlling factor for the phenylacetaldehyde branch of the benzenoid network. In contrast, control of flux through the β‐oxidative and non‐β‐oxidative pathways is highly distributed.     

Novel recurrent neural network for modelling biological networks: Oscillatory p53 interaction dynamics

Understanding the control of cellular networks consisting of gene and protein interactions and their emergent properties is a central activity of Systems Biology research. For this, continuous, discrete, hybrid, and stochastic methods have been proposed. Currently, the most common approach to modelling accurate temporal dynamics of networks is ordinary differential equations (ODE). However, critical limitations of ODE models are difficulty in kinetic parameter estimation and numerical solution of a large number of equations, making them more suited to smaller systems. In this article, we introduce a novel recurrent artificial neural network (RNN) that addresses above limitations and produces a continuous model that easily estimates parameters from data, can handle a large number of molecular interactions and quantifies temporal dynamics and emergent systems properties. This RNN is based on a system of ODEs representing molecular interactions in a signalling network. Each neuron represents concentration change of one molecule represented by an ODE. Weights of the RNN correspond to kinetic parameters in the system and can be adjusted incrementally during network training. The method is applied to the p53-Mdm2 oscillation system – a crucial component of the DNA damage response pathways activated by a damage signal. Simulation results indicate that the proposed RNN can successfully represent the behaviour of the p53-Mdm2 oscillation system and solve the parameter estimation problem with high accuracy. Furthermore, we presented a modified form of the RNN that estimates parameters and captures systems dynamics from sparse data collected over relatively large time steps. We also investigate the robustness of the p53-Mdm2 system using the trained RNN under various levels of parameter perturbation to gain a greater understanding of the control of the p53-Mdm2 system. Its outcomes on robustness are consistent with the current biological knowledge of this system. As more quantitative data become available on individual proteins, the RNN would be able to refine parameter estimation and mapping of temporal dynamics of individual signalling molecules as well as signalling networks as a system. Moreover, RNN can be used to modularise large signalling networks.     

The Small Heat Shock Protein Hsp27 Affects Assembly Dynamics and?Structure of Keratin Intermediate Filament Networks

The mechanical properties of living cells are essential for many processes. They are defined by the cytoskeleton, a composite network of protein fibers. Thus, the precise control of its architecture is of paramount importance. Our knowledge about the molecular and physical mechanisms defining the network structure remains scarce, especially for the intermediate filament cytoskeleton. Here, we investigate the effect of small heat shock proteins on the keratin 8/18 intermediate filament cytoskeleton using a well-controlled model system of reconstituted keratin networks. We demonstrate that Hsp27 severely alters the structure of such networks by changing their assembly dynamics. Furthermore, the C-terminal tail domain of keratin 8 is shown to be essential for this effect. Combining results from fluorescence and electron microscopy with data from analytical ultracentrifugation reveals the crucial role of kinetic trapping in keratin network formation.     

Network analysis of protein dynamics

The network paradigm is increasingly used to describe the topology and dynamics of complex systems. Here, we review the results of the topological analysis of protein structures as molecular networks describing their small-world character, and the role of hubs and central network elements in governing enzyme activity, allosteric regulation, protein motor function, signal transduction and protein stability. We summarize available data how central network elements are enriched in active centers and ligand binding sites directing the dynamics of the entire protein. We assess the feasibility of conformational and energy networks to simplify the vast complexity of rugged energy landscapes and to predict protein folding and dynamics. Finally, we suggest that modular analysis, novel centrality measures, hierarchical representation of networks and the analysis of network dynamics will soon lead to an expansion of this field.     

Multiple sources and routes of information transmission: Implications for epidemic dynamics

In a recent paper [20], we proposed and analyzed a compartmental ODE-based model describing the dynamics of an infectious disease where the presence of the pathogen also triggers the diffusion of information about the disease. In this paper, we extend this previous work by presenting results based on pairwise and simulation models that are better suited for capturing the population contact structure at a local level. We use the pairwise model to examine the potential of different information generating mechanisms and routes of information transmission to stop disease spread or to minimize the impact of an epidemic. The individual-based simulation is used to better differentiate between the networks of disease and information transmission and to investigate the impact of different basic network topologies and network overlap on epidemic dynamics. The paper concludes with an individual-based semi-analytic calculation of R₀ at the non-trivial disease free equilibrium.     

The Small Heat Shock Protein Hsp27 Affects Assembly Dynamics and Structure of Keratin Intermediate Filament Networks

The mechanical properties of living cells are essential for many processes. They are defined by the cytoskeleton, a composite network of protein fibers. Thus, the precise control of its architecture is of paramount importance. Our knowledge about the molecular and physical mechanisms defining the network structure remains scarce, especially for the intermediate filament cytoskeleton. Here, we investigate the effect of small heat shock proteins on the keratin 8/18 intermediate filament cytoskeleton using a well-controlled model system of reconstituted keratin networks. We demonstrate that Hsp27 severely alters the structure of such networks by changing their assembly dynamics. Furthermore, the C-terminal tail domain of keratin 8 is shown to be essential for this effect. Combining results from fluorescence and electron microscopy with data from analytical ultracentrifugation reveals the crucial role of kinetic trapping in keratin network formation.     

Dynamic proteomics in modeling of the living cell. Protein-protein interactions

This review is devoted to describing, summarizing, and analyzing of dynamic proteomics data obtained over the last few years and concerning the role of protein-protein interactions in modeling of the living cell. Principles of modern high-throughput experimental methods for investigation of protein-protein interactions are described. Systems biology approaches based on integrative view on cellular processes are used to analyze organization of protein interaction networks. It is proposed that finding of some proteins in different protein complexes can be explained by their multi-modular and polyfunctional properties; the different protein modules can be located in the nodes of protein interaction networks. Mathematical and computational approaches to modeling of the living cell with emphasis on molecular dynamics simulation are provided. The role of the network analysis in fundamental medicine is also briefly reviewed.     

The organisational structure of protein networks: revisiting the centrality–lethality hypothesis

Protein networks, describing physical interactions as well as functional associations between proteins, have been unravelled for many organisms in the recent past. Databases such as the STRING provide excellent resources for the analysis of such networks. In this contribution, we revisit the organisation of protein networks, particularly the centrality–lethality hypothesis, which hypothesises that nodes with higher centrality in a network are more likely to produce lethal phenotypes on removal, compared to nodes with lower centrality. We consider the protein networks of a diverse set of 20 organisms, with essentiality information available in the Database of Essential Genes and assess the relationship between centrality measures and lethality. For each of these organisms, we obtained networks of high-confidence interactions from the STRING database, and computed network parameters such as degree, betweenness centrality, closeness centrality and pairwise disconnectivity indices. We observe that the networks considered here are predominantly disassortative. Further, we observe that essential nodes in a network have a significantly higher average degree and betweenness centrality, compared to the network average. Most previous studies have evaluated the centrality–lethality hypothesis for Saccharomyces cerevisiae and Escherichia coli; we here observe that the centrality–lethality hypothesis hold goods for a large number of organisms, with certain limitations. Betweenness centrality may also be a useful measure to identify essential nodes, but measures like closeness centrality and pairwise disconnectivity are not significantly higher for essential nodes.     

How to identify essential genes from molecular networks?

Background  

The prediction of essential genes from molecular networks is a way to test the understanding of essentiality in the context of what is known about the network. However, the current knowledge on molecular network structures is incomplete yet, and consequently the strategies aimed to predict essential genes are prone to uncertain predictions. We propose that simultaneously evaluating different network structures and different algorithms representing gene essentiality (centrality measures) may identify essential genes in networks in a reliable fashion.     

Signal transduction networks: topology, response and biochemical processes

Conventionally, biological signal transduction networks are analysed using experimental and theoretical methods to describe specific protein components, interactions, and biochemical processes and to model network behavior under various conditions. While these studies provide crucial information on specific networks, this information is not easily converted to a broader understanding of signal transduction systems. Here, using a specific model of protein interaction we analyse small network topologies to understand their response and general properties. In particular, we catalogue the response for all possible topologies of a given network size to generate a response distribution, analyse the effects of specific biochemical processes on this distribution, and analyse the robustness and diversity of responses with respect to internal fluctuations or mutations in the network. The results show that even three- and four-protein networks are capable of creating diverse and biologically relevant responses, that the distribution of response types changes drastically as a function of biochemical processes at protein level, and that certain topologies strongly pre-dispose a specific response type while others allow for diverse types of responses. This study sheds light on the response types and properties that could be expected from signal transduction networks, provides possible explanations for the role of certain biochemical processes in signal transduction and suggests novel approaches to interfere with signaling pathways at the molecular level. Furthermore it shows that network topology plays a key role on determining response type and properties and that proper representation of network topology is crucial to discover and understand so-called building blocks of large networks.     

iSCHRUNK – In Silico Approach to Characterization and Reduction of Uncertainty in the Kinetic Models of Genome-scale Metabolic Networks

Accurate determination of physiological states of cellular metabolism requires detailed information about metabolic fluxes, metabolite concentrations and distribution of enzyme states. Integration of fluxomics and metabolomics data, and thermodynamics-based metabolic flux analysis contribute to improved understanding of steady-state properties of metabolism. However, knowledge about kinetics and enzyme activities though essential for quantitative understanding of metabolic dynamics remains scarce and involves uncertainty. Here, we present a computational methodology that allow us to determine and quantify the kinetic parameters that correspond to a certain physiology as it is described by a given metabolic flux profile and a given metabolite concentration vector. Though we initially determine kinetic parameters that involve a high degree of uncertainty, through the use of kinetic modeling and machine learning principles we are able to obtain more accurate ranges of kinetic parameters, and hence we are able to reduce the uncertainty in the model analysis. We computed the distribution of kinetic parameters for glucose-fed E. coli producing 1,4-butanediol and we discovered that the observed physiological state corresponds to a narrow range of kinetic parameters of only a few enzymes, whereas the kinetic parameters of other enzymes can vary widely. Furthermore, this analysis suggests which are the enzymes that should be manipulated in order to engineer the reference state of the cell in a desired way. The proposed approach also sets up the foundations of a novel type of approaches for efficient, non-asymptotic, uniform sampling of solution spaces.     

A mathematical model to study the effects of drugs administration on tumor growth dynamics

A mathematical model for describing the cancer growth dynamics in response to anticancer agents administration in xenograft models is discussed. The model consists of a system of ordinary differential equations involving five parameters (three for describing the untreated growth and two for describing the drug action). Tumor growth in untreated animals is modelled by an exponential growth followed by a linear growth. In treated animals, tumor growth rate is decreased by an additional factor proportional to both drug concentration and proliferating cells. The mathematical analysis conducted in this paper highlights several interesting properties of this tumor growth model. It suggests also effective strategies to design in vivo experiments in animals with potential saving of time and resources. For example, the drug concentration threshold for the tumor eradication, the delay between drug administration and tumor regression, and a time index that measures the efficacy of a treatment are derived and discussed. The model has already been employed in several drug discovery projects. Its application on a data set coming from one of these projects is discussed in this paper.     

The signaling petri net-based simulator: a non-parametric strategy for characterizing the dynamics of cell-specific signaling networks

Reconstructing cellular signaling networks and understanding how they work are major endeavors in cell biology. The scale and complexity of these networks, however, render their analysis using experimental biology approaches alone very challenging. As a result, computational methods have been developed and combined with experimental biology approaches, producing powerful tools for the analysis of these networks. These computational methods mostly fall on either end of a spectrum of model parameterization. On one end is a class of structural network analysis methods; these typically use the network connectivity alone to generate hypotheses about global properties. On the other end is a class of dynamic network analysis methods; these use, in addition to the connectivity, kinetic parameters of the biochemical reactions to predict the network's dynamic behavior. These predictions provide detailed insights into the properties that determine aspects of the network's structure and behavior. However, the difficulty of obtaining numerical values of kinetic parameters is widely recognized to limit the applicability of this latter class of methods. Several researchers have observed that the connectivity of a network alone can provide significant insights into its dynamics. Motivated by this fundamental observation, we present the signaling Petri net, a non-parametric model of cellular signaling networks, and the signaling Petri net-based simulator, a Petri net execution strategy for characterizing the dynamics of signal flow through a signaling network using token distribution and sampling. The result is a very fast method, which can analyze large-scale networks, and provide insights into the trends of molecules' activity-levels in response to an external stimulus, based solely on the network's connectivity. We have implemented the signaling Petri net-based simulator in the PathwayOracle toolkit, which is publicly available at http://bioinfo.cs.rice.edu/pathwayoracle. Using this method, we studied a MAPK1,2 and AKT signaling network downstream from EGFR in two breast tumor cell lines. We analyzed, both experimentally and computationally, the activity level of several molecules in response to a targeted manipulation of TSC2 and mTOR-Raptor. The results from our method agreed with experimental results in greater than 90% of the cases considered, and in those where they did not agree, our approach provided valuable insights into discrepancies between known network connectivities and experimental observations.