Dynamic Modelling under Uncertainty: The Case of Trypanosoma brucei Energy Metabolism

Kinetic models of metabolism require detailed knowledge of kinetic parameters. However, due to measurement errors or lack of data this knowledge is often uncertain. The model of glycolysis in the parasitic protozoan Trypanosoma brucei is a particularly well analysed example of a quantitative metabolic model, but so far it has been studied with a fixed set of parameters only. Here we evaluate the effect of parameter uncertainty. In order to define probability distributions for each parameter, information about the experimental sources and confidence intervals for all parameters were collected. We created a wiki-based website dedicated to the detailed documentation of this information: the SilicoTryp wiki (http://silicotryp.ibls.gla.ac.uk/wiki/Glycolysis). Using information collected in the wiki, we then assigned probability distributions to all parameters of the model. This allowed us to sample sets of alternative models, accurately representing our degree of uncertainty. Some properties of the model, such as the repartition of the glycolytic flux between the glycerol and pyruvate producing branches, are robust to these uncertainties. However, our analysis also allowed us to identify fragilities of the model leading to the accumulation of 3-phosphoglycerate and/or pyruvate. The analysis of the control coefficients revealed the importance of taking into account the uncertainties about the parameters, as the ranking of the reactions can be greatly affected. This work will now form the basis for a comprehensive Bayesian analysis and extension of the model considering alternative topologies

Targeting pathogen metabolism without collateral damage to the host

The development of drugs that can inactivate disease-causing cells (e.g. cancer cells or parasites) without causing collateral damage to healthy or to host cells is complicated by the fact that many proteins are very similar between organisms. Nevertheless, due to subtle, quantitative differences between the biochemical reaction networks of target cell and host, a drug can limit the flux of the same essential process in one organism more than in another. We identified precise criteria for this â €network-based' drug selectivity, which can serve as an alternative or additive to structural differences. We combined computational and experimental approaches to compare energy metabolism in the causative agent of sleeping sickness, Trypanosoma brucei, with that of human erythrocytes, and identified glucose transport and glyceraldehyde-3-phosphate dehydrogenase as the most selective antiparasitic targets. Computational predictions were validated experimentally in a novel parasite-erythrocytes co-culture system. Glucose-transport inhibitors killed trypanosomes without killing erythrocytes, neurons or liver cells

Two Component Systems: Physiological Effect of a Third Component

Signal transduction systems mediate the response and adaptation of organisms to environmental changes. In prokaryotes, this signal transduction is often done through Two Component Systems (TCS). These TCS are phosphotransfer protein cascades, and in their prototypical form they are composed by a kinase that senses the environmental signals (SK) and by a response regulator (RR) that regulates the cellular response. This basic motif can be modified by the addition of a third protein that interacts either with the SK or the RR in a way that could change the dynamic response of the TCS module. In this work we aim at understanding the effect of such an additional protein (which we call “third component”) on the functional properties of a prototypical TCS. To do so we build mathematical models of TCS with alternative designs for their interaction with that third component. These mathematical models are analyzed in order to identify the differences in dynamic behavior inherent to each design, with respect to functionally relevant properties such as sensitivity to changes in either the parameter values or the molecular concentrations, temporal responsiveness, possibility of multiple steady states, or stochastic fluctuations in the system. The differences are then correlated to the physiological requirements that impinge on the functioning of the TCS. This analysis sheds light on both, the dynamic behavior of synthetically designed TCS, and the conditions under which natural selection might favor each of the designs. We find that a third component that modulates SK activity increases the parameter space where a bistable response of the TCS module to signals is possible, if SK is monofunctional, but decreases it when the SK is bifunctional. The presence of a third component that modulates RR activity decreases the parameter space where a bistable response of the TCS module to signals is possible

Competitive and Cooperative Interactions Mediate RNA Transfer from Herpesvirus Saimiri ORF57 to the Mammalian Export Adaptor ALYREF

The essential herpesvirus adaptor protein HVS ORF57, which has homologs in all other herpesviruses, promotes viral mRNA export by utilizing the cellular mRNA export machinery. ORF57 protein specifically recognizes viral mRNA transcripts, and binds to proteins of the cellular transcription-export (TREX) complex, in particular ALYREF. This interaction introduces viral mRNA to the NXF1 pathway, subsequently directing it to the nuclear pore for export to the cytoplasm. Here we have used a range of techniques to reveal the sites for direct contact between RNA and ORF57 in the absence and presence of ALYREF. A binding site within ORF57 was characterized which recognizes specific viral mRNA motifs. When ALYREF is present, part of this ORF57 RNA binding site, composed of an a-helix, binds preferentially to ALYREF. This competitively displaces viral RNA from the a-helix, but contact with RNA is still maintained by a flanking region. At the same time, the flexible N-terminal domain of ALYREF comes into contact with the viral RNA, which becomes engaged in an extensive network of synergistic interactions with both ALYREF and ORF57. Transfer of RNA to ALYREF in the ternary complex, and involvement of individual ORF57 residues in RNA recognition, were confirmed by UV cross-linking and mutagenesis. The atomic-resolution structure of the ORF57-ALYREF interface was determined, which noticeably differed from the homologous ICP27-ALYREF structure. Together, the data provides the first site-specific description of how viral mRNA is locked by a herpes viral adaptor protein in complex with cellular ALYREF, giving herpesvirus access to the cellular mRNA export machinery. The NMR strategy used may be more generally applicable to the study of fuzzy protein-protein-RNA complexes which involve flexible polypeptide regions

Multi-state Modeling of Biomolecules

Multi-state modeling of biomolecules refers to a series of techniques used to represent and compute the behavior of biological molecules or complexes that can adopt a large number of possible functional states. Biological signaling systems often rely on complexes of biological macromolecules that can undergo several functionally significant modifications that are mutually compatible. Thus, they can exist in a very large number of functionally different states. Modeling such multi-state systems poses two problems: the problem of how to describe and specify a multi-state system (the “specification problem”) and the problem of how to use a computer to simulate the progress of the system over time (the “computation problem”). To address the specification problem, modelers have in recent years moved away from explicit specification of all possible states and towards rule-based formalisms that allow for implicit model specification, including the κ-calculus [1], BioNetGen [2]–[5], the Allosteric Network Compiler [6], and others [7], [8]. To tackle the computation problem, they have turned to particle-based methods that have in many cases proved more computationally efficient than population-based methods based on ordinary differential equations, partial differential equations, or the Gillespie stochastic simulation algorithm [9], [10]. Given current computing technology, particle-based methods are sometimes the only possible option. Particle-based simulators fall into two further categories: nonspatial simulators, such as StochSim [11], DYNSTOC [12], RuleMonkey [9], [13], and the Network-Free Stochastic Simulator (NFSim) [14], and spatial simulators, including Meredys [15], SRSim [16], [17], and MCell [18]–[20]. Modelers can thus choose from a variety of tools, the best choice depending on the particular problem. Development of faster and more powerful methods is ongoing, promising the ability to simulate ever more complex signaling processes in the future

Leveraging Modeling Approaches: Reaction Networks and Rules

We have witnessed an explosive growth in research involving mathematical models and computer simulations of intracellular molecular interactions, ranging from metabolic pathways to signaling and gene regulatory networks. Many software tools have been developed to aid in the study of such biological systems, some of which have a wealth of features for model building and visualization, and powerful capabilities for simulation and data analysis. Novel high resolution and/or high throughput experimental techniques have led to an abundance of qualitative and quantitative data related to the spatio-temporal distribution of molecules and complexes, their interactions kinetics, and functional modifications. Based on this information, computational biology researchers are attempting to build larger and more detailed models. However, this has proved to be a major challenge. Traditionally, modeling tools require the explicit specification of all molecular species and interactions in a model, which can quickly become a major limitation in the case of complex networks – the number of ways biomolecules can combine to form multimolecular complexes can be combinatorially large. Recently, a new breed of software tools has been created to address the problems faced when building models marked by combinatorial complexity. These have a different approach for model specification, using reaction rules and species patterns. Here we compare the traditional modeling approach with the new rule-based methods. We make a case for combining the capabilities of conventional simulation software with the unique features and flexibility of a rule-based approach in a single software platform for building models of molecular interaction networks