Ranking network mechanisms by how they fit diverse experiments and deciding on E. coli's ammonium transport and assimilation network

Biological ammonium transporters from the Amt/Mep/Rh superfamily: mechanism, energetics, and technical limitations

The exchange of ammonium across cellular membranes is a fundamental process in all domains of life and is facilitated by the ubiquitous Amt/Mep/Rh transporter superfamily. Remarkably, despite a high structural conservation in all domains of life, these proteins have gained various biological functions during evolution. It is tempting to hypothesise that the physiological functions gained by these proteins may be explained at least in part by differences in the energetics of their translocation mechanisms. Therefore, in this review, we will explore our current knowledge of energetics of the Amt/Mep/Rh family, discuss variations in observations between different organisms, and highlight some technical drawbacks which have hampered effects at mechanistic characterisation. Through the review, we aim to provide a comprehensive overview of current understanding of the mechanism of transport of this unique and extraordinary Amt/Mep/Rh superfamily of ammonium transporters.

On paradoxes between optimal growth, metabolic control analysis, and flux balance analysis

In Microbiology it is often assumed that growth rate is maximal. This may be taken to suggest that the dependence of the growth rate on every enzyme activity is at the top of an inverse-parabolic function, i.e. that all flux control coefficients should equal zero. This might seem to imply that the sum of these flux control coefficients equals zero. According to the summation law of Metabolic Control Analysis (MCA) the sum of flux control coefficients should equal 1 however. And in Flux Balance Analysis (FBA) catabolism is often limited by a hard bound, causing catabolism to fully control the fluxes, again in apparent contrast with a flux control coefficient of zero. Here we resolve these paradoxes (apparent contradictions) in an analysis that uses the 'Edinburgh pathway', the 'Amsterdam pathway', as well as a generic metabolic network providing the building blocks or Gibbs energy for microbial growth. We review and show that (i) optimization depends on so-called enzyme control coefficients rather than the 'catalytic control coefficients' of MCA's summation law, (ii) when optimization occurs at fixed total protein, the former differ from the latter to the extent that they may all become equal to zero in the optimum state, (iii) in more realistic scenarios of optimization where catalytically inert biomass is compensating or maintenance metabolism is taken into consideration, the optimum enzyme concentrations should not be expected to equal those that maximize the specific growth rate, (iv) optimization may be in terms of yield rather than specific growth rate, which resolves the paradox because the sum of catalytic control coefficients on yield equals 0, (v) FBA effectively maximizes growth yield, and for yield the summation law states 0 rather than 1, thereby removing the paradox, (vi) furthermore, FBA then comes more often to a 'hard optimum' defined by a maximum catabolic flux and a catabolic-enzyme control coefficient of 1. The trade-off between maintenance metabolism and growth is highlighted as worthy of further analysis.

'Social' versus 'asocial' cells-dynamic competition flux balance analysis

In multicellular organisms cells compete for resources or growth factors. If any one cell type wins, the co-existence of diverse cell types disappears. Existing dynamic Flux Balance Analysis (dFBA) does not accommodate changes in cell density caused by competition. Therefore we here develop 'dynamic competition Flux Balance Analysis' (dcFBA). With total biomass synthesis as objective, lower-growth-yield cells were outcompeted even when cells synthesized mutually required nutrients. Signal transduction between cells established co-existence, which suggests that such 'socialness' is required for multicellularity. Whilst mutants with increased specific growth rate did not outgrow the other cell types, loss of social characteristics did enable a mutant to outgrow the other cells. We discuss that 'asocialness' rather than enhanced growth rates, i.e., a reduced sensitivity to regulatory factors rather than enhanced growth rates, may characterize cancer cells and organisms causing ecological blooms. Therapies reinforcing cross-regulation may therefore be more effective than those targeting replication rates.

Construction, Characterization, and Application of an Ammonium Transporter (AmtB) Deletion Mutant of the Nitrogen-Fixing Bacterium Kosakonia radicincitans GXGL-4A in Cucumis sativus L. Seedlings

Nitrogen is an important factor affecting crop yield, but excessive use of chemical nitrogen fertilizer has caused decline in nitrogen utilization and soil and water pollution. Reducing the utilization of chemical nitrogen fertilizers by biological nitrogen fixation (BNF) is feasible for green production of crops. However, there are few reports on how to have more ammonium produced by nitrogen-fixing bacteria (NFB) flow outside the cell. In the present study, the amtB gene encoding an ammonium transporter (AmtB) in the genome of NFB strain Kosakonia radicincitans GXGL-4A was deleted and the △amtB mutant was characterized. The results showed that deletion of the amtB gene had no influence on the growth of bacterial cells. The extracellular ammonium nitrogen (NH₄⁺) content of the △amtB mutant under nitrogen-free culture conditions was significantly higher than that of the wild-type strain GXGL-4A (WT-GXGL-4A), suggesting disruption of NH₄⁺ transport. Meanwhile, the plant growth-promoting effect in cucumber seedlings was visualized after fertilization using cells of the △amtB mutant. NFB fertilization continuously increased the cucumber rhizosphere soil pH. The nitrate nitrogen (NO₃^-) content in soil in the △amtB treatment group was significantly higher than that in the WT-GXGL-4A treatment group in the short term but there was no difference in soil NH₄⁺ contents between groups. Soil enzymatic activities varied during a 45-day assessment period, indicating that △amtB fertilization influenced soil nitrogen cycling in the cucumber rhizosphere. The results will provide a solid foundation for developing the NFB GXGL-4A into an efficient biofertilizer agent.

Recommendations for performing measurements of apparent equilibrium constants of enzyme-catalyzed reactions and for reporting the results of these measurements

The measurement of values of apparent equilibrium constants K' for enzyme-catalyzed reactions involve a substantial number of critical details, neglect of which could lead to systematic errors. Here, interferences, impurities in the substances used, and failure to achieve equilibrium are matters of substantial consequence. Careful reporting of results is of great importance if the results are to have archival value. Thus, attention must be paid to the identification of the substances, specification of the reaction(s), the conditions of reaction, the definition of the equilibrium constant(s) and standard states, the use of standard nomenclature, symbols, and units, and uncertainties. This document contains a general discussion of various aspects of these equilibrium measurements as well as STRENDA (Standards for Reporting Enzymology Data) recommendations regarding the measurements and the reporting of results.

MLAGO: machine learning-aided global optimization for Michaelis constant estimation of kinetic modeling

BACKGROUND

Kinetic modeling is a powerful tool for understanding the dynamic behavior of biochemical systems. For kinetic modeling, determination of a number of kinetic parameters, such as the Michaelis constant (K_m), is necessary, and global optimization algorithms have long been used for parameter estimation. However, the conventional global optimization approach has three problems: (i) It is computationally demanding. (ii) It often yields unrealistic parameter values because it simply seeks a better model fitting to experimentally observed behaviors. (iii) It has difficulty in identifying a unique solution because multiple parameter sets can allow a kinetic model to fit experimental data equally well (the non-identifiability problem).

RESULTS

To solve these problems, we propose the Machine Learning-Aided Global Optimization (MLAGO) method for K_m estimation of kinetic modeling. First, we use a machine learning-based K_m predictor based only on three factors: EC number, KEGG Compound ID, and Organism ID, then conduct a constrained global optimization-based parameter estimation by using the machine learning-predicted K_m values as the reference values. The machine learning model achieved relatively good prediction scores: RMSE = 0.795 and R² = 0.536, making the subsequent global optimization easy and practical. The MLAGO approach reduced the error between simulation and experimental data while keeping K_m values close to the machine learning-predicted values. As a result, the MLAGO approach successfully estimated K_m values with less computational cost than the conventional method. Moreover, the MLAGO approach uniquely estimated K_m values, which were close to the measured values.

CONCLUSIONS

MLAGO overcomes the major problems in parameter estimation, accelerates kinetic modeling, and thus ultimately leads to better understanding of complex cellular systems. The web application for our machine learning-based K_m predictor is accessible at https://sites.google.com/view/kazuhiro-maeda/software-tools-web-apps , which helps modelers perform MLAGO on their own parameter estimation tasks.

The role of nitrogen-responsive regulators in controlling inorganic polyphosphate synthesis in Escherichia coli

Inorganic polyphosphate (polyP) is synthesized by bacteria under stressful environmental conditions and acts by a variety of mechanisms to promote cell survival. While the kinase that synthesizes polyP (PPK, encoded by the ppk gene) is well known, ppk transcription is not activated by environmental stress and little is understood about how environmental stress signals lead to polyP accumulation. Previous work has shown that the transcriptional regulators DksA, RpoN (σ54) and RpoE (σ24) positively regulate polyP production, but not ppk transcription, in Escherichia coli. In this work, we examine the role of the alternative sigma factor RpoN and nitrogen starvation stress response pathways in controlling polyP synthesis. We show that the RpoN enhancer binding proteins GlnG and GlrR impact polyP production, and uncover a new role for the nitrogen phosphotransferase regulator PtsN (EIIANtr) as a positive regulator of polyP production, acting upstream of DksA, downstream of RpoN and apparently independently of RpoE. However, neither these regulatory proteins nor common nitrogen metabolites appear to act directly on PPK, and the precise mechanism(s) by which polyP production is modulated after stress remain(s) unclear. Unexpectedly, we also found that the genes that impact polyP production vary depending on the composition of the rich media in which the cells were grown before exposure to polyP-inducing stress. These results constitute progress towards deciphering the regulatory networks driving polyP production under stress, and highlight the remarkable complexity of this regulation and its connections to a broad range of stress-sensing pathways.

Virtual metabolic human dynamic model for pathological analysis and therapy design for diabetes

A virtual metabolic human model is a valuable complement to experimental biology and clinical studies, because in vivo research involves serious ethical and technical problems. I have proposed a multi-organ and multi-scale kinetic model that formulates the reactions of enzymes and transporters with the regulation of hormonal actions at postprandial and postabsorptive states. The computational model consists of 202 ordinary differential equations for metabolites with 217 reaction rates and 1,140 kinetic parameter constants. It is the most comprehensive, largest, and highly predictive model of the whole-body metabolism. Use of the model revealed the mechanisms by which individual disorders, such as steatosis, β cell dysfunction, and insulin resistance, were combined to cause diabetes. The model predicted a glycerol kinase inhibitor to be an effective medicine for type 2 diabetes, which not only decreased hepatic triglyceride but also reduced plasma glucose. The model also enabled us to rationally design combination therapy.

•

A standard of virtual metabolic human dynamic models is proposed

•

It integrates the three scales of molecules, organs, and whole body

•

It gets insight into pathological mechanisms of type 1 and type 2 diabetes

•

It enables the computer-aided design of medication treatment for diabetes

ROS networks: designs, aging, Parkinson's disease and precision therapies

How the network around ROS protects against oxidative stress and Parkinson's disease (PD), and how processes at the minutes timescale cause disease and aging after decades, remains enigmatic. Challenging whether the ROS network is as complex as it seems, we built a fairly comprehensive version thereof which we disentangled into a hierarchy of only five simpler subnetworks each delivering one type of robustness. The comprehensive dynamic model described in vitro data sets from two independent laboratories. Notwithstanding its five-fold robustness, it exhibited a relatively sudden breakdown, after some 80 years of virtually steady performance: it predicted aging. PD-related conditions such as lack of DJ-1 protein or increased α-synuclein accelerated the collapse, while antioxidants or caffeine retarded it. Introducing a new concept (aging-time-control coefficient), we found that as many as 25 out of 57 molecular processes controlled aging. We identified new targets for "life-extending interventions": mitochondrial synthesis, KEAP1 degradation, and p62 metabolism.

A two-lane mechanism for selective biological ammonium transport

The transport of charged molecules across biological membranes faces the dual problem of accommodating charges in a highly hydrophobic environment while maintaining selective substrate translocation. This has been the subject of a particular controversy for the exchange of ammonium across cellular membranes, an essential process in all domains of life. Ammonium transport is mediated by the ubiquitous Amt/Mep/Rh transporters that includes the human Rhesus factors. Here, using a combination of electrophysiology, yeast functional complementation and extended molecular dynamics simulations, we reveal a unique two-lane pathway for electrogenic NH₄⁺ transport in two archetypal members of the family, the transporters AmtB from Escherichia coli and Rh50 from Nitrosomonas europaea. The pathway underpins a mechanism by which charged H⁺ and neutral NH₃ are carried separately across the membrane after NH₄⁺ deprotonation. This mechanism defines a new principle of achieving transport selectivity against competing ions in a biological transport process.

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.

Self-replenishment cycles generate a threshold response

Many metabolic cycles, including the tricarboxylic acid cycle, glyoxylate cycle, Calvin cycle, urea cycle, coenzyme recycling, and substrate cycles, are well known to catabolize and anabolize different metabolites for efficient energy and mass conversion. In terms of stoichiometric structure, this study explicitly identifies two types of metabolic cycles. One is the well-known, elementary cycle that converts multiple substrates into different products and recycles one of the products as a substrate, where the recycled substrate is supplied from the outside to run the cycle. The other is the self-replenishment cycle that merges multiple substrates into two or multiple identical products and reuses one of the products as a substrate. The substrates are autonomously supplied within the cycle. This study first defines the self-replenishment cycles that many scientists have overlooked despite its functional importance. Theoretical analysis has revealed the design principle of the self-replenishment cycle that presents a threshold response without any bistability nor cooperativity. To verify the principle, three detailed kinetic models of self-replenishment cycles embedded in an E. coli metabolic system were simulated. They presented the threshold response or digital switch-like function that steeply shift metabolic status.