Genotype networks,innovation, and robustness in sulfur metabolism

Background  

A metabolism is a complex network of chemical reactions. This network synthesizes multiple small precursor molecules of biomass from chemicals that occur in the environment. The metabolic network of any one organism is encoded by a metabolic genotype, defined as the set of enzyme-coding genes whose products catalyze the network's reactions. Each metabolic genotype has a metabolic phenotype. We define this metabolic phenotype as the spectrum of different sources of a chemical element that a metabolism can use to synthesize biomass. We here focus on the element sulfur. We study properties of the space of all possible metabolic genotypes in sulfur metabolism by analyzing random metabolic genotypes that are viable on different numbers of sulfur sources.     

Trade-Offs Predicted by Metabolic Network Structure Give Rise to Evolutionary Specialization and Phenotypic Diversification

Mitigating trade-offs between different resource-utilization functions is key to an organism’s ecological and evolutionary success. These trade-offs often reflect metabolic constraints with a complex molecular underpinning; therefore, their consequences for evolutionary processes have remained elusive. Here, we investigate how metabolic architecture induces resource-utilization constraints and how these constraints, in turn, elicit evolutionary specialization and diversification. Guided by the metabolic network structure of the bacterium Lactococcus cremoris, we selected two carbon sources (fructose and galactose) with predicted coutilization constraints. By evolving L. cremoris on either fructose, galactose, or a mix of both sugars, we imposed selection favoring divergent metabolic specializations or coutilization of both resources, respectively. Phenotypic characterization revealed the evolution of either fructose or galactose specialists in the single-sugar treatments. In the mixed-sugar regime, we observed adaptive diversification: both specialists coexisted, and no generalist evolved. Divergence from the ancestral phenotype occurred at key pathway junctions in the central carbon metabolism. Fructose specialists evolved mutations in the fbp and pfk genes that appear to balance anabolic and catabolic carbon fluxes. Galactose specialists evolved increased expression of pgmA (the primary metabolic bottleneck of galactose metabolism) and silencing of ptnABCD (the main glucose transporter) and ldh (regulator/enzyme of downstream carbon metabolism). Overall, our study shows how metabolic network architecture and historical contingency serve to predict targets of selection and inform the functional interpretation of evolved mutations. The elucidation of the relationship between molecular constraints and phenotypic trade-offs contributes to an integrative understanding of evolutionary specialization and diversification.     

Biomolecular hydrogels formed and degraded via site-specific enzymatic reactions

We present polymeric hydrogel biomaterials that are biomimetic both in their synthesis and degradation. The design of oligopeptide building blocks with dual enzymatic responsiveness allows us to create polymer networks that are formed and functionalized via enzymatic reactions and are degradable via other enzymatic reactions, both occurring under physiological conditions. The activated transglutaminase enzyme factor XIIIa was utilized for site-specific coupling of prototypical cell adhesion ligands and for simultaneous cross-linking of hydrogel networks from factor XIIIa substrate-modified multiarm poly(ethylene glycol) macromers. Ligand incorporation is nearly quantitative and thus controllable, and does not alter the network's macroscopic properties over a concentration range that elicits specific cell adhesion. Living mammalian cells can be encapsulated in the gels without any noticeable decrease in viability. The degradation of gels can be engineered to occur, for example, via cell-secreted matrix metalloproteinases, thus rendering these gels interesting for biomedical applications such as drug delivery systems or smart implants for in situ tissue engineering.     

Effect of regulatory mechanism on hyperbolic reaction network properties

The notion that the regulated and flux-controlling enzyme in a metabolic network need not correspond suggests that the purpose of regulation may not be flux homeostasis under all physiological circumstances. Additionally, the fact that diversity in the function of intact metabolic networks exists suggests that in addition to time constant separation, other kinetic structure/regulatory mechanism patterns exist. In order to compliment and expand prior work on identifying kinetic structure-property relationships in networks, the present work explores in a general way how the control, dynamic, and energetic properties of metabolic networks depend on operating point, kinetic structure, and regulatory mechanism. The basic feature of trade-offs between properties is illustrated and used as a basis for indicating how particular subsets of structure, regulatory mechanism, and operating point emphasize certain properties that can be associated with a physiological function. Examples of scavenging trace metabolites and amphibolite coordination are proposed. Microstructure logic in terms of turnover number distributions as well as a potential mixed polynomial network analysis approach are also discussed.     

Attractor Metabolic Networks

Background

The experimental observations and numerical studies with dissipative metabolic networks have shown that cellular enzymatic activity self-organizes spontaneously leading to the emergence of a Systemic Metabolic Structure in the cell, characterized by a set of different enzymatic reactions always locked into active states (metabolic core) while the rest of the catalytic processes are only intermittently active. This global metabolic structure was verified for Escherichia coli, Helicobacter pylori and Saccharomyces cerevisiae, and it seems to be a common key feature to all cellular organisms. In concordance with these observations, the cell can be considered a complex metabolic network which mainly integrates a large ensemble of self-organized multienzymatic complexes interconnected by substrate fluxes and regulatory signals, where multiple autonomous oscillatory and quasi-stationary catalytic patterns simultaneously emerge. The network adjusts the internal metabolic activities to the external change by means of flux plasticity and structural plasticity.

Methodology/Principal Findings

In order to research the systemic mechanisms involved in the regulation of the cellular enzymatic activity we have studied different catalytic activities of a dissipative metabolic network under different external stimuli. The emergent biochemical data have been analysed using statistical mechanic tools, studying some macroscopic properties such as the global information and the energy of the system. We have also obtained an equivalent Hopfield network using a Boltzmann machine. Our main result shows that the dissipative metabolic network can behave as an attractor metabolic network.

Conclusions/Significance

We have found that the systemic enzymatic activities are governed by attractors with capacity to store functional metabolic patterns which can be correctly recovered from specific input stimuli. The network attractors regulate the catalytic patterns, modify the efficiency in the connection between the multienzymatic complexes, and stably retain these modifications. Here for the first time, we have introduced the general concept of attractor metabolic network, in which this dynamic behavior is observed.     

Beyond topology: coevolution of structure and flux in metabolic networks

Interactions between the structure of a metabolic network and its functional properties underlie its evolutionary diversification, but the mechanism by which such interactions arise remains elusive. Particularly unclear is whether metabolic fluxes that determine the concentrations of compounds produced by a metabolic network, are causally linked to a network's structure or emerge independently of it. A direct empirical study of populations where both structural and functional properties vary among individuals’ metabolic networks is required to establish whether changes in structure affect the distribution of metabolic flux. In a population of house finches (Haemorhous mexicanus), we reconstructed full carotenoid metabolic networks for 442 individuals and uncovered 11 structural variants of this network with different compounds and reactions. We examined the consequences of this structural diversity for the concentrations of plumage‐bound carotenoids produced by flux in these networks. We found that concentrations of metabolically derived, but not dietary carotenoids, depended on network structure. Flux was partitioned similarly among compounds in individuals of the same network structure: within each network, compound concentrations were closely correlated. The highest among‐individual variation in flux occurred in networks with the strongest among‐compound correlations, suggesting that changes in the magnitude, but not the distribution of flux, underlie individual differences in compound concentrations on a static network structure. These findings indicate that the distribution of flux in carotenoid metabolism closely follows network structure. Thus, evolutionary diversification and local adaptations in carotenoid metabolism may depend more on the gain or loss of enzymatic reactions than on changes in flux within a network structure.     

Candidate metabolic network states in human mitochondria. Impact of diabetes, ischemia, and diet

The human mitochondrial metabolic network was recently reconstructed based on proteomic and biochemical data. Linear programming and uniform random sampling were applied herein to identify candidate steady states of the metabolic network that were consistent with the imposed physico-chemical constraints and available experimental data. The activity of the mitochondrion was studied under four metabolic conditions: normal physiologic, diabetic, ischemic, and dietetic. Pairwise correlations between steady-state reaction fluxes were calculated in each condition to evaluate the dependence among the reactions in the network. Applying constraints on exchange fluxes resulted in predictions for intracellular fluxes that agreed with experimental data. Analyses of the steady-state flux distributions showed that the experimentally observed reduced activity of pyruvate dehydrogenase in vivo could be a result of stoichiometric constraints and therefore would not necessarily require enzymatic inhibition. The observed changes in the energy metabolism of the mitochondrion under diabetic conditions were used to evaluate the impact of previously suggested treatments. The results showed that neither normalized glucose uptake nor decreased ketone body uptake have a positive effect on the mitochondrial energy metabolism or network flexibility. Taken together, this study showed that sampling of the steady-state flux space is a powerful method to investigate network properties under different conditions and provides a basis for in silico evaluations of effects of potential disease treatments.     

Functional and topological characterization of protein interaction networks

The elucidation of the cell's large-scale organization is a primary challenge for post-genomic biology, and understanding the structure of protein interaction networks offers an important starting point for such studies. We compare four available databases that approximate the protein interaction network of the yeast, Saccharomyces cerevisiae, aiming to uncover the network's generic large-scale properties and the impact of the proteins' function and cellular localization on the network topology. We show how each database supports a scale-free, topology with hierarchical modularity, indicating that these features represent a robust and generic property of the protein interactions network. We also find strong correlations between the network's structure and the functional role and subcellular localization of its protein constituents, concluding that most functional and/or localization classes appear as relatively segregated subnetworks of the full protein interaction network. The uncovered systematic differences between the four protein interaction databases reflect their relative coverage for different functional and localization classes and provide a guide for their utility in various bioinformatics studies.     

Properties of metabolic graphs: biological organization or representation artifacts?

Background  

Standard graphs, where each edge links two nodes, have been extensively used to represent the connectivity of metabolic networks. It is based on this representation that properties of metabolic networks, such as hierarchical and small-world structures, have been elucidated and null models have been proposed to derive biological organization hypotheses. However, these graphs provide a simplistic model of a metabolic network's connectivity map, since metabolic reactions often involve more than two reactants. In other words, this map is better represented as a hypergraph. Consequently, a question that naturally arises in this context is whether these properties truly reflect biological organization or are merely an artifact of the representation.     

Monte-Carlo Modeling of the Central Carbon Metabolism of Lactococcus lactis: Insights into Metabolic Regulation

Metabolic pathways are complex dynamic systems whose response to perturbations and environmental challenges are governed by multiple interdependencies between enzyme properties, reactions rates, and substrate levels. Understanding the dynamics arising from such a network can be greatly enhanced by the construction of a computational model that embodies the properties of the respective system. Such models aim to incorporate mechanistic details of cellular interactions to mimic the temporal behavior of the biochemical reaction system and usually require substantial knowledge of kinetic parameters to allow meaningful conclusions. Several approaches have been suggested to overcome the severe data requirements of kinetic modeling, including the use of approximative kinetics and Monte-Carlo sampling of reaction parameters. In this work, we employ a probabilistic approach to study the response of a complex metabolic system, the central metabolism of the lactic acid bacterium Lactococcus lactis, subject to perturbations and brief periods of starvation. Supplementing existing methodologies, we show that it is possible to acquire a detailed understanding of the control properties of a corresponding metabolic pathway model that is directly based on experimental observations. In particular, we delineate the role of enzymatic regulation to maintain metabolic stability and metabolic recovery after periods of starvation. It is shown that the feedforward activation of the pyruvate kinase by fructose-1,6-bisphosphate qualitatively alters the bifurcation structure of the corresponding pathway model, indicating a crucial role of enzymatic regulation to prevent metabolic collapse for low external concentrations of glucose. We argue that similar probabilistic methodologies will help our understanding of dynamic properties of small-, medium- and large-scale metabolic networks models.     

A probabilistic framework for the exploration of enzymatic capabilities based on feasible kinetics and control analysis

BackgroundAnalysis of limiting steps within enzyme-catalyzed reactions is fundamental to understand their behavior and regulation. Methods capable of unravelling control properties and exploring kinetic capabilities of enzymatic reactions would be particularly useful for protein and metabolic engineering. While single-enzyme control analysis formalism has previously been applied to well-studied enzymatic mechanisms, broader application of this formalism is limited in practice by the limited amount of kinetic data and the difficulty of describing complex allosteric mechanisms.MethodsTo overcome these limitations, we present here a probabilistic framework enabling control analysis of previously unexplored mechanisms under uncertainty. By combining a thermodynamically consistent parameterization with an efficient Sequential Monte Carlo sampler embedded in a Bayesian setting, this framework yields insights into the capabilities of enzyme-catalyzed reactions with modest kinetic information, provided that the catalytic mechanism and a thermodynamic reference point are defined.ResultsThe framework was used to unravel the impact of thermodynamic affinity, substrate saturation levels and effector concentrations on the flux control and response coefficients of a diverse set of enzymatic reactions.ConclusionsOur results highlight the importance of the metabolic context in the control analysis of isolated enzymes as well as the use of statistically sound methods for their interpretation.General SignificanceThis framework significantly expands our current capabilities for unravelling the control properties of general reaction kinetics with limited amount of information. This framework will be useful for both theoreticians and experimentalists in the field.     

Pathway Length and Evolutionary Constraint in Amino Acid Biosynthesis

The evolutionary properties of a metabolic network may be determined by the topology of the network. One attribute of pathways that make up the network is the number of enzymatic steps between initial substrates and final products. To determine the effect of pathway length on evolutionary lability of pathway structure, we examined amino acid biosynthetic pathways across 48 sequenced organisms. We demonstrate that longer pathways exhibit lower rates of change in pathway structure than shorter pathways. This finding suggests that increasing complexity may increase constraint on evolutionary change. (Matthew T. Rutter and Rebecca A. Zufall) Both authors contributed equally to this work.     

A quantitative theory of solid tumor growth, metabolic rate and vascularization

The relationships between cellular, structural and dynamical properties of tumors have traditionally been studied separately. Here, we construct a quantitative, predictive theory of solid tumor growth, metabolic rate, vascularization and necrosis that integrates the relationships between these properties. To accomplish this, we develop a comprehensive theory that describes the interface and integration of the tumor vascular network and resource supply with the cardiovascular system of the host. Our theory enables a quantitative understanding of how cells, tissues, and vascular networks act together across multiple scales by building on recent theoretical advances in modeling both healthy vasculature and the detailed processes of angiogenesis and tumor growth. The theory explicitly relates tumor vascularization and growth to metabolic rate, and yields extensive predictions for tumor properties, including growth rates, metabolic rates, degree of necrosis, blood flow rates and vessel sizes. Besides these quantitative predictions, we explain how growth rates depend on capillary density and metabolic rate, and why similar tumors grow slower and occur less frequently in larger animals, shedding light on Peto's paradox. Various implications for potential therapeutic strategies and further research are discussed.     

The thermodynamic meaning of metabolic exchange fluxes

Metabolic flux analysis (MFA) deals with the experimental determination of steady-state fluxes in metabolic networks. An important feature of the ¹³C MFA method is its capability to generate information on both directions of bidirectional reaction steps given by exchange fluxes. The biological interpretation of these exchange fluxes and their relation to thermodynamic properties of the respective reaction steps has never been systematically investigated. As a central result, it is shown here that for a general class of enzyme reaction mechanisms the quotients of net and exchange fluxes measured by ¹³C MFA are coupled to Gibbs energies of the reaction steps. To establish this relation the concept of apparent flux ratios of enzymatic isotope-labeling networks is introduced and some computing rules for these flux ratios are given. Application of these rules reveals a conceptional pitfall of ¹³C MFA, which is the inherent dependency of measured exchange fluxes on the chosen tracer atom. However, it is shown that this effect can be neglected for typical biochemical reaction steps under physiological conditions. In this situation, the central result can be formulated as a two-sided inequality relating fluxes, pool sizes, and standard Gibbs energies. This relation has far-reaching consequences for metabolic flux analysis, quantitative metabolomics, and network thermodynamics.     

Parental and Mating Effort: Is There Necessarily a Trade‐Off?

One of the common assumptions in the study of the evolution of parental care is that trade-offs exist between parental investment and other fitness-related traits. In general, this body of work follows the traditional definition that parental investment (in the current offspring) decreases that individual's ability to invest in future reproduction ( Trivers 1972 ). However, examination of the empirical evidence shows that assuming a trade-off between parental and mating effort is not always appropriate. This overemphasis on a trade-off between mating and parental effort has arisen in part because of an oversimplification of female reproductive strategies, a failure to consider interactions between the sexes, and a tendency to consider behaviours as unifunctional, thereby ignoring the more complex relationship between mating and parental effort in many species. Here, we first examine the empirical evidence for trade-offs between mating and parental effort in males and females to ask when trade-offs occur and what pattern they take. By highlighting a number of exemplar species, we then explore how the presence or absence of trade-offs relates to mate choice and sexual selection in both sexes. Finally, we highlight the importance of considering individual variation, which has been particularly overlooked in examinations of female investment, and how preferences in one sex may influence the existence and our interpretation of apparent trade-offs in the other sex.     

Functional Roles of Slow Enzyme Conformational Changes in Network Dynamics

Extensive studies from different fields reveal that many macromolecules, especially enzymes, show slow transitions among different conformations. This phenomenon is named such things as dynamic disorder, heterogeneity, hysteretic or mnemonic enzymes across these different fields, and has been directly demonstrated by single molecule enzymology and NMR studies recently. We analyzed enzyme slow conformational changes in the context of regulatory networks. A single enzymatic reaction with slow conformational changes can filter upstream network noises, and can either resonantly respond to the system stimulus at certain frequencies or respond adaptively for sustained input signals of the network fluctuations. It thus can serve as a basic functional motif with properties that are normally for larger intermolecular networks in the field of systems biology. We further analyzed examples including enzymes functioning against pH fluctuations, metabolic state change of Artemia embryos, and kinetic insulation of fluctuations in metabolic networks. The study also suggests that hysteretic enzymes may be building blocks of synthetic networks with various properties such as narrow-banded filtering. The work fills the missing gap between studies on enzyme biophysics and network level dynamics, and reveals that the coupling between the two is functionally important; it also suggests that the conformational dynamics of some enzymes may be evolutionally selected.