Glycolysis in bloodstream form Trypanosoma brucei can be understood in terms of the kinetics of the glycolytic enzymes

The RNA-binding protein RBP42 regulates cellular energy metabolism in mammalian-infective Trypanosoma brucei

RNA-binding proteins (RBPs) are key players in coordinated post-transcriptional regulation of functionally related genes, defined as RNA regulons. RNA regulons play particularly critical roles in parasitic trypanosomes, which exhibit unregulated co-transcription of long unrelated gene arrays. In this report, we present a systematic analysis of an essential RBP, RBP42, in the mammalian-infective bloodstream form of African trypanosome and show that RBP42 is a key regulator of parasite's central carbon and energy metabolism. Using individual-nucleotide resolution UV cross-linking and immunoprecipitation to identify genome-wide RBP42-RNA interactions, we show that RBP42 preferentially binds within the coding region of mRNAs encoding core metabolic enzymes. Global quantitative transcriptomic and proteomic analyses reveal that loss of RBP42 reduces the abundance of target mRNA-encoded proteins, but not target mRNA, suggesting a positive translational regulatory role of RBP42. Significant changes in central carbon metabolic intermediates, following loss of RBP42, further support its critical role in cellular energy metabolism. Trypanosoma brucei infection, transmitted through the bite of blood-feeding tsetse flies, causes deadly diseases in humans and livestock. This disease, if left untreated, is almost always fatal. Existing therapies are toxic and difficult to administer. During T. brucei's lifecycle in two different host environments, the parasite progresses through distinctive life stages with major morphological and metabolic changes, requiring precise alteration of parasite gene expression program. In the absence of regulated transcription, post-transcriptional processes mediated by RNA-binding proteins play critical roles in T. brucei gene regulation. In this study, we show that the RNA-binding protein RBP42 plays crucial roles in cellular energy metabolic regulation of this important human pathogen. Metabolic dysregulation observed in RBP42 knockdown cells offers a breadth of potential interest to researchers studying parasite biology and can also impact research in general eukaryotic biology.

Personalised modelling of clinical heterogeneity between medium-chain acyl-CoA dehydrogenase patients

BACKGROUND

Monogenetic inborn errors of metabolism cause a wide phenotypic heterogeneity that may even differ between family members carrying the same genetic variant. Computational modelling of metabolic networks may identify putative sources of this inter-patient heterogeneity. Here, we mainly focus on medium-chain acyl-CoA dehydrogenase deficiency (MCADD), the most common inborn error of the mitochondrial fatty acid oxidation (mFAO). It is an enigma why some MCADD patients-if untreated-are at risk to develop severe metabolic decompensations, whereas others remain asymptomatic throughout life. We hypothesised that an ability to maintain an increased free mitochondrial CoA (CoASH) and pathway flux might distinguish asymptomatic from symptomatic patients.

RESULTS

We built and experimentally validated, for the first time, a kinetic model of the human liver mFAO. Metabolites were partitioned according to their water solubility between the bulk aqueous matrix and the inner membrane. Enzymes are also either membrane-bound or in the matrix. This metabolite partitioning is a novel model attribute and improved predictions. MCADD substantially reduced pathway flux and CoASH, the latter due to the sequestration of CoA as medium-chain acyl-CoA esters. Analysis of urine from MCADD patients obtained during a metabolic decompensation showed an accumulation of medium- and short-chain acylcarnitines, just like the acyl-CoA pool in the MCADD model. The model suggested some rescues that increased flux and CoASH, notably increasing short-chain acyl-CoA dehydrogenase (SCAD) levels. Proteome analysis of MCADD patient-derived fibroblasts indeed revealed elevated levels of SCAD in a patient with a clinically asymptomatic state. This is a rescue for MCADD that has not been explored before. Personalised models based on these proteomics data confirmed an increased pathway flux and CoASH in the model of an asymptomatic patient compared to those of symptomatic MCADD patients.

CONCLUSIONS

We present a detailed, validated kinetic model of mFAO in human liver, with solubility-dependent metabolite partitioning. Personalised modelling of individual patients provides a novel explanation for phenotypic heterogeneity among MCADD patients. Further development of personalised metabolic models is a promising direction to improve individualised risk assessment, management and monitoring for inborn errors of metabolism.

SBMLKinetics: a tool for annotation-independent classification of reaction kinetics for SBML models

BACKGROUND

Reaction networks are widely used as mechanistic models in systems biology to reveal principles of biological systems. Reactions are governed by kinetic laws that describe reaction rates. Selecting the appropriate kinetic laws is difficult for many modelers. There exist tools that attempt to find the correct kinetic laws based on annotations. Here, I developed annotation-independent technologies that assist modelers by focusing on finding kinetic laws commonly used for similar reactions.

RESULTS

Recommending kinetic laws and other analyses of reaction networks can be viewed as a classification problem. Existing approaches to determining similar reactions rely heavily on having good annotations, a condition that is often unsatisfied in model repositories such as BioModels. I developed an annotation-independent approach to find similar reactions via reaction classifications. I proposed a two-dimensional kinetics classification scheme (2DK) that analyzed reactions along the dimensions of kinetics type (K type) and reaction type (R type). I identified approximately ten mutually exclusive K types, including zeroth order, mass action, Michaelis-Menten, Hill kinetics, and others. R types were organized by the number of distinct reactants and the number of distinct products in reactions. I constructed a tool, SBMLKinetics, that inputted a collection of SBML models and then calculated reaction classifications as the probability of each 2DK class. The effectiveness of 2DK was evaluated on BioModels, and the scheme classified over 95% of the reactions.

CONCLUSIONS

2DK had many applications. It provided a data-driven annotation-independent approach to recommending kinetic laws by using type common for the kind of models in combination with the R type of the reactions. Alternatively, 2DK could also be used to alert users that a kinetic law was unusual for the K type and R type. Last, 2DK provided a way to analyze groups of models to compare their kinetic laws. I applied 2DK to BioModels to compare the kinetics of signaling networks with the kinetics of metabolic networks and found significant differences in K type distributions.

Delineating transitions during the evolution of specialised peroxisomes: Glycosome formation in kinetoplastid and diplonemid protists

One peculiarity of protists belonging to classes Kinetoplastea and Diplonemea within the phylum Euglenozoa is compartmentalisation of most glycolytic enzymes within peroxisomes that are hence called glycosomes. This pathway is not sequestered in peroxisomes of the third Euglenozoan class, Euglenida. Previous analysis of well-studied kinetoplastids, the ‘TriTryps’ parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania spp., identified within glycosomes other metabolic processes usually not present in peroxisomes. In addition, trypanosomatid peroxins, i.e. proteins involved in biogenesis of these organelles, are divergent from human and yeast orthologues. In recent years, genomes, transcriptomes and proteomes for a variety of euglenozoans have become available. Here, we track the possible evolution of glycosomes by querying these databases, as well as the genome of Naegleria gruberi, a non-euglenozoan, which belongs to the same protist supergroup Discoba. We searched for orthologues of TriTryps proteins involved in glycosomal metabolism and biogenesis. Predicted cellular location(s) of each metabolic enzyme identified was inferred from presence or absence of peroxisomal-targeting signals. Combined with a survey of relevant literature, we refine extensively our previously postulated hypothesis about glycosome evolution. The data agree glycolysis was compartmentalised in a common ancestor of the kinetoplastids and diplonemids, yet additionally indicates most other processes found in glycosomes of extant trypanosomatids, but not in peroxisomes of other eukaryotes were either sequestered in this ancestor or shortly after separation of the two lineages. In contrast, peroxin divergence is evident in all euglenozoans. Following their gain of pathway complexity, subsequent evolution of peroxisome/glycosome function is complex. We hypothesize compartmentalisation in glycosomes of glycolytic enzymes, their cofactors and subsequently other metabolic enzymes provided selective advantage to kinetoplastids and diplonemids during their evolution in changing marine environments. We contend two specific properties derived from the ancestral peroxisomes were key: existence of nonselective pores for small solutes and the possibility of high turnover by pexophagy. Critically, such pores and pexophagy are characterised in extant trypanosomatids. Increasing amenability of free-living kinetoplastids and recently isolated diplonemids to experimental study means our hypothesis and interpretation of bioinformatic data are suited to experimental interrogation.

Synthetic biology tools for engineering Goodwin oscillation in Trypanosoma brucei brucei

Kinetoplastid protozoa possess properties that are highly divergent from the mammalian, yeast and bacterial cells more commonly used in synthetic biology and represent a tantalisingly untapped source of bioengineering potential. Trypanosoma brucei brucei (T. b. brucei), an established model organism for studying the Kinetoplastida, is non-pathogenic to humans and provides an interesting test case for establishing synthetic biology in this phylogenetic class. To demonstrate further the tractability of Kinetoplastida to synthetic biology, we sought to construct and demonstrate a Goodwin oscillator, the simplest oscillatory gene network, in T. b. brucei for the first time. We report one completed iteration of the archetypal synthetic biology Design-Build-Test-Learn (DBTL) cycle; firstly, using Ab initio mathematical modelling of the behaviour a theoretical, oscillatory, trypanosomal synthetic gene network (SGN) to inform the design of a plasmid encoding that network. Once assembled, the plasmid was then used to generate a stable transfectant T. b. brucei cell line. To test the performance of the oscillatory SGN, a novel experimental setup was established to capture images of the fluorescent signal from motion-restricted live cells. Data captured were consistent with oscillatory behaviour of the SGN, with cellular fluorescence observed to oscillate with a period of 50 min, with varying amplitude and linear growth trend. This first DBTL cycle establishes a foundation for future cycles in which the SGN design and experimental monitoring setup can be further refined.

Protein Lactylation Critically Regulates Energy Metabolism in the Protozoan Parasite Trypanosoma brucei

Lysine lactylation has been recognized as a novel post-translational modification occurring on histones. However, lactylation in non-histone proteins, especially in proteins of early branching organisms, is not well understood. Energy metabolism and the histone repertoire in the early diverging protozoan parasite Trypanosoma brucei, the causative agent of African trypanosomiasis, markedly diverge from those of conventional eukaryotes. Here, we present the first exhaustive proteome-wide investigation of lactylated sites in T. brucei. We identified 387 lysine-lactylated sites in 257 proteins of various cellular localizations and biological functions. Further, we revealed that glucose metabolism critically regulates protein lactylation in T. brucei although the parasite lacks lactate dehydrogenase. However, unlike mammals, increasing the glucose concentration reduced the level of lactate, and protein lactylation decreased in T. brucei via a unique lactate production pathway. In addition to providing a valuable resource, these foregoing data reveal the regulatory roles of protein lactylation of trypanosomes in energy metabolism and gene expression.

Glycerol suppresses glucose consumption in trypanosomes through metabolic contest

Microorganisms must make the right choice for nutrient consumption to adapt to their changing environment. As a consequence, bacteria and yeasts have developed regulatory mechanisms involving nutrient sensing and signaling, known as "catabolite repression," allowing redirection of cell metabolism to maximize the consumption of an energy-efficient carbon source. Here, we report a new mechanism named "metabolic contest" for regulating the use of carbon sources without nutrient sensing and signaling. Trypanosoma brucei is a unicellular eukaryote transmitted by tsetse flies and causing human African trypanosomiasis, or sleeping sickness. We showed that, in contrast to most microorganisms, the insect stages of this parasite developed a preference for glycerol over glucose, with glucose consumption beginning after the depletion of glycerol present in the medium. This "metabolic contest" depends on the combination of 3 conditions: (i) the sequestration of both metabolic pathways in the same subcellular compartment, here in the peroxisomal-related organelles named glycosomes; (ii) the competition for the same substrate, here ATP, with the first enzymatic step of the glycerol and glucose metabolic pathways both being ATP-dependent (glycerol kinase and hexokinase, respectively); and (iii) an unbalanced activity between the competing enzymes, here the glycerol kinase activity being approximately 80-fold higher than the hexokinase activity. As predicted by our model, an approximately 50-fold down-regulation of the GK expression abolished the preference for glycerol over glucose, with glucose and glycerol being metabolized concomitantly. In theory, a metabolic contest could be found in any organism provided that the 3 conditions listed above are met.

Basic Biology of Trypanosoma brucei with Reference to the Development of Chemotherapies

Trypanosoma brucei are protozoan parasites that cause the lethal human disease African sleeping sickness and the economically devastating disease of cattle, Nagana. African sleeping sickness, also known as Human African Trypanosomiasis (HAT), threatens 65 million people and animal trypanosomiasis makes large areas of farmland unusable. There is no vaccine and licensed therapies against the most severe, late-stage disease are toxic, impractical and ineffective. Trypanosomes are transmitted by tsetse flies, and HAT is therefore predominantly confined to the tsetse fly belt in sub-Saharan Africa. They are exclusively extracellular and they differentiate between at least seven developmental forms that are highly adapted to host and vector niches. In the mammalian (human) host they inhabit the blood, cerebrospinal fluid (late-stage disease), skin, and adipose fat. In the tsetse fly vector they travel from the tsetse midgut to the salivary glands via the ectoperitrophic space and proventriculus. Trypanosomes are evolutionarily divergent compared with most branches of eukaryotic life. Perhaps most famous for their extraordinary mechanisms of monoallelic gene expression and antigenic variation, they have also been investigated because much of their biology is either highly unconventional or extreme. Moreover, in addition to their importance as pathogens, many researchers have been attracted to the field because trypanosomes have some of the most advanced molecular genetic tools and database resources of any model system. The following will cover just some aspects of trypanosome biology and how its divergent biochemistry has been leveraged to develop drugs to treat African sleeping sickness. This is by no means intended to be a comprehensive survey of trypanosome features. Rather, I hope to present trypanosomes as one of the most fascinating and tractable systems to do discovery biology.

Carbohydrate metabolism in trypanosomatids: New insights revealing novel complexity, diversity and species-unique features

The human pathogenic trypanosomatid species collectively called the "TriTryp parasites" - Trypanosoma brucei, Trypanosoma cruzi and Leishmania spp. - have complex life cycles, with each of these parasitic protists residing in a different niche during their successive developmental stages where they encounter diverse nutrients. Consequently, they adapt their metabolic network accordingly. Yet, throughout the life cycles, carbohydrate metabolism - involving the glycolytic, gluconeogenic and pentose-phosphate pathways - always plays a central role in the biology of these parasites, whether the available carbon and free energy sources are saccharides, amino acids or lipids. In this paper, we provide an updated review of the carbohydrate metabolism of the TriTryps, highlighting new data about this metabolic network, the interconnection of its pathways and the compartmentalisation of its enzymes within glycosomes, cytosol and mitochondrion. Differences in the expression of the branches of the metabolic network between the successive life-cycle stages of each of these parasitic trypanosomatids are discussed, as well as differences between them. Recent structural and kinetic studies have revealed unique regulatory mechanisms for some of the network's key enzymes with important species-specific variations. Furthermore, reports of multiple post-translational modifications of trypanosomal glycolytic enzymes suggest that additional mechanisms for stage- and/or environmental cues that regulate activity are operational in the parasites. The detailed comparison of the carbohydrate metabolism of the TriTryps has thus revealed multiple differences and a greater complexity, including for the reduced metabolic network in bloodstream-form T. brucei, than previously appreciated. Although these parasites are related, share many cytological and metabolic features and are grouped within a single taxonomic family, the differences highlighted in this review reflect their separate evolutionary tracks from a common ancestor to the extant organisms. These differences are indicative of their adaptation to the different insect vectors and niches occupied in their mammalian hosts.

Depletion of cardiolipin induces major changes in energy metabolism in Trypanosoma brucei bloodstream forms

The mitochondrial inner membrane glycerophospholipid cardiolipin (CL) associates with mitochondrial proteins to regulate their activities and facilitate protein complex and supercomplex formation. Loss of CL leads to destabilized respiratory complexes and mitochondrial dysfunction. The role of CL in an organism lacking a conventional electron transport chain (ETC) has not been elucidated. Trypanosoma brucei bloodstream forms use an unconventional ETC composed of glycerol-3-phosphate dehydrogenase and alternative oxidase (AOX), while the mitochondrial membrane potential (ΔΨm) is generated by the hydrolytic action of the F_o F₁ -ATP synthase (aka F_o F₁ -ATPase). We now report that the inducible depletion of cardiolipin synthase (TbCls) is essential for survival of T brucei bloodstream forms. Loss of CL caused a rapid drop in ATP levels and a decline in the ΔΨm. Unbiased proteomic analyses revealed a reduction in the levels of many mitochondrial proteins, most notably of F_o F₁ -ATPase subunits and AOX, resulting in a strong decline of glycerol-3-phosphate-stimulated oxygen consumption. The changes in cellular respiration preceded the observed decrease in F_o F₁ -ATPase stability, suggesting that the AOX-mediated ETC is the first pathway responding to the decline in CL. Select proteins and pathways involved in glucose and amino acid metabolism were upregulated to counteract the CL depletion-induced drop in cellular ATP.

Complex Stability and an Irrevertible Transition Reverted by Peptide and Fibroblasts in a Dynamic Model of Innate Immunity

We here apply a control analysis and various types of stability analysis to an in silico model of innate immunity that addresses the management of inflammation by a therapeutic peptide. Motivation is the observation, both in silico and in experiments, that this therapy is not robust. Our modeling results demonstrate how (1) the biological phenomena of acute and chronic modes of inflammation may reflect an inherently complex bistability with an irrevertible flip between the two modes, (2) the chronic mode of the model has stable, sometimes unique, steady states, while its acute-mode steady states are stable but not unique, (3) as witnessed by TNF levels, acute inflammation is controlled by multiple processes, whereas its chronic-mode inflammation is only controlled by TNF synthesis and washout, (4) only when the antigen load is close to the acute mode's flipping point, many processes impact very strongly on cells and cytokines, (5) there is no antigen exposure level below which reduction of the antigen load alone initiates a flip back to the acute mode, and (6) adding healthy fibroblasts makes the transition from acute to chronic inflammation revertible, although (7) there is a window of antigen load where such a therapy cannot be effective. This suggests that triple therapies may be essential to overcome chronic inflammation. These may comprise (1) anti-immunoglobulin light chain peptides, (2) a temporarily reduced antigen load, and (3a) fibroblast repopulation or (3b) stem cell strategies.

The kinetic characteristics of human and trypanosomatid phosphofructokinases for the reverse reaction

Eukaryotic ATP-dependent phosphofructokinases (PFKs) are often considered unidirectional enzymes catalysing the transfer of a phospho moiety from ATP to fructose 6-phosphate to produce ADP and fructose 1,6-bisphosphate. The reverse reaction is not generally considered to occur under normal conditions and has never been demonstrated for any eukaryotic ATP-dependent PFKs, though it does occur in inorganic pyrophosphate-dependent PFKs and has been experimentally shown for bacterial ATP-dependent PFKs. The evidence is provided via two orthogonal assays that all three human PFK isoforms can catalyse the reverse reaction in vitro, allowing determination of kinetic properties. Additionally, the reverse reaction was shown possible for PFKs from three clinically important trypanosomatids; these enzymes are contained within glycosomes in vivo. This compartmentalisation may facilitate reversal, given the potential for trypanosomatids to have an altered ATP/ADP ratio in glycosomes compared with the cytosol. The kinetic properties of each trypanosomatid PFK were determined, including the response to natural and artificial modulators of enzyme activity. The possible physiological relevance of the reverse reaction in trypanosomatid and human PFKs is discussed.

Glycerol supports growth of the Trypanosoma brucei bloodstream forms in the absence of glucose: Analysis of metabolic adaptations on glycerol-rich conditions

The bloodstream forms of Trypanosoma brucei (BSF), the parasite protist causing sleeping sickness, primarily proliferate in the blood of their mammalian hosts. The skin and adipose tissues were recently identified as additional major sites for parasite development. Glucose was the only carbon source known to be used by bloodstream trypanosomes to feed their central carbon metabolism, however, the metabolic behaviour of extravascular tissue-adapted parasites has not been addressed yet. Since the production of glycerol is an important primary function of adipocytes, we have adapted BSF trypanosomes to a glucose-depleted but glycerol-rich culture medium (CMM_Glyc/GlcNAc) and compared their metabolism and proteome to those of parasites grown in standard glucose-rich conditions (CMM_Glc). BSF were shown to consume 2-folds more oxygen per consumed carbon unit in CMM_Glyc/GlcNAc and were 11.5-times more sensitive to SHAM, a specific inhibitor of the plant-like alternative oxidase (TAO), which is the only mitochondrial terminal oxidase expressed in BSF. This is consistent with (i) the absolute requirement of the mitochondrial respiratory activity to convert glycerol into dihydroxyacetone phosphate, as deduced from the updated metabolic scheme and (ii) with the 1.8-fold increase of the TAO expression level compared to the presence of glucose. Proton NMR analysis of excreted end products from glycerol and glucose metabolism showed that these two carbon sources are metabolised through the same pathways, although the contributions of the acetate and succinate branches are more important in the presence of glycerol than glucose (10.2% versus 3.4% of the excreted end products, respectively). In addition, metabolomic analyses by mass spectrometry showed that, in the absence of glucose, ¹³C-labelled glycerol was incorporated into hexose phosphates through gluconeogenesis. As expected, RNAi-mediated down-regulation of glycerol kinase expression abolished glycerol metabolism and was lethal for BSF grown in CMM_Glyc/GlcNAc. Interestingly, BSF have adapted their metabolism to grow in CMM_Glyc/GlcNAc by concomitantly increasing their rate of glycerol consumption and decreasing that of glucose. However, the glycerol kinase activity was 7.8-fold lower in CMM_Glyc/GlcNAc, as confirmed by both western blotting and proteomic analyses. This suggests that the huge excess in glycerol kinase that is not absolutely required for glycerol metabolism, might be used for another yet undetermined non-essential function in glucose rich-conditions. Altogether, these data demonstrate that BSF trypanosomes are well-adapted to glycerol-rich conditions that could be encountered by the parasite in extravascular niches, such as the skin and adipose tissues.

Until very recently, the bloodstream forms (BSF) of the Trypanosoma brucei group species have been considered to propagate exclusively in the mammalian fluids, including the blood, the lymphatic network and the cerebrospinal fluid. All these fluids are rich in glucose, which is widely considered by the scientific community as the only carbon source used by the parasite to feed its central carbon metabolism and its ATP production. Here, we show for the first time that the BSF trypanosomes efficiently grow in glucose-free conditions as long as glycerol is supplied. The raison d'être of this capacity developed by BSF trypanosomes to grow in glycerol-rich conditions regardless of the glucose concentration, including in glucose-free conditions, is not yet understood. However, the recent discovery that trypanosomes colonize and proliferate in the skin and the adipose tissues of their mammalian hosts may provide a rational explanation for the development of a glycerol-based metabolism in BSF. Indeed, the adipocytes composing adipose tissues and also abundantly present in subcutaneous layers excrete large amounts of glycerol produced from the catabolism of glucose and triglycerides. We also show that BSF trypanosomes adapted to glucose-depleted conditions activate gluconeogenesis to produce the essential hexose phosphates from glycerol metabolism. Interestingly, the constitutive expression of the key gluconeogenic enzyme fructose-1,6-bisphosphatase, which is not used for glycolysis, suggests that BSF trypanosomes maintained in the standard glucose-rich medium are pre-adapted to glucose-depleted conditions. This further strengthens the new paradigm that BSF trypanosomes can use glycerol in tissues producing this carbon source, such as the skin the adipose tissues.

Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Unicellular flagellates of the family Trypanosomatidae are obligatory parasites of invertebrates, vertebrates and plants. Dixenous species are aetiological agents of a number of diseases in humans, domestic animals and plants. Their monoxenous relatives are restricted to insects. Because of the high biological diversity, adaptability to dramatically different environmental conditions, and omnipresence, these protists have major impact on all biotic communities that still needs to be fully elucidated. In addition, as these organisms represent a highly divergent evolutionary lineage, they are strikingly different from the common 'model system' eukaryotes, such as some mammals, plants or fungi. A number of excellent reviews, published over the past decade, were dedicated to specialized topics from the areas of trypanosomatid molecular and cell biology, biochemistry, host-parasite relationships or other aspects of these fascinating organisms. However, there is a need for a more comprehensive review that summarizing recent advances in the studies of trypanosomatids in the last 30 years, a task, which we tried to accomplish with the current paper.

Complexity of Biochemical and Genetic Responses Reduced Using Simple Theoretical Models

Living systems are known to behave in a complex and sometimes unpredictable manner. Humans, for a very long time, have been intrigued by nature, and have attempted to understand biological processes and mechanisms using numerous experimental and mathematical techniques. In this chapter, we will look at simple theoretical models, using both linear and nonlinear differential equations, that realistically capture complex biochemical and genetic responses of living cells. Even for cases where cellular behaviors are stochastic, as for single-cell responses, randomness added to well-defined deterministic models has elegantly been shown to be useful. The data collectively present evidence for further exploration of the self-organizing rules and laws of living matter.

The kinetic properties of hexokinases in African trypanosomes of the subgenus Trypanozoon match the blood glucose levels of mammal hosts

We hypothesize that the hexokinases of trypanosomes of the subgenus Trypanozoon match the blood glucose levels of hosts. We studied the kinetic properties of purified hexokinase in T. equiperdum (specific activity=302U/mg), and compare with other members of Trypanozoon. With ATP (K_m=104.7μM) as phosphate donor, hexokinase catalyzes the phosphorylation of glucose (K_m=24.9μM) and mannose (K_m=8.8μM). With respect to glucose, mannose and inorganic pyrophosphate respectively are a competitive, and a mixed inhibitor of hexokinase. With respect to ATP, both are mixed inhibitors of this enzyme. In T. equiperdum, hexokinase shows a high affinity for glucose. Pleomorphism-transformation of trypanosomes from a multiplicative to a non-multiplicative form-results in a self-limited growth stabilizing glucose consumption. It delays the death of the host, thus prolonging its exposure to tsetse flies. When glucose levels descend, top-down regulation allows trypanosomes to survive through the expression of alternative metabolic pathways. It accelerates the death of the host, but helps trypanosome density to increase enough to ensure transmission without tsetse flies. Pleomorphism, and a hexokinase with a high affinity for glucose, are two main adaptive traits of T. b. brucei. The latter trait, and a strong top-down regulation, are two main adaptive traits of T. equiperdum. For trypanosomes living in glucose-rich blood, a hexokinase with a high affinity for glucose would unnecessarily harm hosts. This may explain why the human parasites, T. b. gambiense and T. b. rhodesiense, possess hexokinases with a low affinity for glucose.

A computational approach to extinction events in chemical reaction networks with discrete state spaces

Recent work of Johnston et al. has produced sufficient conditions on the structure of a chemical reaction network which guarantee that the corresponding discrete state space system exhibits an extinction event. The conditions consist of a series of systems of equalities and inequalities on the edges of a modified reaction network called a domination-expanded reaction network. In this paper, we present a computational implementation of these conditions written in Python and apply the program on examples drawn from the biochemical literature. We also run the program on 458 models from the European Bioinformatics Institute's BioModels Database and report our results.

Aquaglyceroporin-null trypanosomes display glycerol transport defects and respiratory-inhibitor sensitivity

Aquaglyceroporins (AQPs) transport water and glycerol and play important roles in drug-uptake in pathogenic trypanosomatids. For example, AQP2 in the human-infectious African trypanosome, Trypanosoma brucei gambiense, is responsible for melarsoprol and pentamidine-uptake, and melarsoprol treatment-failure has been found to be due to AQP2-defects in these parasites. To further probe the roles of these transporters, we assembled a T. b. brucei strain lacking all three AQP-genes. Triple-null aqp1-2-3 T. b. brucei displayed only a very moderate growth defect in vitro, established infections in mice and recovered effectively from hypotonic-shock. The aqp1-2-3 trypanosomes did, however, display glycerol uptake and efflux defects. They failed to accumulate glycerol or to utilise glycerol as a carbon-source and displayed increased sensitivity to salicylhydroxamic acid (SHAM), octyl gallate or propyl gallate; these inhibitors of trypanosome alternative oxidase (TAO) can increase intracellular glycerol to toxic levels. Notably, disruption of AQP2 alone generated cells with glycerol transport defects. Consistent with these findings, AQP2-defective, melarsoprol-resistant clinical isolates were sensitive to the TAO inhibitors, SHAM, propyl gallate and ascofuranone, relative to melarsoprol-sensitive reference strains. We conclude that African trypanosome AQPs are dispensable for viability and osmoregulation but they make important contributions to drug-uptake, glycerol-transport and respiratory-inhibitor sensitivity. We also discuss how the AQP-dependent inverse sensitivity to melarsoprol and respiratory inhibitors described here might be exploited.

Forward design of a complex enzyme cascade reaction

Enzymatic reaction networks are unique in that one can operate a large number of reactions under the same set of conditions concomitantly in one pot, but the nonlinear kinetics of the enzymes and the resulting system complexity have so far defeated rational design processes for the construction of such complex cascade reactions. Here we demonstrate the forward design of an in vitro 10-membered system using enzymes from highly regulated biological processes such as glycolysis. For this, we adapt the characterization of the biochemical system to the needs of classical engineering systems theory: we combine online mass spectrometry and continuous system operation to apply standard system theory input functions and to use the detailed dynamic system responses to parameterize a model of sufficient quality for forward design. This allows the facile optimization of a 10-enzyme cascade reaction for fine chemical production purposes.

Building multi-component enzymatic processes in one pot is challenged by the inherent complexity of each biochemical system. Here, the authors use online mass spectroscopy and engineering systems theory to achieve forward design of a ten-membered reaction cascade.

Identification of Novel Chemical Scaffolds Inhibiting Trypanothione Synthetase from Pathogenic Trypanosomatids

Background

The search for novel chemical entities targeting essential and parasite-specific pathways is considered a priority for neglected diseases such as trypanosomiasis and leishmaniasis. The thiol-dependent redox metabolism of trypanosomatids relies on bis-glutathionylspermidine [trypanothione, T(SH)₂], a low molecular mass cosubstrate absent in the host. In pathogenic trypanosomatids, a single enzyme, trypanothione synthetase (TryS), catalyzes trypanothione biosynthesis, which is indispensable for parasite survival. Thus, TryS qualifies as an attractive drug target candidate.

Methodology/Principal Finding

A library composed of 144 compounds from 7 different families and several singletons was screened against TryS from three major pathogen species (Trypanosoma brucei, Trypanosoma cruzi and Leishmania infantum). The screening conditions were adjusted to the TryS´ kinetic parameters and intracellular concentration of substrates corresponding to each trypanosomatid species, and/or to avoid assay interference. The screening assay yielded suitable Z’ and signal to noise values (≥0.85 and ~3.5, respectively), and high intra-assay reproducibility. Several novel chemical scaffolds were identified as low μM and selective tri-tryp TryS inhibitors. Compounds displaying multi-TryS inhibition (N,N'-bis(3,4-substituted-benzyl) diamine derivatives) and an N⁵-substituted paullone (MOL2008) halted the proliferation of infective Trypanosoma brucei (EC₅₀ in the nM range) and Leishmania infantum promastigotes (EC₅₀ = 12 μM), respectively. A bis-benzyl diamine derivative and MOL2008 depleted intracellular trypanothione in treated parasites, which confirmed the on-target activity of these compounds.

Conclusions/Significance

Novel molecular scaffolds with on-target mode of action were identified as hit candidates for TryS inhibition. Due to the remarkable species-specificity exhibited by tri-tryp TryS towards the compounds, future optimization and screening campaigns should aim at designing and detecting, respectively, more potent and broad-range TryS inhibitors.

Parasites from the genus Trypanosoma and Leishmania are etiologic agents for a group of neglected diseases with high morbidity and mortality rates in the developing world. Inasmuch as vaccine development is hampered by the successful mechanisms employed by the pathogens to evade the host immune response, chemotherapy remains as a safe option to fight these diseases. However, new drugs with better pharmacological performance (i.e. safety, efficacy and ease of administration) than those in current use are urgently needed. The thiol-redox metabolism of trypanosomatids offers an excellent opportunity for the development of more selective and efficacious medicines because it depends on a molecule, trypanothione (a bis-glutathionyl derivative of spermidine), unique and indispensable to the pathogens. Here we report the identification of novel inhibitors of trypanothione synthetase from three major trypanosomatid species of medical and veterinary relevance. Although highly conserved in sequence, trypanothione synthetases display significant species-specifity towards compounds, pointing to structural differences as determinants of ligand selectivity. Most of the active compounds presented two-digit μM inhibitory activity and serve as primary scaffolds to develop more potent inhibitors. Among them, N,N'-bis(benzyl)-substituted diamine and paullone derivatives are interesting candidates because of their potent and/or selective anti-trypanosomal and anti-trypanothione synthetase activity.

Targeting glycolysis in the malaria parasite Plasmodium falciparum

Glycolysis is the main pathway for ATP production in the malaria parasite Plasmodium falciparum and essential for its survival. Following a sensitivity analysis of a detailed kinetic model for glycolysis in the parasite, the glucose transport reaction was identified as the step whose activity needed to be inhibited to the least extent to result in a 50% reduction in glycolytic flux. In a subsequent inhibitor titration with cytochalasin B, we confirmed the model analysis experimentally and measured a flux control coefficient of 0.3 for the glucose transporter. In addition to the glucose transporter, the glucokinase and phosphofructokinase had high flux control coefficients, while for the ATPase a small negative flux control coefficient was predicted. In a broader comparative analysis of glycolytic models, we identified a weakness in the P. falciparum pathway design with respect to stability towards perturbations in the ATP demand.

DATABASE

The mathematical model described here has been submitted to the JWS Online Cellular Systems Modelling Database and can be accessed at http://jjj.bio.vu.nl/database/vanniekerk1. The SEEK-study including the experimental data set is available at DOI 10.15490/seek.1.

INVESTIGATION

56 (http://dx.doi.org/10.15490/seek.1.

INVESTIGATION

56).

Drug target identification through systems biology

To rationalise drug target selection, we should understand the role of putative targets in biological pathways quantitatively. We review here how experimental and computational network-based approaches can aid more rational drug target selection and illustrate this with results obtained for microbes and for cancer. Comparison of the drug response of biochemical networks in target cells and (healthy) host cells can reveal network-selective targets.

Number and regulation of protozoan aquaporins reflect environmental complexity

Protozoa are a diverse group of unicellular eukaryotes. Evidence has accumulated that protozoan aquaporin water and solute channels (AQP) contribute to adaptation in changing environments. Intracellular protozoan parasites live a well-sheltered life. Plasmodium spp. express a single AQP, Toxoplasma gondii two, while Trypanosoma cruzi and Leishamnia spp. encode up to five AQPs. Their AQPs are thought to import metabolic precursors and simultaneously to dispose of waste and to help parasites survive osmotic stress during transmission to and from the insect vector or during kidney passages. Trypanosoma brucei is a protozoan parasite that swims freely in the human blood. Expression and intracellular localization of the three T. brucei AQPs depend on the stage of differentiation during the life cycle, suggesting distinct roles in energy generation, metabolism, and cell motility. Free-living amoebae are in direct contact with the environment, encountering severe and sudden changes in the availability of nutrition, and in the osmotic conditions due to rainfall or drought. Amoeba proteus expresses a single AQP that is present in the contractile vacuole complex required for osmoregulation, whereas Dictyostelium discoideum expresses four AQPs, of which two are present in the single-celled amoeboidal stage and two more in the later multicellular stages preceding spore formation. The number and regulation of protozoan aquaporins may reflect environmental complexity. We highlight the gated AqpB from D. discoideum as an example of how life in the wild is challenged by a complex AQP structure-function relationship.

Construction and validation of a detailed kinetic model of glycolysis in Plasmodium falciparum

The enzymes in the Embden-Meyerhof-Parnas pathway of Plasmodium falciparum trophozoites were kinetically characterized and their integrated activities analyzed in a mathematical model. For validation of the model, we compared model predictions for steady-state fluxes and metabolite concentrations of the hexose phosphates with experimental values for intact parasites. The model, which is completely based on kinetic parameters that were measured for the individual enzymes, gives an accurate prediction of the steady-state fluxes and intermediate concentrations. This is the first detailed kinetic model for glucose metabolism in P. falciparum, one of the most prolific malaria-causing protozoa, and the high predictive power of the model makes it a strong tool for future drug target identification studies. The modelling workflow is transparent and reproducible, and completely documented in the SEEK platform, where all experimental data and model files are available for download.

DATABASE

The mathematical models described in the present study have been submitted to the JWS Online Cellular Systems Modelling Database (http://jjj.bio.vu.nl/database/penkler). The investigation and complete experimental data set is available on SEEK (10.15490/seek.1.

INVESTIGATION

56).

Probing the metabolic network in bloodstream-form Trypanosoma brucei using untargeted metabolomics with stable isotope labelled glucose

Metabolomics coupled with heavy-atom isotope-labelled glucose has been used to probe the metabolic pathways active in cultured bloodstream form trypomastigotes of Trypanosoma brucei, a parasite responsible for human African trypanosomiasis. Glucose enters many branches of metabolism beyond glycolysis, which has been widely held to be the sole route of glucose metabolism. Whilst pyruvate is the major end-product of glucose catabolism, its transamination product, alanine, is also produced in significant quantities. The oxidative branch of the pentose phosphate pathway is operative, although the non-oxidative branch is not. Ribose 5-phosphate generated through this pathway distributes widely into nucleotide synthesis and other branches of metabolism. Acetate, derived from glucose, is found associated with a range of acetylated amino acids and, to a lesser extent, fatty acids; while labelled glycerol is found in many glycerophospholipids. Glucose also enters inositol and several sugar nucleotides that serve as precursors to macromolecule biosynthesis. Although a Krebs cycle is not operative, malate, fumarate and succinate, primarily labelled in three carbons, were present, indicating an origin from phosphoenolpyruvate via oxaloacetate. Interestingly, the enzyme responsible for conversion of phosphoenolpyruvate to oxaloacetate, phosphoenolpyruvate carboxykinase, was shown to be essential to the bloodstream form trypanosomes, as demonstrated by the lethal phenotype induced by RNAi-mediated downregulation of its expression. In addition, glucose derivatives enter pyrimidine biosynthesis via oxaloacetate as a precursor to aspartate and orotate.

In this work we have followed the distribution of carbon derived from glucose in bloodstream form trypanosomes, the causative agent of African trypanosomiasis, revealing it to enter a diverse range of metabolites. The work involved using ¹³C-labelled glucose and following the fate of the labelled carbon with an LC-MS based metabolomics platform. Beyond glycolysis and the oxidative branch of the pentose phosphate pathway the label entered lipid biosynthesis both through glycerol 3-phosphate and also acetate. Glucose derived carbon also entered nucleotide synthesis through ribose and pyrimidine synthesis through oxaloacetate-derived aspartate. Appreciable quantities of the carboxylic acids succinate and malate were identified, although labelling patterns indicate they are not TCA cycle derived. Amino sugars and sugar nucleotides were also labelled as was inositol used in protein modification but not in inositol phospholipid headgroup production. We confirm active and essential oxaloacetate production in bloodstream form trypanosomes and show that phosphoenolpyruvate carboxykinase is essential to these parasites using RNA interference. The amount of glucose entering these metabolites is minor compared to the quantity that enters pyruvate excreted from the cell, but the observation that enzymes contributing to the metabolism of glucose beyond glycolysis can be essential offers potential new targets for chemotherapy against trypanosomiasis.

25 years of African trypanosome research: From description to molecular dissection and new drug discovery

The Molecular Parasitology conference was first held at the Marine Biological laboratory, Woods Hole, USA 25 years ago. Since that first meeting, the conference has evolved and expanded but has remained the showcase for the latest research developments in molecular parasitology. In this perspective, I reflect on the scientific discoveries focussed on African trypanosomes (Trypanosoma brucei spp.) that have occurred since the inaugural MPM meeting and discuss the current and future status of research on these parasites.

LC-MS-based absolute metabolite quantification: application to metabolic flux measurement in trypanosomes

Human African trypanosomiasis is a neglected tropical disease caused by the protozoan parasite, Trypanosoma brucei. In the mammalian bloodstream, the trypanosome's metabolism differs significantly from that of its host. For example, the parasite relies exclusively on glycolysis for energy source. Recently, computational and mathematical models of trypanosome metabolism have been generated to assist in understanding the parasite metabolism with the aim of facilitating drug development. Optimisation of these models requires quantitative information, including metabolite concentrations and/or metabolic fluxes that have been hitherto unavailable on a large scale. Here, we have implemented an LC-MS-based method that allows large scale quantification of metabolite levels by using U-¹³C-labelled E.coli extracts as internal standards. Known amounts of labelled E. coli extract were added into the parasite samples, as well as calibration standards, and used to obtain calibration curves enabling us to convert intensities into concentrations. This method allowed us to reliably quantify the changes of 43 intracellular metabolites and 32 extracellular metabolites in the medium over time. Based on the absolute quantification, we were able to compute consumption and production fluxes. These quantitative data can now be used to optimise computational models of parasite metabolism.

The silicon trypanosome: a test case of iterative model extension in systems biology

The African trypanosome, Trypanosoma brucei, is a unicellular parasite causing African Trypanosomiasis (sleeping sickness in humans and nagana in animals). Due to some of its unique properties, it has emerged as a popular model organism in systems biology. A predictive quantitative model of glycolysis in the bloodstream form of the parasite has been constructed and updated several times. The Silicon Trypanosome is a project that brings together modellers and experimentalists to improve and extend this core model with new pathways and additional levels of regulation. These new extensions and analyses use computational methods that explicitly take different levels of uncertainty into account. During this project, numerous tools and techniques have been developed for this purpose, which can now be used for a wide range of different studies in systems biology.

Comparative ribosome profiling reveals extensive translational complexity in different Trypanosoma brucei life cycle stages

While gene expression is a fundamental and tightly controlled cellular process that is regulated at multiple steps, the exact contribution of each step remains unknown in any organism. The absence of transcription initiation regulation for RNA polymerase II in the protozoan parasite Trypanosoma brucei greatly simplifies the task of elucidating the contribution of translation to global gene expression. Therefore, we have sequenced ribosome-protected mRNA fragments in T. brucei, permitting the genome-wide analysis of RNA translation and translational efficiency. We find that the latter varies greatly between life cycle stages of the parasite and ∼100-fold between genes, thus contributing to gene expression to a similar extent as RNA stability. The ability to map ribosome positions at sub-codon resolution revealed extensive translation from upstream open reading frames located within 5' UTRs and enabled the identification of hundreds of previously un-annotated putative coding sequences (CDSs). Evaluation of existing proteomics and genome-wide RNAi data confirmed the translation of previously un-annotated CDSs and suggested an important role for >200 of those CDSs in parasite survival, especially in the form that is infective to mammals. Overall our data show that translational control plays a prevalent and important role in different parasite life cycle stages of T. brucei.

Handling uncertainty in dynamic models: the pentose phosphate pathway in Trypanosoma brucei

Dynamic models of metabolism can be useful in identifying potential drug targets, especially in unicellular organisms. A model of glycolysis in the causative agent of human African trypanosomiasis, Trypanosoma brucei, has already shown the utility of this approach. Here we add the pentose phosphate pathway (PPP) of T. brucei to the glycolytic model. The PPP is localized to both the cytosol and the glycosome and adding it to the glycolytic model without further adjustments leads to a draining of the essential bound-phosphate moiety within the glycosome. This phosphate “leak” must be resolved for the model to be a reasonable representation of parasite physiology. Two main types of theoretical solution to the problem could be identified: (i) including additional enzymatic reactions in the glycosome, or (ii) adding a mechanism to transfer bound phosphates between cytosol and glycosome. One example of the first type of solution would be the presence of a glycosomal ribokinase to regenerate ATP from ribose 5-phosphate and ADP. Experimental characterization of ribokinase in T. brucei showed that very low enzyme levels are sufficient for parasite survival, indicating that other mechanisms are required in controlling the phosphate leak. Examples of the second type would involve the presence of an ATP:ADP exchanger or recently described permeability pores in the glycosomal membrane, although the current absence of identified genes encoding such molecules impedes experimental testing by genetic manipulation. Confronted with this uncertainty, we present a modeling strategy that identifies robust predictions in the context of incomplete system characterization. We illustrate this strategy by exploring the mechanism underlying the essential function of one of the PPP enzymes, and validate it by confirming the model predictions experimentally.

Mathematical models have been valuable tools for investigating the complex behaviors of metabolism. Due to incomplete knowledge of biological systems, these models contain inevitable uncertainty. This uncertainty is present in the measured or estimated parameter values, but also in the structure of the metabolic network. In this paper we increase the coverage of a particularly well studied model of glucose metabolism in Trypanosoma brucei, a tropical parasite that causes African sleeping sickness, by extending it with an additional pathway in two compartments. During this modeling process we highlighted uncertainties in parameter values and network structure and used these to formulate new hypotheses which were subsequently tested experimentally. The models were improved with the experimentally derived data, but uncertainty remained concerning the exact topology of the system. These models allowed us to investigate the effects of the loss of one enzyme, 6-phosphogluconate dehydrogenase. By taking uncertainty into account, the models demonstrated that the loss of this enzyme is lethal to the parasite by a mechanism different than that in other organisms. Our methodology shows how formally introducing uncertainty into model building provides robust model behavior that is independent of the network structure or parameter values.

New types of experimental data shape the use of enzyme kinetics for dynamic network modeling

Since the publication of Leonor Michaelis and Maude Menten's paper on the reaction kinetics of the enzyme invertase in 1913, molecular biology has evolved tremendously. New measurement techniques allow in vivo characterization of the whole genome, proteome or transcriptome of cells, whereas the classical enzyme essay only allows determination of the two Michaelis-Menten parameters V and K(m). Nevertheless, Michaelis-Menten kinetics are still commonly used, not only in the in vitro context of enzyme characterization but also as a rate law for enzymatic reactions in larger biochemical reaction networks. In this review, we give an overview of the historical development of kinetic rate laws originating from Michaelis-Menten kinetics over the past 100 years. Furthermore, we briefly summarize the experimental techniques used for the characterization of enzymes, and discuss web resources that systematically store kinetic parameters and related information. Finally, describe the novel opportunities that arise from using these data in dynamic mathematical modeling. In this framework, traditional in vitro approaches may be combined with modern genome-scale measurements to foster thorough understanding of the underlying complex mechanisms.

Explicit consideration of topological and parameter uncertainty gives new insights into a well-established model of glycolysis

Previous models of glycolysis in the sleeping sickness parasite Trypanosoma brucei assumed that the core part of glycolysis in this unicellular parasite is tightly compartimentalized within an organelle, the glycosome, which had previously been shown to contain most of the glycolytic enzymes. The glycosomes were assumed to be largely impermeable, and exchange of metabolites between the cytosol and the glycosome was assumed to be regulated by specific transporters in the glycosomal membrane. This tight compartmentalization was considered to be essential for parasite viability. Recently, size-specific metabolite pores were discovered in the membrane of glycosomes. These channels are proposed to allow smaller metabolites to diffuse across the membrane but not larger ones. In light of this new finding, we re-analyzed the model taking into account uncertainty about the topology of the metabolic system in T. brucei, as well as uncertainty about the values of all parameters of individual enzymatic reactions. Our analysis shows that these newly-discovered nonspecific pores are not necessarily incompatible with our current knowledge of the glycosomal metabolic system, provided that the known cytosolic activities of the glycosomal enzymes play an important role in the regulation of glycolytic fluxes and the concentration of metabolic intermediates of the pathway.

Proteins and lipids of glycosomal membranes from Leishmania tarentolae and Trypanosoma brucei

In kinetoplastid protists, several metabolic pathways, including glycolysis and purine salvage, are located in glycosomes, which are microbodies that are evolutionarily related to peroxisomes. With the exception of some potential transporters for fatty acids, and one member of the mitochondrial carrier protein family, proteins that transport metabolites across the glycosomal membrane have yet to be identified. We show here that the phosphatidylcholine species composition of Trypanosoma brucei glycosomal membranes resembles that of other cellular membranes, which means that glycosomal membranes are expected to be impermeable to small hydrophilic molecules unless transport is facilitated by specialized membrane proteins. Further, we identified 464 proteins in a glycosomal membrane preparation from Leishmania tarentolae. The proteins included approximately 40 glycosomal matrix proteins, and homologues of peroxisomal membrane proteins - PEX11, GIM5A and GIM5B; PXMP4, PEX2 and PEX16 - as well as the transporters GAT1 and GAT3. There were 27 other proteins that could not be unambiguously assigned to other compartments, and that had predicted trans-membrane domains. However, no clear candidates for transport of the major substrates and intermediates of energy metabolism were found. We suggest that, instead, these metabolites are transported via pores formed by the known glycosomal membrane proteins.

Trade-off of dynamic fragility but not of robustness in metabolic pathways in silico

Selective robustness is a key feature of biochemical networks. It confers a fitness benefit to organisms living in dynamic environments. The (in-)sensitivity of a network to external perturbations results from the interplay between network dynamics, structure and enzyme kinetics. In this work, we focus on the subtle interplay between robustness and control (fragility). We describe a quantitative method for defining the fragility and robustness of system fluxes to perturbations. We find that for many mathematical models of metabolic pathways, the robustness of fluxes vis-à-vis perturbations of all the enzyme activities is captured by a broad distribution of the robustness coefficients. We find that in cases where a metabolic pathway flux is made less robust with respect to the perturbation of a particular network step, the average robustness may still be increased. We then show that fragility is conserved upon a perturbation of network processes and equate fragility with control as defined in metabolic control analysis. This highlights the non-intuitive nature of the interplay between fragility and robustness and the need for a dynamic network understanding.

A new family of covalent inhibitors block nucleotide binding to the active site of pyruvate kinase

PYK (pyruvate kinase) plays a central role in the metabolism of many organisms and cell types, but the elucidation of the details of its function in a systems biology context has been hampered by the lack of specific high-affinity small-molecule inhibitors. High-throughput screening has been used to identify a family of saccharin derivatives which inhibit LmPYK (Leishmania mexicana PYK) activity in a time- (and dose-) dependent manner, a characteristic of irreversible inhibition. The crystal structure of DBS {4-[(1,1-dioxo-1,2-benzothiazol-3-yl)sulfanyl]benzoic acid} complexed with LmPYK shows that the saccharin moiety reacts with an active-site lysine residue (Lys335), forming a covalent bond and sterically hindering the binding of ADP/ATP. Mutation of the lysine residue to an arginine residue eliminated the effect of the inhibitor molecule, providing confirmation of the proposed inhibitor mechanism. This lysine residue is conserved in the active sites of the four human PYK isoenzymes, which were also found to be irreversibly inhibited by DBS. X-ray structures of PYK isoforms show structural differences at the DBS-binding pocket, and this covalent inhibitor of PYK provides a chemical scaffold for the design of new families of potentially isoform-specific irreversible inhibitors.

Determining enzyme kinetics for systems biology with nuclear magnetic resonance spectroscopy

Enzyme kinetics for systems biology should ideally yield information about the enzyme’s activity under in vivo conditions, including such reaction features as substrate cooperativity, reversibility and allostery, and be applicable to enzymatic reactions with multiple substrates. A large body of enzyme-kinetic data in the literature is based on the uni-substrate Michaelis-Menten equation, which makes unnatural assumptions about enzymatic reactions (e.g., irreversibility), and its application in systems biology models is therefore limited. To overcome this limitation, we have utilised NMR time-course data in a combined theoretical and experimental approach to parameterize the generic reversible Hill equation, which is capable of describing enzymatic reactions in terms of all the properties mentioned above and has fewer parameters than detailed mechanistic kinetic equations; these parameters are moreover defined operationally. Traditionally, enzyme kinetic data have been obtained from initial-rate studies, often using assays coupled to NAD(P)H-producing or NAD(P)H-consuming reactions. However, these assays are very labour-intensive, especially for detailed characterisation of multi-substrate reactions. We here present a cost-effective and relatively rapid method for obtaining enzyme-kinetic parameters from metabolite time-course data generated using NMR spectroscopy. The method requires fewer runs than traditional initial-rate studies and yields more information per experiment, as whole time-courses are analyzed and used for parameter fitting. Additionally, this approach allows real-time simultaneous quantification of all metabolites present in the assay system (including products and allosteric modifiers), which demonstrates the superiority of NMR over traditional spectrophotometric coupled enzyme assays. The methodology presented is applied to the elucidation of kinetic parameters for two coupled glycolytic enzymes from Escherichia coli (phosphoglucose isomerase and phosphofructokinase). ³¹P-NMR time-course data were collected by incubating cell extracts with substrates, products and modifiers at different initial concentrations. NMR kinetic data were subsequently processed using a custom software module written in the Python programming language, and globally fitted to appropriately modified Hill equations.

Flux Analysis of the Trypanosoma brucei Glycolysis Based on a Multiobjective-Criteria Bioinformatic Approach

Trypanosoma brucei is a protozoan parasite of major of interest in discovering new genes for drug targets. This parasite alternates its life cycle between the mammal host(s) (bloodstream form) and the insect vector (procyclic form), with two divergent glucose metabolism amenable to in vitro culture. While the metabolic network of the bloodstream forms has been well characterized, the flux distribution between the different branches of the glucose metabolic network in the procyclic form has not been addressed so far. We present a computational analysis (called Metaboflux) that exploits the metabolic topology of the procyclic form, and allows the incorporation of multipurpose experimental data to increase the biological relevance of the model. The alternatives resulting from the structural complexity of networks are formulated as an optimization problem solved by a metaheuristic where experimental data are modeled in a multiobjective function. Our results show that the current metabolic model is in agreement with experimental data and confirms the observed high metabolic flexibility of glucose metabolism. In addition, Metaboflux offers a rational explanation for the high flexibility in the ratio between final products from glucose metabolism, thsat is, flux redistribution through the malic enzyme steps.

Mass conserved elementary kinetics is sufficient for the existence of a non-equilibrium steady state concentration

Living systems are forced away from thermodynamic equilibrium by exchange of mass and energy with their environment. In order to model a biochemical reaction network in a non-equilibrium state one requires a mathematical formulation to mimic this forcing. We provide a general formulation to force an arbitrary large kinetic model in a manner that is still consistent with the existence of a non-equilibrium steady state. We can guarantee the existence of a non-equilibrium steady state assuming only two conditions; that every reaction is mass balanced and that continuous kinetic reaction rate laws never lead to a negative molecule concentration. These conditions can be verified in polynomial time and are flexible enough to permit one to force a system away from equilibrium. With expository biochemical examples we show how reversible, mass balanced perpetual reaction(s), with thermodynamically infeasible kinetic parameters, can be used to perpetually force various kinetic models in a manner consistent with the existence of a steady state. Easily testable existence conditions are foundational for efforts to reliably compute non-equilibrium steady states in genome-scale biochemical kinetic models.

Kinetic modelling of GlmU reactions - prioritization of reaction for therapeutic application

Mycobacterium tuberculosis(Mtu), a successful pathogen, has developed resistance against the existing anti-tubercular drugs necessitating discovery of drugs with novel action. Enzymes involved in peptidoglycan biosynthesis are attractive targets for antibacterial drug discovery. The bifunctional enzyme mycobacterial GlmU (Glucosamine 1-phosphate N-acetyltransferase/ N-acetylglucosamine-1-phosphate uridyltransferase) has been a target enzyme for drug discovery. Its C- and N- terminal domains catalyze acetyltransferase (rxn-1) and uridyltransferase (rxn-2) activities respectively and the final product is involved in peptidoglycan synthesis. However, the bifunctional nature of GlmU poses difficulty in deciding which function to be intervened for therapeutic advantage. Genetic analysis showed this as an essential gene but it is still unclear whether any one or both of the activities are critical for cell survival. Often enzymatic activity with suitable high-throughput assay is chosen for random screening, which may not be the appropriate biological function inhibited for maximal effect. Prediction of rate-limiting function by dynamic network analysis of reactions could be an option to identify the appropriate function. With a view to provide insights into biochemical assays with appropriate activity for inhibitor screening, kinetic modelling studies on GlmU were undertaken. Kinetic model of Mtu GlmU-catalyzed reactions was built based on the available kinetic data on Mtu and deduction from Escherichia coli data. Several model variants were constructed including coupled/decoupled, varying metabolite concentrations and presence/absence of product inhibitions. This study demonstrates that in coupled model at low metabolite concentrations, inhibition of either of the GlmU reactions cause significant decrement in the overall GlmU rate. However at higher metabolite concentrations, rxn-2 showed higher decrement. Moreover, with available intracellular concentration of the metabolites and in vivo variant of model, uncompetitive inhibition of rxn-2 caused highest decrement. Thus, at physiologically relevant metabolite concentrations, targeting uridyltranferase activity of Mtu GlmU would be a better choice for therapeutic intervention.

Translocation of solutes and proteins across the glycosomal membrane of trypanosomes; possibilities and limitations for targeting with trypanocidal drugs

Glycosomes are specialized peroxisomes found in all kinetoplastid organisms. The organelles are unique in harbouring most enzymes of the glycolytic pathway. Matrix proteins, synthesized in the cytosol, cofactors and metabolites have to be transported across the membrane. Recent research on Trypanosoma brucei has provided insight into how these translocations across the membrane occur, although many details remain to be elucidated. Proteins are imported by a cascade of reactions performed by specialized proteins, called peroxins, in which a cytosolic receptor with bound matrix protein inserts itself in the membrane to deliver its cargo into the organelle and is subsequently retrieved from the glycosome to perform further rounds of import. Bulky solutes, such as cofactors and acyl-CoAs, seem to be translocated by specific transporter molecules, whereas smaller solutes such as glycolytic intermediates probably cross the membrane through pore-forming channels. The presence of such channels is in apparent contradiction with previous results that suggested a low permeability of the glycosomal membrane. We propose 3 possible, not mutually exclusive, solutions for this paradox. Glycosomal glycolytic enzymes have been validated as drug targets against trypanosomatid-borne diseases. We discuss the possible implications of the new data for the design of drugs to be delivered into glycosomes.

Understanding complexity in neurodegenerative diseases: in silico reconstruction of emergence

Healthy functioning is an emergent property of the network of interacting biomolecules that comprise an organism. It follows that disease (a network shift that causes malfunction) is also an emergent property, emerging from a perturbation of the network. On the one hand, the biomolecular network of every individual is unique and this is evident when similar disease-producing agents cause different individual pathologies. Consequently, a personalized model and approach for every patient may be required for therapies to become effective across mankind. On the other hand, diverse combinations of internal and external perturbation factors may cause a similar shift in network functioning. We offer this as an explanation for the multi-factorial nature of most diseases: they are "systems biology diseases," or "network diseases." Here we use neurodegenerative diseases, like Parkinson's disease (PD), as an example to show that due to the inherent complexity of these networks, it is difficult to understand multi-factorial diseases with simply our "naked brain." When describing interactions between biomolecules through mathematical equations and integrating those equations into a mathematical model, we try to reconstruct the emergent properties of the system in silico. The reconstruction of emergence from interactions between huge numbers of macromolecules is one of the aims of systems biology. Systems biology approaches enable us to break through the limitation of the human brain to perceive the extraordinarily large number of interactions, but this also means that we delegate the understanding of reality to the computer. We no longer recognize all those essences in the system's design crucial for important physiological behavior (the so-called "design principles" of the system). In this paper we review evidence that by using more abstract approaches and by experimenting in silico, one may still be able to discover and understand the design principles that govern behavioral emergence.