A mathematical model of the folate cycle: new insights into folate homeostasis

Systems-wide analysis of A. fumigatus using kinetic modeling of metabolic pathways to identify putative drug targets

Aspergillosis is a major causative factor for morbidity in those with impaired immune systems, often caused by Aspergillus fumigatus. The diagnosis and treatment are difficult due to the diversity of individuals and risk factors and still pose a challenge for medical professionals. To understand the pathogenicity of any organism, it is critical to identify the significant metabolic pathways that are involved. Our work focused on developing kinetic models of critical pathways crucial for the survival of A. fumigatus using COPASI. While focusing on the folate biosynthesis, ergosterol biosynthesis and glycolytic pathway; sensitivity, time-course and steady-state analysis were performed to find the proteins/enzymes that are essential in the pathway and can be considered as potential drug targets. For further analysis of the interaction of drug targets identified, a protein-protein interaction (PPI) network was built, and hub nodes were identified using the Cytohubba package from Cytoscape. Based on the findings, dihydropteroate-synthase, dihydrofolate-reductase, 4-amino-4-deoxychorismate synthase, HMG-CoA-reductase, PG-isomerase and hexokinase could act as potential drug targets. Further, molecular docking and MM-GBSA analysis were performed with ligands chosen from DrugBank, and PubChem, and validated by experimental evidence and existing literature based on results from kinetic modeling and PPI network analysis. Based on docking scores and MM-GBSA results, molecular simulations were carried out for 1AJ2-dapsone, 1DIS-sulfamethazine, 1T02-lovastatin and 70YL-3-bromopyruvic acid complexes, which validated our findings. Our study provides a deeper insight into the mechanisms of A. fumigatus's metabolism to reveal dapsone, sulfamethazine, lovastatin and 3-bromopyruvic acid as potential drugs for the treatment of Aspergillosis.

The genetic landscape of a metabolic interaction

While much prior work has explored the constraints on protein sequence and evolution induced by physical protein-protein interactions, the sequence-level constraints emerging from non-binding functional interactions in metabolism remain unclear. To quantify how variation in the activity of one enzyme constrains the biochemical parameters and sequence of another, we focus on dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), a pair of enzymes catalyzing consecutive reactions in folate metabolism. We use deep mutational scanning to quantify the growth rate effect of 2696 DHFR single mutations in 3 TYMS backgrounds under conditions selected to emphasize biochemical epistasis. Our data are well-described by a relatively simple enzyme velocity to growth rate model that quantifies how metabolic context tunes enzyme mutational tolerance. Together our results reveal the structural distribution of epistasis in a metabolic enzyme and establish a foundation for the design of multi-enzyme systems.

The association of preconceptional folic acid supplementation with gestational age and preterm birth risk in twin pregnancies: a prospective cohort study in China

PURPOSE

This study aims to explore the association of maternal preconceptional folic acid (FA) supplementation with gestational age and preterm birth in twin pregnancies, and whether the association varies by chorionicity or conception mode.

METHODS

From November 2018 to December 2021, the information of FA supplementation and pregnancy outcomes were collected in twin pregnant women. The linear regression models and the logistic regression were used to test the association of preconceptional FA supplementation with gestational age at delivery and preterm birth and premature rupture of membranes (PROM).

RESULTS

A total of 416 twin pregnancies were included. Compared with no use in twins, maternal preconceptional FA use was associated with a 0.385-week longer gestational age (95% CI 0.019-0.751) and lower risk of preterm birth < 36 weeks (adjusted OR 0.519; 95% CI 0.301-0.895) and PROM (adjusted OR 0.426; 95% CI 0.215-0.845). The protective effect on preterm birth < 36 weeks and PROM is similar whether taking FA supplements alone or multivitamins. However, the associations varied by chorionicity and conception mode of twins or compliance with supplementation. The positive associations between preconceptional FA use and gestational age only remained significant among twins via assisted reproductive technology or dichorionic diamniotic twins. Significant protective effects on preterm birth < 36 weeks and PROM were only found among women who took FA at least 4 times a week before conception.

CONCLUSION

Maternal preconceptional FA supplementation was associated with longer gestation duration and lower risk of preterm birth < 36 weeks and PROM in twin pregnancies. To improve the success of their pregnancies, reproductive women should start taking FA supplements well before conception and with good compliance.

A new mathematical model of folate homeostasis in E. coli highlights the potential importance of the folinic acid futile cycle in cell growth

Folate (vitamin B9) plays a central role in one-carbon metabolism in prokaryotes and eukaryotes. This pathway mediates the transfer of one-carbon units, playing a crucial role in nucleotide synthesis, methylation, and amino acid homeostasis. The folinic acid futile cycle adds a layer of intrigue to this pathway, due to its associations with metabolism, cell growth, and dormancy. It also introduces additional complexity to folate metabolism. A logical way to deal with such complexity is to examine it by using mathematical modelling. This work describes the construction and analysis of a model of folate metabolism, which includes the folinic acid futile cycle. This model was tested under three in silico growth conditions. Model simulations revealed: 1) the folate cycle behaved as a stable biochemical system in three growth states (slow, standard, and rapid); 2) the initial concentration of serine had the greatest impact on metabolite concentrations; 3) 5-formyltetrahydrofolate cyclo-ligase (5-FCL) activity had a significant impact on the levels of the 7 products that carry the one-carbon donated from folates, and the redox couple NADP/NADPH; this was particularly evident in the rapid growth state; 4) 5-FCL may be vital to the survival of the cells by maintaining low levels of homocysteine, as high levels can induce toxicity; and 5) the antifolate therapeutic trimethoprim had a greater impact on folate metabolism with higher nutrient availability. These results highlight the important role of 5-FCL in intracellular folate homeostasis and mass generation under different metabolic scenarios.

Structural and functional analysis of two SHMT8 variants associated with soybean cyst nematode resistance

Two amino acid variants in soybean serine hydroxymethyltransferase 8 (SHMT8) are associated with resistance to the soybean cyst nematode (SCN), a devastating agricultural pathogen with worldwide economic impacts on soybean production. SHMT8 is a cytoplasmic enzyme that catalyzes the pyridoxal 5-phosphate-dependent conversion of serine and tetrahydrofolate (THF) to glycine and 5,10-methylenetetrahydrofolate. A previous study of the P130R/N358Y double variant of SHMT8, identified in the SCN-resistant soybean cultivar (cv.) Forrest, showed profound impairment of folate binding affinity and reduced THF-dependent enzyme activity, relative to the highly active SHMT8 in cv. Essex, which is susceptible to SCN. Given the importance of SCN-resistance in soybean agriculture, we report here the biochemical and structural characterization of the P130R and N358Y single variants to elucidate their individual effects on soybean SHMT8. We find that both single variants have reduced THF-dependent catalytic activity relative to Essex SHMT8 (10- to 50-fold decrease in kcat /Km ) but are significantly more active than the P130R/N368Y double variant. The kinetic data also show that the single variants lack THF-substrate inhibition as found in Essex SHMT8, an observation with implications for regulation of the folate cycle. Five crystal structures of the P130R and N358Y variants in complex with various ligands (resolutions from 1.49 to 2.30 Å) reveal distinct structural impacts of the mutations and provide new insights into allosterism. Our results support the notion that the P130R/N358Y double variant in Forrest SHMT8 produces unique and unexpected effects on the enzyme, which cannot be easily predicted from the behavior of the individual variants.

A Bloody Feast-Nutritional Regulation of Hematopoiesis

Hematopoietic stem cells provide us with a lifelong supply of blood cells. Hence, their proper function is absolutely essential for life, and their dysfunction can lead to infectious and malignant diseases. These cells have specific metabolic requirements to enable their lifelong function and blood-producing capacity. With the words of the Roman poet Juvenal "a healthy mind in a healthy body" in mind, it is intriguing to understand the connection between our daily diet and the quality of our blood, with the hope that through specific dietary adjustments we can improve our hematopoietic stem cell function and prevent disease. Nowadays, dietary supplements are an expanding market filled with potential and promises for better health. However, the link between many of those supplements and human physiology is obscure. Several groups have begun to shed light on this by investigating the metabolic regulation of hematopoiesis by specific nutrients. Beyond the link to dietary supplementation, these studies have also significantly improved our understanding of basic hematopoietic stem cell biology. Herein we summarize recent knowledge on the effect of specific vitamins and amino acids, which might be considered as dietary supplements, on normal hematopoiesis and hematopoietic stem cell function. We propose that improving our understanding of the link between nutrition in general and blood physiology can ultimately lead to the optimization of health-care policies, protocols, and standards of care.

The Genetic Landscape of a Metabolic Interaction

Enzyme abundance, catalytic activity, and ultimately sequence are all shaped by the need of growing cells to maintain metabolic flux while minimizing accumulation of deleterious intermediates. While much prior work has explored the constraints on protein sequence and evolution induced by physical protein-protein interactions, the sequence-level constraints emerging from non-binding functional interactions in metabolism remain unclear. To quantify how variation in the activity of one enzyme constrains the biochemical parameters and sequence of another, we focused on dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), a pair of enzymes catalyzing consecutive reactions in folate metabolism. We used deep mutational scanning to quantify the growth rate effect of 2,696 DHFR single mutations in 3 TYMS backgrounds under conditions selected to emphasize biochemical epistasis. Our data are well-described by a relatively simple enzyme velocity to growth rate model that quantifies how metabolic context tunes enzyme mutational tolerance. Together our results reveal the structural distribution of epistasis in a metabolic enzyme and establish a foundation for the design of multi-enzyme systems.

Homeostasis and injectivity: a reaction network perspective

Homeostasis represents the idea that a feature may remain invariant despite changes in some external parameters. We establish a connection between homeostasis and injectivity for reaction network models. In particular, we show that a reaction network cannot exhibit homeostasis if a modified version of the network (which we call homeostasis-associated network) is injective. We provide examples of reaction networks which can or cannot exhibit homeostasis by analyzing the injectivity of their homeostasis-associated networks.

One-carbon metabolism during the menstrual cycle and pregnancy

Many enzymes in one-carbon metabolism (OCM) are up- or down-regulated by the sex hormones which vary diurnally and throughout the menstrual cycle. During pregnancy, estradiol and progesterone levels increase tremendously to modulate physiological changes in the reproductive system. In this work, we extend and improve an existing mathematical model of hepatic OCM to understand the dynamic metabolic changes that happen during the menstrual cycle and pregnancy due to estradiol variation. In particular, we add the polyamine drain on S-adenosyl methionine and the direct effects of estradiol on the enzymes cystathionine β-synthase (CBS), thymidylate synthase (TS), and dihydrofolate reductase (DHFR). We show that the homocysteine concentration varies inversely with estradiol concentration, discuss the fluctuations in 14 other one-carbon metabolites and velocities throughout the menstrual cycle, and draw comparisons with the literature. We then use the model to study the effects of vitamin B12, vitamin B6, and folate deficiencies and explain why homocysteine is not a good biomarker for vitamin deficiencies. Additionally, we compute homocysteine throughout pregnancy, and compare the results with experimental data. Our mathematical model explains how numerous homeostatic mechanisms in OCM function and provides new insights into how homocysteine and its deleterious effects are influenced by estradiol. The mathematical model can be used by others for further in silico experiments on changes in one-carbon metabolism during the menstrual cycle and pregnancy.

Melatonin modulates metabolic remodeling in HNSCC by suppressing MTHFD1L-formate axis

Metabolic remodeling is now widely recognized as a hallmark of cancer, yet its role in head and neck squamous cell carcinoma (HNSCC) remains largely unknown. In this study, metabolomic analysis of melatonin-treated HNSCC cell lines revealed that exogenous melatonin inhibited many important metabolic pathways including folate cycle in HNSCC cells. Methylenetetrahydrofolate dehydrogenase 1 like (MTHFD1L), a metabolic enzyme of the folate cycle regulating the production of formate, was identified as a downstream target of melatonin. MTHFD1L was found to be markedly upregulated in HNSCC, and MTHFD1L overexpression was significantly associated with unfavorable clinical outcome of HNSCC patients. In addition, MTHFD1L promoted HNSCC progression in vitro and in vivo and reversed the oncostatic effects of exogenous melatonin. More importantly, the malignant phenotypes suppressed by knockdown of MTHFD1L or exogenous melatonin could be partially rescued by formate. Furthermore, we found that melatonin inhibited the expression of MTHFD1L in HNSCC cells through the downregulation of cyclic AMP-responsive element-binding protein 1 (CREB1) phosphorylation. Lastly, this novel regulatory axis of melatonin-p-CREB1-MTHFD1L-formate was also verified in HNSCC tissues. Collectively, our findings have demonstrated that MTHFD1L-formate axis promotes HNSCC progression and melatonin inhibits HNSCC progression through CREB1-mediated downregulation of MTHFD1L and formate. These findings have revealed new metabolic mechanisms in HNSCC and may provide novel insights on the therapeutic intervention of HNSCC.

The Gut Microbiome, Metformin, and Aging

Metformin has been extensively used for the treatment of type 2 diabetes, and it may also promote healthy aging. Despite its widespread use and versatility, metformin's mechanisms of action remain elusive. The gut typically harbors thousands of bacterial species, and as the concentration of metformin is much higher in the gut as compared to plasma, it is plausible that microbiome-drug-host interactions may influence the functions of metformin. Detrimental perturbations in the aging gut microbiome lead to the activation of the innate immune response concomitant with chronic low-grade inflammation. With the effectiveness of metformin in diabetes and antiaging varying among individuals, there is reason to believe that the gut microbiome plays a role in the efficacy of metformin. Metformin has been implicated in the promotion and maintenance of a healthy gut microbiome and reduces many age-related degenerative pathologies. Mechanistic understanding of metformin in the promotion of a healthy gut microbiome and aging will require a systems-level approach. Expected final online publication date for the Annual Review of Pharmacology and Toxicology, Volume 62 is January 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

One-carbon metabolism and global DNA methylation in mothers of individuals with Down syndrome

Down syndrome (DS) is the most common chromosomal disorder, resulting from the failure of normal chromosome 21 segregation. Studies have suggested that impairments within the one-carbon metabolic pathway can be of relevance for the global genome instability observed in mothers of individuals with DS. Based on the association between global DNA hypomethylation, genome instability, and impairments within the one-carbon metabolic pathway, the present study aimed to identify possible predictors, within the one-carbon metabolism, of global DNA methylation, measured by methylation patterns of LINE-1 and Alu repetitive sequences, in mothers of individuals with DS and mothers of individuals without the syndrome. In addition, we investigated one-carbon genetic polymorphisms and metabolites as maternal predisposing factors for the occurrence of trisomy 21 in children. Eighty-three samples of mothers of children with DS with karyotypically confirmed free trisomy 21 (case group) and 84 of mothers who had at least one child without DS or any other aneuploidy were included in the study. Pyrosequencing assays were performed to access global methylation. The results showed that group affiliation (case or control), betaine-homocysteine methyltransferase (BHMT) G742A and transcobalamin 2 (TCN2) C776G polymorphisms, and folate concentration were identified as predictors of global Alu DNA methylation values. In addition, thymidylate synthase (TYMS) 28-bp repeats 2R/3R or 3R/3R genotypes are independent maternal predisposing factors for having a child with DS. This study adds evidence that supports the association of impairments in the one-carbon metabolism, global DNA methylation, and the possibility of having a child with DS.

Metformin Is a Pyridoxal-5'-phosphate (PLP)-Competitive Inhibitor of SHMT2

Simple Summary

The mitochondrial enzyme serine hydroxymethyltransferase (SHMT2), which converts serine into glycine and generates 1C units for cell growth, is one of the most consistently overexpressed metabolic enzymes in cancer. Here, we reveal that the anti-diabetic biguanide metformin operates as a novel class of non-catalytic SHMT2 inhibitor that disrupts the pyridoxal-5′-phosphate (PLP)-dependent SHMT2 oligomerization process and ultimately SHMT2 activity. As SHMT2 inhibitors have not yet reached the clinic, these findings may aid the rational design of PLP-competitive SHMT2 inhibitors based on the biguanide skeleton of metformin.

Abstract

The anticancer actions of the biguanide metformin involve the functioning of the serine/glycine one-carbon metabolic network. We report that metformin directly and specifically targets the enzymatic activity of mitochondrial serine hydroxymethyltransferase (SHMT2). In vitro competitive binding assays with human recombinant SHMT1 and SHMT2 isoforms revealed that metformin preferentially inhibits SHMT2 activity by a non-catalytic mechanism. Computational docking coupled with molecular dynamics simulation predicted that metformin could occupy the cofactor pyridoxal-5′-phosphate (PLP) cavity and destabilize the formation of catalytically active SHMT2 oligomers. Differential scanning fluorimetry-based biophysical screening confirmed that metformin diminishes the capacity of PLP to promote the conversion of SHMT2 from an inactive, open state to a highly ordered, catalytically competent closed state. CRISPR/Cas9-based disruption of SHMT2, but not of SHMT1, prevented metformin from inhibiting total SHMT activity in cancer cell lines. Isotope tracing studies in SHMT1 knock-out cells confirmed that metformin decreased the SHMT2-channeled serine-to-formate flux and restricted the formate utilization in thymidylate synthesis upon overexpression of the metformin-unresponsive yeast equivalent of mitochondrial complex I (mCI). While maintaining its capacity to inhibit mitochondrial oxidative phosphorylation, metformin lost its cytotoxic and antiproliferative activity in SHMT2-null cancer cells unable to produce energy-rich NADH or FADH₂ molecules from tricarboxylic acid cycle (TCA) metabolites. As currently available SHMT2 inhibitors have not yet reached the clinic, our current data establishing the structural and mechanistic bases of metformin as a small-molecule, PLP-competitive inhibitor of the SHMT2 activating oligomerization should benefit future discovery of biguanide skeleton-based novel SHMT2 inhibitors in cancer prevention and treatment.

The structure of infinitesimal homeostasis in input-output networks

Homeostasis refers to a phenomenon whereby the output [Formula: see text] of a system is approximately constant on variation of an input [Formula: see text]. Homeostasis occurs frequently in biochemical networks and in other networks of interacting elements where mathematical models are based on differential equations associated to the network. These networks can be abstracted as digraphs [Formula: see text] with a distinguished input node [Formula: see text], a different distinguished output node o, and a number of regulatory nodes [Formula: see text]. In these models the input-output map [Formula: see text] is defined by a stable equilibrium [Formula: see text] at [Formula: see text]. Stability implies that there is a stable equilibrium [Formula: see text] for each [Formula: see text] near [Formula: see text] and infinitesimal homeostasis occurs at [Formula: see text] when [Formula: see text]. We show that there is an [Formula: see text] homeostasis matrix [Formula: see text] for which [Formula: see text] if and only if [Formula: see text]. We note that the entries in H are linearized couplings and [Formula: see text] is a homogeneous polynomial of degree [Formula: see text] in these entries. We use combinatorial matrix theory to factor the polynomial [Formula: see text] and thereby determine a menu of different types of possible homeostasis associated with each digraph [Formula: see text]. Specifically, we prove that each factor corresponds to a subnetwork of [Formula: see text]. The factors divide into two combinatorially defined classes: structural and appendage. Structural factors correspond to feedforward motifs and appendage factors correspond to feedback motifs. Finally, we discover an algorithm for determining the homeostasis subnetwork motif corresponding to each factor of [Formula: see text] without performing numerical simulations on model equations. The algorithm allows us to classify low degree factors of [Formula: see text]. There are two types of degree 1 homeostasis (negative feedback loops and kinetic or Haldane motifs) and there are two types of degree 2 homeostasis (feedforward loops and a degree two appendage motif).

Folylpoly-ɣ-glutamate synthetase association to the cytoskeleton: Implications to folate metabolon compartmentalization

Folates are essential for nucleotide biosynthesis, amino acid metabolism and cellular proliferation. Following carrier-mediated uptake, folates are polyglutamylated by folylpoly-ɣ-glutamate synthetase (FPGS), resulting in their intracellular retention. FPGS appears as a long isoform, directed to mitochondria via a leader sequence, and a short isoform reported as a soluble cytosolic protein (cFPGS). However, since folates are labile and folate metabolism is compartmentalized, we herein hypothesized that cFPGS is associated with the cytoskeleton, to couple folate uptake and polyglutamylation and channel folate polyglutamates to metabolon compartments. We show that cFPGS is a cytoskeleton-microtubule associated protein: Western blot analysis revealed that endogenous cFPGS is associated with the insoluble cellular fraction, i.e., cytoskeleton and membranes, but not with the cytosol. Mass spectrometry analysis identified the putative cFPGS interactome primarily consisting of microtubule subunits and cytoskeletal motor proteins. Consistently, immunofluorescence microscopy with cytosol-depleted cells demonstrated the association of cFPGS with the cytoskeleton and unconventional myosin-1c. Furthermore, since anti-microtubule, anti-actin cytoskeleton, and coatomer dissociation-inducing agents yielded perinuclear pausing of cFPGS, we propose an actin- and microtubule-dependent transport of cFPGS between the ER-Golgi and the plasma membrane. These novel findings support the coupling of folate transport with polyglutamylation and folate channeling to intracellular metabolon compartments. SIGNIFICANCE: FPGS, an essential enzyme catalyzing intracellular folate polyglutamylation and efficient retention, was described as a soluble cytosolic enzyme in the past 40 years. However, based on the lability of folates and the compartmentalization of folate metabolism and nucleotide biosynthesis, we herein hypothesized that cytoplasmic FPGS is associated with the cytoskeleton, to couple folate transport and polyglutamylation as well as channel folate polyglutamates to biosynthetic metabolon compartments. Indeed, using complementary techniques including Mass-spectrometry proteomics and fluorescence microscopy, we show that cytoplasmic FPGS is associated with the cytoskeleton and unconventional myosin-1c. This novel cytoskeletal localization of cytoplasmic FPGS supports the dynamic channeling of polyglutamylated folates to metabolon compartments to avoid oxidation and intracellular dilution of folates, while enhancing folate-dependent de novo biosynthesis of nucleotides and DNA/protein methylation.

S-adenosylmethionine: A metabolite critical to the regulation of autophagy

Autophagy is a mechanism that enables cells to maintain cellular homeostasis by removing damaged materials and mobilizing energy reserves in conditions of starvation. Although nutrient availability strongly impacts the process of autophagy, the specific metabolites that regulate autophagic responses have not yet been determined. Recent results indicate that S-adenosylmethionine (SAM) represents a critical inhibitor of methionine starvation-induced autophagy. SAM is primarily involved in four key metabolic pathways: transmethylation, transsulphuration, polyamine synthesis and 5'-deoxyadenosyl 5'-radical-mediated biochemical transformations. SAM is the sole methyl group donor involved in the methylation of DNA, RNA and histones, modulating the autophagic process by mediating epigenetic effects. Moreover, the metabolites of SAM, such as homocysteine, glutathione, decarboxylated SAM and spermidine, also exert important influences on the regulation of autophagy. From our perspective, nuclear-cytosolic SAM is a conserved metabolic inhibitor that connects cellular metabolic status and the regulation of autophagy. In the future, SAM might be a new target of autophagy regulators and be widely used in the treatment of various diseases.

Influence of Dietary Supplementation for Hyperhomocysteinemia Treatments

Hyperhomocysteinemia is recognized as risk factor for cardiovascular and age-associated diseases. Folic acid supplementation efficiently lowers plasma homocysteine (Hcy) levels, but high intake may negatively affect health because of unnatural levels of unmetabolized folic acid in the systemic circulation. Oxoproline (Oxo) provides by glutamic acid production an increase of intracellular folic acid trapping. Aim of this study was to compare the efficacy of three supplementation protocols: (1) traditional therapy (5-methyl-tetrahydrofolate: 15 mg/day); (2) 5 mL/day of Oxo with 300 μg folic acid (oxifolic); (3) 5 mL/day of Oxo alone (magnesio+) in a 90 days randomized trial on thirty-two moderate hyperhomocysteinemic (18.6 ± 2.4 μmol.L-1) patients (age 48 ± 14 yrs). Thiols: cysteine (Cys), cysteinylglycine (Cys-Gly) and glutathione levels were assessed too. Every supplementation induced significant (p range <0.05-0.0001) reductions of Hcy level and Cys concentration after the three protocols adopted. Otherwise glutathione concentration significantly increased after oxifolic (p < 0.01) and traditional (p < 0.05) supplementation. The integration of Oxo resulted an interesting alternative to traditional therapy because absence or minimal number of folates in the integrator eliminates any chance of excess that can constitute a long-term risk.

The association of prenatal folate and vitamin B12 levels with postnatal neurodevelopment varies by maternal MTHFR 677C>T genotype

Prenatal folate and vitamin B12 status have been linked to child neuropsychological development, but less is known about maternal genetic influences on this association. We conducted an exploratory longitudinal study of 181 mother–child pairs to assess whether maternal MTHFR 677C>T genotype modifies the association between maternal plasma folate and vitamin B12 in the first trimester of pregnancy and child neuropsychological development. Maternal plasma folate and vitamin B12 were determined by radioimmunoassay, and MTHFR 677C>T genotypes by PCR. We evaluated child neuropsychological development at 1, 3, 6, 12, 18, 24, and 30 month old using the Bayley Scales of Infant Development II. We analyzed the data using mixed-effects multivariate linear regression. The MTHFR 677C>T genotype distribution among the mothers was 18.2% CC, 49.8% CT, and 32.0% TT. The Mental Development Index (MDI) was inversely associated with maternal plasma folate among offspring of MTHFR 677CC mothers (β = -2.18 per twofold increase, 95% CI -4.07; -0.30, corrected P value = 0.02); no significant associations were observed among children born to women of other genotypes. The Motor Development Index (PDI) was not significantly associated with maternal plasma folate in any maternal MTHFR 677C>T genotype group, nor were MDI or PDI significantly associated with maternal plasma vitamin B12 in any maternal MTHFR 677C>T genotype group. This study suggests that maternal MTHFR 677CC genotype interacts with first-trimester plasma folate to influence offspring mental development.

Infinitesimal homeostasis in three-node input-output networks

Homeostasis occurs in a system where an output variable is approximately constant on an interval on variation of an input variable [Formula: see text]. Homeostasis plays an important role in the regulation of biological systems, cf.  Ferrell (Cell Syst 2:62-67, 2016), Tang and McMillen (J Theor Biol 408:274-289, 2016), Nijhout et al. (BMC Biol 13:79, 2015), and Nijhout et al. (Wiley Interdiscip Rev Syst Biol Med 11:e1440, 2018). A method for finding homeostasis in mathematical models is given in the control theory literature as points where the derivative of the output variable with respect to [Formula: see text] is identically zero. Such points are called perfect homeostasis or perfect adaptation. Alternatively, Golubitsky and Stewart (J Math Biol 74:387-407, 2017) use an infinitesimal notion of homeostasis (namely, the derivative of the input-output function is zero at an isolated point) to introduce singularity theory into the study of homeostasis. Reed et al. (Bull Math Biol 79(9):1-24, 2017) give two examples of infinitesimal homeostasis in three-node chemical reaction systems: feedforward excitation and substrate inhibition. In this paper we show that there are 13 different three-node networks leading to 78 three-node input-output network configurations, under the assumption that there is one input node, one output node, and they are distinct. The different configurations are based on which node is the input node and which node is the output node. We show nonetheless that there are only three basic mechanisms for three-node input-output networks that lead to infinitesimal homeostasis and we call them structural homeostasis, Haldane homeostasis, and null-degradation homeostasis. Substantial parts of this classification are given in Ma et al. (Cell 138:760-773, 2009) and Ferrell (2016) among others. Our contributions include giving a complete classification using general admissible systems (Golubitsky and Stewart in Bull Am Math Soc 43:305-364, 2006) rather than specific biochemical models, relating the types of infinitesimal homeostasis to the graph theoretic existence of simple paths, and providing the basis to use singularity theory to study higher codimension homeostasis singularities such as the chair singularities introduced in Nijhout and Reed (Integr Comp Biol 54(2):264-275, 2014. https://doi.org/10.1093/icb/icu010) and Nijhout et al. (Math Biosci 257:104-110, 2014). See Golubitsky and Stewart (2017). The first two of these mechanisms are illustrated by feedforward excitation and substrate inhibition. Structural homeostasis occurs only when the network has a feedforward loop as a subnetwork; that is, when there are two distinct simple paths connecting the input node to the output node. Moreover, when the network is just the feedforward loop motif itself, one of the paths must be excitatory and one inhibitory to support infinitesimal homeostasis. Haldane homeostasis occurs when there is a single simple path from the input node to the output node and then only when one of the couplings along this path has strength 0. Null-degradation homeostasis is illustrated by a biochemical example from Ma et al. (2009); this kind of homeostasis can occur only when the degradation constant of the third node is 0. The paper ends with an analysis of Haldane homeostasis infinitesimal chair singularities.

MTHFD1L-Mediated Redox Homeostasis Promotes Tumor Progression in Tongue Squamous Cell Carcinoma

Background: Routine changes in cell metabolism can drive tumor development, as the cellular program develops to promote glycolysis and redox homeostasis during tumor progression; however, the associated mechanisms in tongue squamous cell carcinoma (TSCC) remain unclear.

Methods: We investigated methylenetetrahydrofolate dehydrogenase 1-like (MTHFD1L) expression, its clinical relevance, redox modification, and molecular mechanisms using TSCC cells and tissues. The anti-tumor effects of MTHFD1L knockdown on TSCC tumorigenesis were evaluated in vitro and in vivo. Kaplan-Meier curves and the log-rank test were used to analyze disease-free survival and overall survival.

Results: TSCC patients with high expression levels of MTHFD1L had shorter overall survival (P < 0.05) and disease-free survival (P < 0.05). Knockdown of MTHFD1L reduced nicotinamide adenine dinucleotide phosphate (NADPH) levels and increased reactive oxygen species (ROS), which accelerated cell death under oxidative stress, such as hypoxia or glucose deprivation. Additionally, inhibition of MTHFD1L suppressed TSCC cell growth and delayed the cell cycle, including in xenograft experiments.

Conclusions: MTHFD1L confers redox homeostasis and promotes TSCC cell growth, which provides a great opportunity to study tumor metabolism in head and neck cancer. The mTORC1-4EBP1-eIF4E axis may affect the expression of MTHFD1L in TSCC. Inhibition of the expression of MTHFD1L may be an actionable and effective therapeutic target in TSCC.

Folic Acid and Risk of Preterm Birth: A Meta-Analysis

The results from epidemiologic studies linking blood folate concentrations, folic acid supplementation, or dietary folate to the risk of preterm birth are inconsistent. In this study, we aimed to summarize the available evidence on these associations. A systematic search of the PubMed/MEDLINE, Google Scholar, Web of Science, and Cochrane Library databases up to October 20, 2018 was performed and reference lists of retrieved articles were screened. Pooled odds ratios (ORs) and 95% confidence intervals (CIs) for the highest vs. the lowest levels of folate concentrations, folic acid supplementation, and dietary folate were calculated using random-effects models. Subgroup analyses and univariate meta-regression were performed to explore the sources of heterogeneity. Ten studies (six prospective cohort studies and four case-control studies) were included on folate concentrations, 13 cohort studies were included about folic acid supplementation, and 4 cohort studies were included regarding dietary folate intake. Higher maternal folate levels were associated with a 28% reduction in the risk of preterm birth (OR 0.72, 95% CI 0.56–0.93). Higher folic acid supplementation was associated with 10% lower risk of preterm birth (OR 0.90, 95% CI 0.85–0.95). In addition, a significant negative association was observed between dietary folate intake and the risk of preterm birth (OR 0.68, 95% CI 0.55–0.84), but no significant relation was seen between dietary folate and the risk of spontaneous preterm birth (OR 0.89, 95% CI 0.57–1.41). In the subgroup analysis, higher maternal folate levels in the third trimester were associated with a lower risk of preterm birth (OR 0.58, 95% CI 0.36–0.94). To initiate taking folic acid supplementation early before conception was adversely associated with preterm birth risk (OR 0.89, 95% CI 0.83–0.95). In conclusion, higher maternal folate levels and folic acid supplementation were significantly associated with a lower risk of preterm birth. The limited data currently available suggest that dietary folate is associated with a significantly decreased risk of preterm birth.

Genome-wide association study of circulating folate one-carbon metabolites

Experimental, observational, and clinical trials support a critical role of folate one-carbon metabolism (FOCM) in colorectal cancer (CRC) development. In this report, we focus on understanding the relationship between common genetic variants and metabolites of FOCM. We conducted a genome-wide association study of FOCM biomarkers among 1,788 unaffected (without CRC) individuals of European ancestry from the Colon Cancer Family Registry. Twelve metabolites, including 5-methyltetrahydrofolate, vitamin B₂ (flavin mononucleotide and riboflavin), vitamin B₆ (4-pyridoxic acid, pyridoxal, and pyridoxamine), total homocysteine, methionine, S-adenosylmethionine, S-adenosylhomocysteine, cystathionine, and creatinine were measured from plasma using liquid chromatography-mass spectrometry (LC-MS) or LC-MS/MS. For each individual biomarker, we estimated genotype array-specific associations followed by a fixed-effect meta-analysis. We identified the variant rs35976024 (at 2p11.2 and intronic of ATOH8) associated with total homocysteine (p = 4.9 × 10^-8 ). We found a group of six highly correlated variants on chromosome 15q14 associated with cystathionine (all p < 5 × 10^-8 ), with the most significant variant rs28391580 (p = 2.8 × 10^-8 ). Two variants (rs139435405 and rs149119426) on chromosome 14q13 showed significant (p < 5 × 10^-8 ) associations with S-adenosylhomocysteine. These three biomarkers with significant associations are closely involved in homocysteine metabolism. Furthermore, when assessing the principal components (PCs) derived from seven individual biomarkers, we identified the variant rs12665366 (at 6p25.3 and intronic of EXOC2) associated with the first PC (p = 2.3 × 10^-8 ). Our data suggest that common genetic variants may play an important role in FOCM, particularly in homocysteine metabolism.

Rat liver folate metabolism can provide an independent functioning of associated metabolic pathways

Folate metabolism in mammalian cells is essential for multiple vital processes, including purine and pyrimidine synthesis, histidine catabolism, methionine recycling, and utilization of formic acid. It remains unknown, however, whether these processes affect each other via folate metabolism or can function independently based on cellular needs. We addressed this question using a quantitative mathematical model of folate metabolism in rat liver cytoplasm. Variation in the rates of metabolic processes associated with folate metabolism (i.e., purine and pyrimidine synthesis, histidine catabolism, and influxes of formate and methionine) in the model revealed that folate metabolism is organized in a striking manner that enables activation or inhibition of each individual process independently of the metabolic fluxes in others. In mechanistic terms, this independence is based on the high activities of a group of enzymes involved in folate metabolism, which efficiently maintain close-to-equilibrium ratios between substrates and products of enzymatic reactions.

The Importance of Maternal Folate Status for Brain Development and Function of Offspring

The importance of an adequate periconceptional maternal folate status to prevent fetal neural tube defects has been well demonstrated and resulted in the recommendation for women to use folic acid supplements during the periconception period. The importance of maternal folate status for offspring neurodevelopment and brain health is less well described. We reviewed the current evidence linking maternal folate status before conception and during pregnancy with neurodevelopment and cognition of the offspring. We discuss both animal and human studies. Preclinical research revealed the importance of maternal folate status for several key processes required for normal neurodevelopment and brain functioning in the offspring, including DNA synthesis, regulation of gene expression, synthesis of phospholipids and neurotransmitters, and maintenance of healthy plasma homocysteine concentrations. Human observational studies are inconclusive; about half have shown a positive association between maternal folate status and cognitive performance of offspring. Whereas some studies suggest a positive association between maternal folate intake and cognition of offspring during childhood, data from interventional studies are too limited to conclude that there is a direct effect. Future preclinical studies are needed to help us characterize the behavioral effects, understand the underlying mechanisms, and to establish an optimal dosage and time window of folate supplementation. Moreover, more conclusive data from well-designed human observational studies and randomized controlled trials are needed to determine whether current recommendations for folic acid supplementation during pregnancy cover the needs for normal cognitive development in the offspring.

Computational modelling folate metabolism and DNA methylation: implications for understanding health and ageing

Dietary folates have a key role to play in health, as deficiencies in the intake of these B vitamins have been implicated in a wide variety of clinical conditions. The reason for this is folates function as single carbon donors in the synthesis of methionine and nucleotides. Moreover, folates have a vital role to play in the epigenetics of mammalian cells by supplying methyl groups for DNA methylation reactions. Intriguingly, a growing body of experimental evidence suggests that DNA methylation status could be a central modulator of the ageing process. This has important health implications because the methylation status of the human genome could be used to infer age-related disease risk. Thus, it is imperative we further our understanding of the processes which underpin DNA methylation and how these intersect with folate metabolism and ageing. The biochemical and molecular mechanisms, which underpin these processes, are complex. However, computational modelling offers an ideal framework for handling this complexity. A number of computational models have been assembled over the years, but to date, no model has represented the full scope of the interaction between the folate cycle and the reactions, which governs the DNA methylation cycle. In this review, we will discuss several of the models, which have been developed to represent these systems. In addition, we will present a rationale for developing a combined model of folate metabolism and the DNA methylation cycle.

Crowders Steal Dihydrofolate Reductase Ligands through Quinary Interactions

Dihydrofolate reductase (DHFR) reduces dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. Due to its role in one carbon metabolism, chromosomal DHFR is the target of the antibacterial drug, trimethoprim. Resistance to trimethoprim has resulted in a type II DHFR that is not structurally related to the chromosomal enzyme target. Because of its metabolic significance, understanding DHFR kinetics and ligand binding behavior in more cell-like conditions, where the total macromolecule concentration can be as great as 300 mg/mL, is important. The progress-curve kinetics and ligand binding properties of the drug target (chromosomal E. coli DHFR) and the drug resistant (R67 DHFR) enzymes were studied in the presence of macromolecular cosolutes. There were varied effects on NADPH oxidation and binding to the two DHFRs, with some cosolutes increasing affinity and others weakening binding. However, DHF binding and reduction in both DHFRs decreased in the presence of all cosolutes. The decreased binding of ligands is mostly attributed to weak associations with the macromolecules, as opposed to crowder effects on the DHFRs. Computer simulations found weak, transient interactions for both ligands with several proteins. The net charge of protein cosolutes correlated with effects on NADP⁺ binding, with near neutral and positively charged proteins having more detrimental effects on binding. For DHF binding, effects correlated more with the size of binding pockets on the protein crowders. These nonspecific interactions between DHFR ligands and proteins predict that the in vivo efficiency of DHFRs may be much lower than expected from their in vitro rates.

Sex differences in hepatic one-carbon metabolism

Background

There are large differences between men and women of child-bearing age in the expression level of 5 key enzymes in one-carbon metabolism almost certainly caused by the sex hormones. These male-female differences in one-carbon metabolism are greatly accentuated during pregnancy. Thus, understanding the origin and consequences of sex differences in one-carbon metabolism is important for precision medicine.

Results

We have created a mathematical model of hepatic one-carbon metabolism based on the underlying physiology and biochemistry. We use the model to investigate the consequences of sex differences in gene expression. We give a mechanistic understanding of observed concentration differences in one-carbon metabolism and explain why women have lower S-andenosylmethionine, lower homocysteine, and higher choline and betaine. We give a new explanation of the well known phenomenon that folate supplementation lowers homocysteine and we show how to use the model to investigate the effects of vitamin deficiencies, gene polymorphisms, and nutrient input changes.

Conclusions

Our model of hepatic one-carbon metabolism is a useful platform for investigating the mechanistic reasons that underlie known associations between metabolites. In particular, we explain how gene expression differences lead to metabolic differences between males and females.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-018-0621-7) contains supplementary material, which is available to authorized users.

In Vivo Titration of Folate Pathway Enzymes

How enzymes behave in cells is likely different from how they behave in the test tube. Previous in vitro studies find that osmolytes interact weakly with folate. Removal of the osmolyte from the solvation shell of folate is more difficult than removal of water, which weakens binding of folate to its enzyme partners. To examine if this phenomenon occurs in vivo, osmotic stress titrations were performed with Escherichia coli Two strategies were employed: resistance to an antibacterial drug and complementation of a knockout strain by the appropriate gene cloned into a plasmid that allows tight control of expression levels as well as labeling by a degradation tag. The abilities of the knockout and complemented strains to grow under osmotic stress were compared. Typically, the knockout strain could grow to high osmolalities on supplemented medium, while the complemented strain stopped growing at lower osmolalities on minimal medium. This pattern was observed for an R67 dihydrofolate reductase clone rescuing a ΔfolA strain, for a methylenetetrahydrofolate reductase clone rescuing a ΔmetF strain, and for a serine hydroxymethyltransferase clone rescuing a ΔglyA strain. Additionally, an R67 dihydrofolate reductase clone allowed E. coli DH5α to grow in the presence of trimethoprim until an osmolality of ∼0.81 is reached, while cells in a control titration lacking antibiotic could grow to 1.90 osmol.IMPORTANCEE. coli can survive in drought and flooding conditions and can tolerate large changes in osmolality. However, the cell processes that limit bacterial growth under high osmotic stress conditions are not known. In this study, the dose of four different enzymes in E. coli was decreased by using deletion strains complemented by the gene carried in a tunable plasmid. Under conditions of limiting enzyme concentration (lower than that achieved by chromosomal gene expression), cell growth can be blocked by osmotic stress conditions that are normally tolerated. These observations indicate that E. coli has evolved to deal with variations in its osmotic environment and that normal protein levels are sufficient to buffer the cell from environmental changes. Additional factors involved in the osmotic pressure response may include altered protein concentration/activity levels, weak solute interactions with ligands which can make it more difficult for proteins to bind their substrates/inhibitors/cofactors in vivo, and/or viscosity effects.

Identification of a co-target for enhancing efficacy of sorafenib in HCC through a quantitative modeling approach

Sorafenib (SFB), a multi-kinase inhibitor, is the only approved drug for treating hepatocellular carcinoma (HCC). However, SFB shows low efficacy in many cases. HCC related mortality therefore remains to be high worldwide. SFB, a multi-kinase inhibitor is also known to modulate the redox homeostasis in cancer cells. To understand the effect of SFB on the redox status, a quantitative understanding of the system is necessary. Kinetic modeling of the relevant pathways is a useful approach for obtaining a quantitative understanding of the pathway dynamics and to rank the individual factors based on the extent of influence they wield on the pathway. Here, we report a comprehensive model of the glutathione reaction network (GSH_net ), consisting of four modules and includes SFB-induced redox stress. We compared GSH_net simulations for HCC of six different etiologies with healthy liver, and correctly identified the expected variations in cancer. Next, we studied alterations induced in the system upon SFB treatment and observed differential H₂ O₂ dynamics in all the conditions. Using metabolic control analysis, we identified glutathione S-transferase (GST) as the enzyme with the highest selective control coefficient, making it an attractive co-target for potentiating the action of SFB across all six etiologies. As a proof-of-concept, we selected ethacrynic acid (EA), a known inhibitor of GST, and verified ex vivo that EA synergistically potentiates the cytotoxic effect of SFB. Being an FDA approved drug, EA is a promising candidate for repurposing as a combination therapy with SFB for HCC treatment.

Modelling and Simulation of Biochemical Processes Using Petri Nets

Systems composed of many components which interact with each other and lead to unpredictable global behaviour, are considered as complex systems. In a biological context, complex systems represent living systems composed of a large number of interacting elements. In order to study these systems, a precise mathematical modelling was typically used in this context. However, this modelling has limitations in the structural understanding and the behavioural study. In this sense, formal computational modelling is an approach that allows to model and to simulate dynamical properties of these particular systems. In this paper, we use Hybrid Functional Petri Net (HFPN), a Petri net extension dedicated to study and verify biopathways, to model and study the Methionine metabolic pathway. Methionine and its derivatives play significant roles in human bodies. We propose a set of simulations for the purpose of studying and analysing the Methionine pathway’s behaviour. Our simulation results have shown that several important abnormalities in this pathway are related to sever diseases such as Alzheimer’s disease, cardiovascular disease, cancers and others.

Selected B vitamins and their possible link to the aetiology of age-related sarcopenia: relevance of UK dietary recommendations

The possible roles of selected B vitamins in the development and progression of sarcopenia are reviewed. Age-related declines in muscle mass and function are associated with huge and increasing costs to healthcare providers. Falls and loss of mobility and independence due to declining muscle mass/function are associated with poor clinical outcomes and their prevention and management are attractive research targets. Nutritional status appears a key modifiable and affordable intervention. There is emerging evidence of sarcopenia being the result not only of diminished anabolic activity but also of declining neurological integrity in older age, which is emerging as an important aspect of the development of age-related decline in muscle mass/function. In this connection, several B vitamins can be viewed as not only cofactors in muscle synthetic processes, but also as neurotrophic agents with involvements in both bioenergetic and trophic pathways. The B vitamins thus selected are examined with respect to their relevance to multiple aspects of neuromuscular function and evidence is considered that requirements, intakes or absorption may be altered in the elderly. In addition, the evidence base for recommended intakes (UK recommended daily allowance) is examined with particular reference to original datasets and their relevance to older individuals. It is possible that inconsistencies in the literature with respect to the nutritional management of sarcopenia may, in part at least, be the result of compromised micronutrient status in some study participants. It is suggested that in order, for example, for intervention with amino acids to be successful, underlying micronutrient deficiencies must first be addressed/eliminated.

Joint effects of prenatal air pollutant exposure and maternal folic acid supplementation on risk of autism spectrum disorder

Independent studies report that periconceptional folic acid (FA) may decrease the risk of autism spectrum disorder (ASD) while exposure to air pollution may increase ASD risk. We examined the joint effects of gestational FA and air pollution exposures in association with ASD. We studied 346 ASD cases and 260 typically developing controls from the CHARGE case-control study. Self-reported FA intake for each month of pregnancy was quantified. Estimates of exposure to near roadway air pollution (NRP) and criteria air pollutant measures were assigned based on maternal residential history. Among mothers with high FA intake (>800 μg) in the first pregnancy month, exposure to increasing levels of all air pollutants, except ozone, during the first trimester was associated with decreased ASD risk, while increased ASD risk was observed for the same pollutant among mothers with low FA intake (≤800 μg). This difference was statistically significant for NO2 (e.g., NO2 and low FA intake: OR = 1.53 (0.91, 2.56) vs NO2 and high FA intake: OR = 0.74 (0.46, 1.19), P-interaction = 0.04). Mothers exposed to higher levels (≥ median) of any air pollutant during the first trimester of pregnancy and who reported low FA intake were at a higher ASD risk compared to mothers exposed to lower levels of that air pollutant and who reported high first month FA intake. Joint effects showed significant (alpha < 0.10) departures from expected interaction for NRP and NO2 . Our results suggest that periconceptional FA intake may reduce ASD risk in those with high prenatal air pollution exposure. Further study is needed to replicate these findings in larger sample sizes and to understand mechanisms of this potential relationship.. Autism Res 2018, 11: 69-80. © 2017 International Society for Autism Research, Wiley Periodicals, Inc.

LAY SUMMARY

We examined interactions between periconceptional folic acid (FA) and air pollution exposure on risk of ASD. Mothers exposed to higher levels of air pollution during the first trimester of pregnancy and who reported low supplemental FA intake during the first pregnancy month were at a higher ASD risk compared to mothers exposed to lower levels of air pollution and who reported high first month FA intake. Our results suggest that periconceptional FA intake may reduce ASD risk in those with high prenatal air pollution exposure.

Analysis of Homeostatic Mechanisms in Biochemical Networks

Cell metabolism is an extremely complicated dynamical system that maintains important cellular functions despite large changes in inputs. This "homeostasis" does not mean that the dynamical system is rigid and fixed. Typically, large changes in external variables cause large changes in some internal variables so that, through various regulatory mechanisms, certain other internal variables (concentrations or velocities) remain approximately constant over a finite range of inputs. Outside that range, the mechanisms cease to function and concentrations change rapidly with changes in inputs. In this paper we analyze four different common biochemical homeostatic mechanisms: feedforward excitation, feedback inhibition, kinetic homeostasis, and parallel inhibition. We show that all four mechanisms can occur in a single biological network, using folate and methionine metabolism as an example. Golubitsky and Stewart have proposed a method to find homeostatic nodes in networks. We show that their method works for two of these mechanisms but not the other two. We discuss the many interesting mathematical and biological questions that emerge from this analysis, and we explain why understanding homeostatic control is crucial for precision medicine.

Folate dietary insufficiency and folic acid supplementation similarly impair metabolism and compromise hematopoiesis

While dietary folate deficiency is associated with increased risk for birth defects and other diseases, evidence suggests that supplementation with folic acid can contribute to predisposition to some diseases, including immune dysfunction and cancer. Herein, we show that diets supplemented with folic acid both below and above the recommended levels led to significantly altered metabolism in multiple tissues in mice. Surprisingly, both low and excessive dietary folate induced similar metabolic changes, which were particularly evident for nucleotide biosynthetic pathways in B-progenitor cells. Diet-induced metabolic changes in these cells partially phenocopied those observed in mice treated with anti-folate drugs, suggesting that both deficiency and excessive levels of dietary folic acid compromise folate-dependent biosynthetic pathways. Both folate deficiency and excessive dietary folate levels compromise hematopoiesis, resulting in defective cell cycle progression, persistent DNA damage, and impaired production of lymphocytes. These defects reduce the reconstitution potential in transplantation settings and increase radiation-induced mortality. We conclude that excessive folic acid supplementation can metabolically mimic dietary folate insufficiency, leading to similar functional impairment of hematopoiesis.

Serum Vitamins and Minerals at Diagnosis and Follow-up in Children With Celiac Disease

OBJECTIVES

Children with celiac disease (CD) may experience deficiencies of several micronutrients. The objectives of the present study were to determine the prevalence of micronutrient deficiencies in children with CD at diagnosis, 6 months, and 18 months after the start of a gluten-free diet (GFD), and examine any correlation between micronutrient deficiencies, serum tissue transglutaminase (TtG) immunoglobulin A (IgA) antibody titers, and the degree of mucosal damage at diagnosis.

METHODS

Children (<17 years) with CD had their serum vitamins, minerals, and anti-TtG IgA antibodies measured at diagnosis, 6 and 18 months after starting a GFD. Histopathological changes of duodenal biopsies at diagnosis were documented using modified MARSH classification.

RESULTS

The medical records of 140 children (mean age at diagnosis 7.8 ± 4.01 years, 87 girls [621%]) with CD were examined. At diagnosis, serum vitamin D was the most commonly deficient vitamin in 70% of children. Serum ferritin was subnormal in 34.5% with zinc in 18.6% children but only 12 (10.9%) children had iron deficiency anemia. There was no correlation between micronutrient deficiencies at diagnosis and serum TtG IgA antibody titers or the degree of villous atrophy. The majority of serum levels of measured micronutrients had normalized after 6 months of starting GFD except for vitamin D, which improved but remained subnormal.

CONCLUSIONS

At diagnosis, most children with CD have vitamin D deficiency. The degree of micronutrient deficiencies does not correlate with the degree of villous atrophy or serum titers of anti-TtG IgA antibodies.

Robustness, flexibility, and sensitivity in a multifunctional motor control model

Motor systems must adapt to perturbations and changing conditions both within and outside the body. We refer to the ability of a system to maintain performance despite perturbations as "robustness," and the ability of a system to deploy alternative strategies that improve fitness as "flexibility." Different classes of pattern-generating circuits yield dynamics with differential sensitivities to perturbations and parameter variation. Depending on the task and the type of perturbation, high sensitivity can either facilitate or hinder robustness and flexibility. Here we explore the role of multiple coexisting oscillatory modes and sensory feedback in allowing multiphasic motor pattern generation to be both robust and flexible. As a concrete example, we focus on a nominal neuromechanical model of triphasic motor patterns in the feeding apparatus of the marine mollusk Aplysia californica. We find that the model can operate within two distinct oscillatory modes and that the system exhibits bistability between the two. In the "heteroclinic mode," higher sensitivity makes the system more robust to changing mechanical loads, but less robust to internal parameter variations. In the "limit cycle mode," lower sensitivity makes the system more robust to changes in internal parameter values, but less robust to changes in mechanical load. Finally, we show that overall performance on a variable feeding task is improved when the system can flexibly transition between oscillatory modes in response to the changing demands of the task. Thus, our results suggest that the interplay of sensory feedback and multiple oscillatory modes can allow motor systems to be both robust and flexible in a variable environment.

Methylfolate Trap Promotes Bacterial Thymineless Death by Sulfa Drugs

The methylfolate trap, a metabolic blockage associated with anemia, neural tube defects, Alzheimer’s dementia, cardiovascular diseases, and cancer, was discovered in the 1960s, linking the metabolism of folate, vitamin B₁₂, methionine and homocysteine. However, the existence or physiological significance of this phenomenon has been unknown in bacteria, which synthesize folate de novo. Here we identify the methylfolate trap as a novel determinant of the bacterial intrinsic death by sulfonamides, antibiotics that block de novo folate synthesis. Genetic mutagenesis, chemical complementation, and metabolomic profiling revealed trap-mediated metabolic imbalances, which induced thymineless death, a phenomenon in which rapidly growing cells succumb to thymine starvation. Restriction of B₁₂ bioavailability, required for preventing trap formation, using an “antivitamin B₁₂” molecule, sensitized intracellular bacteria to sulfonamides. Since boosting the bactericidal activity of sulfonamides through methylfolate trap induction can be achieved in Gram-negative bacteria and mycobacteria, it represents a novel strategy to render these pathogens more susceptible to existing sulfonamides.

Sulfonamides were the first agents to successfully treat bacterial infections, but their use later declined due to the emergence of resistant organisms. Restoration of these drugs may be achieved through inactivation of molecular mechanisms responsible for resistance. A chemo-genomic screen first identified 50 chromosomal loci representing the whole-genome antifolate resistance determinants in Mycobacterium smegmatis. Interestingly, many determinants resembled components of the methylfolate trap, a metabolic blockage exclusively described in mammalian cells. Targeted mutagenesis, genetic and chemical complementation, followed by chemical analyses established the methylfolate trap as a novel mechanism of sulfonamide sensitivity, ubiquitously present in mycobacteria and Gram-negative bacterial pathogens. Furthermore, metabolomic analyses revealed trap-mediated interruptions in folate and related metabolic pathways. These metabolic imbalances induced thymineless death, which was reversible with exogenous thymine supplementation. Chemical restriction of vitamin B₁₂, an important molecule required for prevention of the methylfolate trap, sensitized intracellular bacteria to sulfonamides. Thus, pharmaceutical promotion of the methylfolate trap represents a novel folate antagonistic strategy to render pathogenic bacteria more susceptible to available, clinically approved sulfonamides.

Homeostasis, singularities, and networks

Homeostasis occurs in a biological or chemical system when some output variable remains approximately constant as an input parameter [Formula: see text] varies over some interval. We discuss two main aspects of homeostasis, both related to the effect of coordinate changes on the input-output map. The first is a reformulation of homeostasis in the context of singularity theory, achieved by replacing 'approximately constant over an interval' by 'zero derivative of the output with respect to the input at a point'. Unfolding theory then classifies all small perturbations of the input-output function. In particular, the 'chair' singularity, which is especially important in applications, is discussed in detail. Its normal form and universal unfolding [Formula: see text] is derived and the region of approximate homeostasis is deduced. The results are motivated by data on thermoregulation in two species of opossum and the spiny rat. We give a formula for finding chair points in mathematical models by implicit differentiation and apply it to a model of lateral inhibition. The second asks when homeostasis is invariant under appropriate coordinate changes. This is false in general, but for network dynamics there is a natural class of coordinate changes: those that preserve the network structure. We characterize those nodes of a given network for which homeostasis is invariant under such changes. This characterization is determined combinatorially by the network topology.

Maternal folic acid supplementation modulates DNA methylation and gene expression in the rat offspring in a gestation period-dependent and organ-specific manner

Maternal folic acid supplementation can alter DNA methylation and gene expression in the developing fetus, which may confer disease susceptibility later in life. We determined which gestation period and organ were most sensitive to the modifying effect of folic acid supplementation during pregnancy on DNA methylation and gene expression in the offspring. Pregnant rats were randomized to a control diet throughout pregnancy; folic acid supplementation at 2.5× the control during the 1st, 2nd or 3rd week of gestation only; or folic acid supplementation throughout pregnancy. The brain, liver, kidney and colon from newborn pups were analyzed for folate concentrations, global DNA methylation and gene expression of the Igf2, Er-α, Gr, Ppar-α and Ppar-γ genes. Folic acid supplementation during the 2nd or 3rd week gestation or throughout pregnancy significantly increased brain folate concentrations (P<.001), while only folic acid supplementation throughout pregnancy significantly increased liver folate concentrations (P=.005), in newborn pups. Brain global DNA methylation incrementally decreased from early to late gestational folic acid supplementation and was the lowest with folic acid supplementation throughout pregnancy (P=.026). Folic acid supplementation in late gestation or throughout pregnancy significantly decreased Er-α, Gr and Ppar-α gene expression in the liver (P<.05). The kidney and colon were resistant to the effect of folic acid supplementation. Maternal folic acid supplementation affects tissue folate concentrations, DNA methylation and gene expression in the offspring in a gestation-period-dependent and organ-specific manner.

Mathematical Modeling of Glutathione Status in Type 2 Diabetics with Vitamin B12 Deficiency

Deficiencies in vitamin B12 and glutathione (GSH) are associated with a number of diseases including type 2 diabetes mellitus. We tested newly diagnosed Indian diabetic patients for correlation between their vitamin B12 and GSH, and found it to be weak. Here we seek to examine the theoretical dependence of GSH on vitamin B12 with a mathematical model of 1-carbon metabolism due to Reed and co-workers. We study the methionine cycle of the Reed-Nijhout model by developing a simple "stylized model" that captures its essential topology and whose kinetics are analytically tractable. The analysis shows-somewhat counter-intuitively-that the flux responsible for the homeostasis of homocysteine is, in fact, peripheral to the methionine cycle. Elevation of homocysteine arises from reduced activity of methionine synthase, a vitamin B12-dependent enzyme, however, this does not increase GSH biosynthesis. The model suggests that the lack of vitamin B12-GSH correlation is explained by suppression of activity in the trans-sulfuration pathway that limits the synthesis of cysteine and GSH from homocysteine. We hypothesize this "cysteine-block" is an essential consequence of vitamin B12 deficiency. It can be clinically relevant to appreciate that these secondary effects of vitamin B12 deficiency could be central to its pathophysiology.

A mathematical model of microbial folate biosynthesis and utilisation: implications for antifolate development

The metabolic biochemistry of folate biosynthesis and utilisation has evolved into a complex network of reactions. Although this complexity represents challenges to the field of folate research it has also provided a renewed source for antimetabolite targets. A range of improved folate chemotherapy continues to be developed and applied particularly to cancer and chronic inflammatory diseases. However, new or better antifolates against infectious diseases remain much more elusive. In this paper we describe the assembly of a generic deterministic mathematical model of microbial folate metabolism. Our aim is to explore how a mathematical model could be used to explore the dynamics of this inherently complex set of biochemical reactions. Using the model it was found that: (1) a particular small set of folate intermediates are overrepresented, (2) inhibitory profiles can be quantified by the level of key folate products, (3) using the model to scan for the most effective combinatorial inhibitions of folate enzymes we identified specific targets which could complement current antifolates, and (4) the model substantiates the case for a substrate cycle in the folinic acid biosynthesis reaction. Our model is coded in the systems biology markup language and has been deposited in the BioModels Database (MODEL1511020000), this makes it accessible to the community as a whole.

Mathematical analysis of the regulation of competing methyltransferases

BACKGROUND

Methyltransferase (MT) reactions, in which methyl groups are attached to substrates, are fundamental to many aspects of cell biology and human physiology. The universal methyl donor for these reactions is S-adenosylmethionine (SAM) and this presents the cell with an important regulatory problem. If the flux along one pathway is changed then the SAM concentration will change affecting all the other MT pathways, so it is difficult for the cell to regulate the pathways independently.

METHODS

We created a mathematical model, based on the known biochemistry of the folate and methionine cycles, to study the regulatory mechanisms that enable the cell to overcome this difficulty. Some of the primary mechanisms are long-range allosteric interactions by which substrates in one part of the biochemical network affect the activity of enzymes at distant locations in the network (not distant in the cell). Because of these long-range allosteric interactions, the dynamic behavior of the network is very complicated, and so mathematical modeling is a useful tool for investigating the effects of the regulatory mechanisms and understanding the complicated underlying biochemistry and cell biology.

RESULTS

We study the allosteric binding of 5-methyltetrahydrofolate (5 mTHF) to glycine-N-methyltransferase (GNMT) and explain why data in the literature implies that when one molecule binds, GNMT retains half its activity. Using the model, we quantify the effects of different regulatory mechanisms and show how cell processes would be different if the regulatory mechanisms were eliminated. In addition, we use the model to interpret and understand data from studies in the literature. Finally, we explain why a full understanding of how competing MTs are regulated is important for designing intervention strategies to improve human health.

CONCLUSIONS

We give strong computational evidence that once bound GNMT retains half its activity. The long-range allosteric interactions enable the cell to regulate the MT reactions somewhat independently. The low K m values of many MTs also play a role because the reactions then run near saturation and changes in SAM have little effect. Finally, the inhibition of the MTs by the product S-adenosylhomocysteine also stabilizes reaction rates against changes in SAM.