A mathematical model of the methionine cycle

Detailed resume of S-methyltransferases: Categories, structures, biological functions and research advancements in related pathophysiology and pharmacotherapy

Methylation is a vital chemical reaction in the metabolism of many drugs, neurotransmitters, hormones, and exogenous compounds. Among them, S-methylation plays a significant role in the biotransformation of sulfur-containing compounds, particularly chemicals with sulfhydryl groups. Currently, only three S-methyltransferases have been reported: thiopurine methyltransferase (TPMT), thiol methyltransferase (TMT), and thioether methyltransferase (TEMT). These enzymes are involved in various biological processes such as gene regulation, signal transduction, protein repair, tumor progression, and biosynthesis and degradation reactions in animals, plants, and microorganisms. Furthermore, they play pivotal roles in the metabolic pathways of essential drugs and contribute to the advancement of diseases such as tumors. This paper reviews the research progress on relevant structural features, metabolic mechanisms, inhibitor development, and influencing factors (gene polymorphism, S-adenosylmethionine level, race, sex, age, and disease) of S-methyltransferases. We hope that a better comprehension of S-methyltransferases will help to provide a reference for the development of novel strategies for related disorders and improve long-term efficacy.

Whole genome identification, molecular docking and expression analysis of enzymes involved in the selenomethionine cycle in Cardamine hupingshanensis

BACKGROUND

The selenomethionine cycle (SeMTC) is a crucial pathway for the metabolism of selenium. The basic bioinformatics and functions of four enzymes involved in the cycle including S-adenosyl-methionine synthase (MAT), SAM-dependent methyltransferase (MTase), S-adenosyl-homocysteine hydrolase (SAHH) and methionine synthase (MTR), have been extensively reported in many eukaryotes. The identification and functional analyses of SeMTC genes/proteins in Cardamine hupingshanensis and their response to selenium stress have not yet been reported.

RESULTS

In this study, 45 genes involved in SeMTC were identified in the C. hupingshanensis genome. Phylogenetic analysis showed that seven genes from ChMAT were clustered into four branches, twenty-seven genes from ChCOMT were clustered into two branches, four genes from ChSAHH were clustered into two branches, and seven genes from ChMTR were clustered into three branches. These genes were resided on 16 chromosomes. Gene structure and homologous protein modeling analysis illustrated that proteins in the same family are relatively conserved and have similar functions. Molecular docking showed that the affinity of SeMTC enzymes for selenium metabolites was higher than that for sulfur metabolites. The key active site residues identified for ChMAT were Ala269 and Lys273, while Leu221/231 and Gly207/249 were determined as the crucial residues for ChCOMT. For ChSAHH, the essential active site residues were found to be Asn87, Asp139 and Thr206/207/208/325. Ile204, Ser111/329/377, Asp70/206/254, and His329/332/380 were identified as the critical active site residues for ChMTR. In addition, the results of the expression levels of four enzymes under selenium stress revealed that ChMAT3-1 genes were upregulated approximately 18-fold, ChCOMT9-1 was upregulated approximately 38.7-fold, ChSAHH1-2 was upregulated approximately 11.6-fold, and ChMTR3-2 genes were upregulated approximately 28-fold. These verified that SeMTC enzymes were involved in response to selenium stress to varying degrees.

CONCLUSIONS

The results of this research are instrumental for further functional investigation of SeMTC in C. hupingshanensis. This also lays a solid foundation for deeper investigations into the physiological and biochemical mechanisms underlying selenium metabolism in plants.

Mass spectrometry as a lens into molecular human nutrition and health

Mass spectrometry (MS) has developed over the last decades into the most informative and versatile analytical technology in molecular and structural biology (). The platform enables discovery, identification, and characterisation of non-volatile biomolecules, such as proteins, peptides, DNA, RNA, nutrients, metabolites, and lipids at both speed and scale and can elucidate their interactions and effects. The versatility, robustness, and throughput have rendered MS a major research and development platform in molecular human health and biomedical science. More recently, MS has also been established as the central tool for 'Molecular Nutrition', enabling comprehensive and rapid identification and characterisation of macro- and micronutrients, bioactives, and other food compounds. 'Molecular Nutrition' thereby helps understand bioaccessibility, bioavailability, and bioefficacy of macro- and micronutrients and related health effects. Hence, MS provides a lens through which the fate of nutrients can be monitored along digestion via absorption to metabolism. This in turn provides the bioanalytical foundation for 'Personalised Nutrition' or 'Precision Nutrition' in which design and development of diets and nutritional products is tailored towards consumer and patient groups sharing similar genetic and environmental predisposition, health/disease conditions and lifestyles, and/or objectives of performance and wellbeing. The next level of integrated nutrition science is now being built as 'Systems Nutrition' where public and personal health data are correlated with life condition and lifestyle factors, to establish directional relationships between nutrition, lifestyle, environment, and health, eventually translating into science-based public and personal heath recommendations and actions. This account provides a condensed summary of the contributions of MS to a precise, quantitative, and comprehensive nutrition and health science and sketches an outlook on its future role in this fascinating and relevant field.

Identification and Mapping of QTLs for Sulfur-Containing Amino Acids in Soybean (Glycine max L.)

Soybean is a major source of high-quality protein for humans and animals. The content of sulfur-containing amino acids (SAA) in soybean is insufficient, which has become the main factor limiting soybean nutrition. In this study, we used the high-density genetic maps derived from Guizao 1 and Brazil 13 to evaluate the quantitative trait loci of cysteine (Cys), methionine (Met), SAA, glycinin (7S), β-conglycinin (11S), ratio of glycinin to β-conglycinin (RGC), and protein content (PC). In genetic map linkage analysis, the major and stable 44 QTLs were detected, which shared nine bin intervals. Among them, the bin interval (bin157-bin160) on chromosome 5 was detected in multiple environments as a stable QTL, which was linked to 11S, 7S, RGC, and SSA. Based on the analysis of bioinformatics and RNA-sequencing data, 16 differentially expressed genes (DEGs) within these QTLs were selected as candidate genes. These results will help to elucidate the genetic mechanism of soybean SAA-related traits and provide the basis for the gene mining of sulfur-containing amino acids.

The Association of Homocysteine and Diabetic Retinopathy in Homocysteine Cycle in Chinese Patients With Type 2 Diabetes

OBJECTIVE

This study aimed to explore the relationship between homocysteine (Hcy) and diabetic retinopathy (DR) and the impacts of the Hcy pathway on this relationship against this background.

METHODS

This study retrieved 1979 patients with type 2 diabetes (T2D) from the First Affiliated Hospital of Liaoning Medical University in Jinzhou, Liaoning Province, China. Multiple logistic regression was used to analyze the effects of Hcy cycle on the relationship between Hcy and DR. Spearman's rank correlation analysis was used to analyze the correlation between risk factors related to DR progression and Hcy. Finally, the results of logistic regression were supplemented by mediation analysis.

RESULTS

We found there was a negative correlation between low concentration of Hcy and DR (OR : 0.83, 95%CI: 0.69-1). After stratifying all patients by cysteine (Cys) or Methionine (Met), this relationship remained significant only in low concentration of Cys (OR: 0.75, 95%CI: 0.61-0.94). Through the RCS curve, we found that the effect of Hcy on DR presents a U-shaped curve relationship. Mediating effect in Met and Hcy cycles was also significant [Total effect c (OR: 0.968, 95%CI: 0.938-0.998), Direct effect path c' (OR: 0.969, 95%CI: 0.940-0.999), Path a (OR: 1.047, 95%CI: 1.004-1.091), Path b (OR: 0.964, 95%CI: 0.932-0.998)].

CONCLUSIONS

The relationship between Hcy and DR presents a U-shaped curve and the homocysteine cycle pathway has an impact on it. And too low concentration of Hcy indicates a lack of other substances, such as vitamins. It is suggested that the progression of DR is the result of a combination of many risk factors. Further prospective studies are needed to determine the role of Hcy in the pathogenesis of DR.

S-adenosyl-L-homocysteine hydrolase FgSah1 is required for fungal development and virulence in Fusarium graminearum

The S-adenosyl-L-homocysteine hydrolase (Sah1) plays a crucial role in methylation and lipid metabolism in yeast and mammals, yet its function remains elusive in filamentous fungi. In this study, we characterized Sah1 in the phytopathogenic fungus F. graminearum by generating knockout and knockout-complemented strains of FgSAH1. We found that the FgSah1-GFP fusion protein was localized to the cytoplasm, and that deletion of FgSAH1 resulted in defects in vegetative growth, asexual and sexual reproduction, stress responses, virulence, lipid metabolism, and tolerance against fungicides. Moreover, the accumulations of S-adenosyl-L-homocysteine (AdoHcy) and S-adenosyl-L-methionine (AdoMet) (the methyl group donor in most methyl transfer reactions) in ΔFgSah1 were seven- and ninefold higher than those in the wild-type strain, respectively. All of these defective phenotypes in ΔFgSah1 mutants were rescued by target gene complementation. Taken together, these results demonstrate that FgSah1 plays essential roles in methylation metabolism, fungal development, full virulence, multiple stress responses, lipid metabolism, and fungicide sensitivity in F. graminearum. To our knowledge, this is the first report on the systematic functional characterization of Sah1 in F. graminearum.

Interspecific Variation in One-Carbon Metabolism within the Ovarian Follicle, Oocyte, and Preimplantation Embryo: Consequences for Epigenetic Programming of DNA Methylation

One-carbon (1C) metabolism provides methyl groups for the synthesis and/or methylation of purines and pyrimidines, biogenic amines, proteins, and phospholipids. Our understanding of how 1C pathways operate, however, pertains mostly to the (rat) liver. Here we report that transcripts for all bar two genes (i.e., BHMT, MAT1A) encoding enzymes in the linked methionine-folate cycles are expressed in all cell types within the ovarian follicle, oocyte, and blastocyst in the cow, sheep, and pig; as well as in rat granulosa cells (GCs) and human KGN cells (a granulosa-like tumor cell line). Betaine-homocysteine methyltransferase (BHMT) protein was absent in bovine theca and GCs, as was activity of this enzyme in GCs. Mathematical modeling predicted that absence of this enzyme would lead to more volatile S-adenosylmethionine-mediated transmethylation in response to 1C substrate (e.g., methionine) or cofactor provision. We tested the sensitivity of bovine GCs to reduced methionine (from 50 to 10 µM) and observed a diminished flux of 1C units through the methionine cycle. We then used reduced-representation bisulfite sequencing to demonstrate that this reduction in methionine during bovine embryo culture leads to genome-wide alterations to DNA methylation in >1600 genes, including a cohort of imprinted genes linked to an abnormal fetal-overgrowth phenotype. Bovine ovarian and embryonic cells are acutely sensitive to methionine, but further experimentation is required to determine the significance of interspecific variation in BHMT expression.

Biomedical applications of methionine-based systems

Methionine (Met), an essential amino acid in the human body, possesses versatile features based on its chemical modification, cell metabolism and metabolic derivatives. Benefitting from its multifunctional properties, Met holds immense potential for biomedical applications. In this review, we systematically summarize the recent progress in Met-based strategies for biomedical applications. First, given the unique structural characteristics of Met, two chemical modification methods are briefly introduced. Subsequently, due to the disordered metabolic state of tumor cells, applications of Met in cancer treatment and diagnosis are summarized in detail. Furthermore, the efficacy of S-adenosylmethionine (SAM), as the most important metabolic derivative of Met, for treating liver diseases is mentioned. Finally, we analyze the current challenges and development trends of Met in the biomedical field, and suggest that Met-restriction therapy might be a promising approach to treat COVID-19.

Effects of guanidinoacetic acid and coated folic acid supplementation on growth performance, nutrient digestion and hepatic gene expression in Angus bulls

To evaluate the impacts of guanidinoacetic acid (GAA) and coated folic acid (CFA) on growth performance, nutrient digestion and hepatic gene expression, fifty-two Angus bulls were assigned to four groups in a 2 × 2 factor experimental design. The CFA of 0 or 6 mg/kg dietary DM folic acid was supplemented in diets with GAA of 0 (GAA-) or 0·6 g/kg DM (GAA+), respectively. Average daily gain (ADG), feed efficiency and hepatic creatine concentration increased with GAA or CFA addition, and the increased magnitude of these parameters was greater for addition of CFA in GAA- diets than in GAA+ diets. Blood creatine concentration increased with GAA or CFA addition, and greater increase was observed when CFA was supplemented in GAA+ diets than in GAA- diets. DM intake was unchanged, but rumen total SCFA concentration and digestibilities of DM, crude protein, neutral-detergent fibre and acid-detergent fibre increased with the addition of GAA or CFA. Acetate:propionate ratio was unaffected by GAA, but increased for CFA addition. Increase in blood concentrations of albumin, total protein and insulin-like growth factor-1 (IGF-1) was observed for GAA or CFA addition. Blood folate concentration was decreased by GAA, but increased with CFA addition. Hepatic expressions of IGF-1, phosphoinositide 3-kinase, protein kinase B, mammalian target of rapamycin and ribosomal protein S6 kinase increased with GAA or CFA addition. Results indicated that the combined supplementation of GAA and CFA could not cause ADG increase more when compared with GAA or CFA addition alone.

Biodegradation of typical BFRs 2,4,6-tribromophenol by an indigenous strain Bacillus sp. GZT isolated from e-waste dismantling area through functional heterologous expression

Legacy wastewater contaminants from e-waste dismantling process such as 2,4,6-tribromophenol (TBP), one of the most widely used brominated flame retardants (BFRs), have raised concern owing to their toxicity and recalcitrance. Our previously isolated Bacillus sp. GZT from river sludge in e-waste dismantling area is a good candidate for bioremediation of BFRs contaminated sites considering its remarkable degradability of TBP and its intermediates. However, there exists a new challenge because bio-degrader cannot produce enough biomass or metabolic activity to cleanup TBP in practice complex environment. Here, we heterologously expressed and functionally characterized the genes and enzymes responsible for TBP degradation to examine the feasibility of enhancing the ability of this microorganism to detoxify TBP. Results demonstrated that five recombinant strains containing functional genes, designated tbpA, tbpB, tbpC, tbpD, and tbpE, become more tolerant toward a wide range of brominated compounds than the nontransgenic strain. Cytochrome P450 reductase encoded by tbpA gene could greatly increase efficiency to remove TBP (98.8%), as compared to wild-type strain GZT (93.2%). Its debromination intermediates 2,4-dibromophenol, 2,6-dibromo-4-methylphenol and 2-bromophenol were significantly metabolized by halophenol dehalogenases encoded by tbpB, tbpC, and tbpD, respectively. Finally, under the function of tbpE gene encoding enzyme, further debrominated product (phenol) was dramatically detoxified. To reduce the risk of these xenobiotics, the expression of these genes can be induced and significantly up-regulated during exposure to them. These results open broad scope for future study in developing genetic engineering technologies for more efficient remediation wastewater of e-waste recycling sites contaminated with TBP, which would certainly be important steps to lower TBP exposures and prevent potential health effects.

Perspective: Advancing Understanding of Population Nutrient-Health Relations via Metabolomics and Precision Phenotypes

Diet and lifestyle are vital to population health, but their true contribution is difficult to quantify using traditional methods. Nutrient-health relations are typically based on epidemiological associations that are assessed at the population level, traditionally using self-reported dietary and lifestyle data. Unfortunately, such measures are inherently inaccurate. New technologies such as metabolomics can measure nutritional and micronutrient profiles in body fluids, providing objective evaluation of nutritional status. A critical step toward accurate health prediction models would be the building of integrated repositories of nutritional measures combining subjective methods of reporting with objective metabolomics profiles and precise phenotypic data. Here we outline a roadmap to achieve this goal and discuss both the advantages and risks of this approach. We also highlight the uncertain associations between the complexity of high-dimensional data generated in 'omics research (along with the public confusion this may engender) and the rapid adoption of 'omics approaches by nutrition and health companies to develop nutritional products and services.

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.

Identification of evolutionary and kinetic drivers of NAD-dependent signaling

Nicotinamide adenine dinucleotide (NAD) provides an important link between metabolism and signal transduction and has emerged as central hub between bioenergetics and all major cellular events. NAD-dependent signaling (e.g., by sirtuins and poly-adenosine diphosphate [ADP] ribose polymerases [PARPs]) consumes considerable amounts of NAD. To maintain physiological functions, NAD consumption and biosynthesis need to be carefully balanced. Using extensive phylogenetic analyses, mathematical modeling of NAD metabolism, and experimental verification, we show that the diversification of NAD-dependent signaling in vertebrates depended on 3 critical evolutionary events: 1) the transition of NAD biosynthesis to exclusive usage of nicotinamide phosphoribosyltransferase (NamPT); 2) the occurrence of nicotinamide N-methyltransferase (NNMT), which diverts nicotinamide (Nam) from recycling into NAD, preventing Nam accumulation and inhibition of NAD-dependent signaling reactions; and 3) structural adaptation of NamPT, providing an unusually high affinity toward Nam, necessary to maintain NAD levels. Our results reveal an unexpected coevolution and kinetic interplay between NNMT and NamPT that enables extensive NAD signaling. This has implications for therapeutic strategies of NAD supplementation and the use of NNMT or NamPT inhibitors in disease treatment.

Transcriptomic and proteomic profiling revealed global changes in Streptococcus thermophilus during pH-controlled batch fermentations

Understanding global changes of physiological processes at the molecular level during the growth of Streptococcus thermophilus is essential for the rational design of cultivation media and the optimization of bioprocesses. Transcriptomics and proteomics were combined to investigate the global changes at the transcript and protein level during the growth of S. thermophilus. The expression of 1396 genes (FDR ≤ 0.001) and 876 proteins (P < 0.05) changed significantly over time. The most remarkable growth phase dependent changes occurred in the late-lag phase and were related to heterofermentation, glycolysis, peptidoglycan biosynthesis, conversion between amino acids and stress response. The present results could provide theoretical guidance for high-cell-density culture, help design cultivation media, and help attain a high biomass of S. thermophilus.

Mitochondrial function in immature bovine oocytes is improved by an increase of cellular cyclic AMP

Although in vitro maturation (IVM) of oocytes is important for assisted reproduction, the rate of development of embryos from IVM oocytes is lower than from their in vivo counterparts. It has been shown that an artificial increase of intracellular cAMP before culture significantly improves oocyte developmental competence in cattle and mice. Here, we revealed that forskolin and 3-isobutyl-1-methylxanthine treatment of prophase-stage oocytes induced the expression of genes required for glycolysis, fatty acid degradation, and the mitochondrial electron transport system and improved mitochondrial functions and ATP levels in oocytes without involving nuclear maturation. We propose the existence of a comprehensive energy-supply system in oocytes under follicle-stimulating hormone stimulation as a potential explanation of how oocytes acquire developmental competence.

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.

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.

Folate promotes S-adenosyl methionine reactions and the microbial methylation cycle and boosts ruminants production and reproduction

Folate has gained significant attention due to its vital role in biological methylation and epigenetic machinery. Folate, or vitamin (B₉), is only produced through a de novo mechanism by plants and micro-organisms in the rumen of mature animals. Although limited research has been conducted on folate in ruminants, it has been noted that ruminal synthesis could not maintain folate levels in high yielding dairy animals. Folate has an essential role in one-carbon metabolism and is a strong antiproliferative agent. Folate increases DNA stability, being crucial for DNA synthesis and repair, the methylation cycle, and preventing oxidation of DNA by free radicals. Folate is also critical for cell division, metabolism of proteins, synthesis of purine and pyrimidine, and increasing the de novo delivery of methyl groups and S-adenosylmethionine. However, in ruminants, metabolism of B₁₂ and B₉ vitamins are closely connected and utilization of folate by cells is significantly affected by B₁₂ vitamin concentration. Supplementation of folate through diet, particularly in early lactation, enhanced metabolic efficiency, lactational performance, and nutritional quality of milk. Impaired absorption, oxidative degradation, or deficient supply of folate in ruminants affects DNA stability, cell division, homocysteine remethylation to methionine, de novo synthesis of S-adenosylmethionine, and increases DNA hypomethylation, uracil misincorporation into DNA, chromosomal damage, abnormal cell growth, oxidative species, premature birth, low calf weight, placental tube defects, and decreases production and reproduction of ruminant animals. However, more studies are needed to overcome these problems and reduce enormous dietary supplement waste and impaired absorption of folate in ruminants. This review was aimed to highlight the vital role of folic acid in ruminants performance.

Prediction of Methionine and Homocysteine levels in Zucker diabetic fatty (ZDF) rats as a T2DM animal model after consumption of a Methionine-rich diet

Background

Although alterations in the methionine metabolism cycle (MMC) have been associated with vascular complications of diabetes, there have not been consistent results about the levels of methionine and homocysteine in type 2 diabetes mellitus (T2DM). The aim of the current study was to predict changes in plasma methionine and homocysteine concentrations after simulated consumption of methionine-rich foods, following the development of a mathematical model for MMC in Zucker Diabetic Fatty (ZDF) rats, as a representative T2DM animal model.

Method

The model building and simulation were performed using NONMEM® (ver. 7.3.0) assisted by Perl-Speaks-NONMEM (PsN, ver. 4.3.0). Model parameters were derived using first-order conditional estimation method with interactions permitted among the parameters (FOCE-INTER). NCA was conducted using Phoenix (ver. 6.4.0). For all tests, we considered a P-value < 0.05 to reflect statistical significance.

Results

Our model featured seven compartments that considered all parts of the cycle by applying non-linear mixed effects model. Conversion of S-adenosyl-L-homocysteine (SAH) to homocysteine increased and the metabolism of homocysteine was reduced under diabetic conditions, and consequently homocysteine accumulated in the elimination phase.

Using our model, we performed simulations to compare the changes in plasma methionine and homocysteine concentrations between ZDF and normal rats, by multiple administrations of the methionine-rich diet of 1 mmol/kg, daily for 60 days. The levels of methionine and homocysteine were elevated approximately two- and three-fold, respectively, in ZDF rats, while there were no changes observed in the normal control rats.

Conclusion

These results can be interpreted to mean that both methionine and homocysteine will accumulate in patients with T2DM, who regularly consume high-methionine foods.

Neuro-fuzzy model of homocysteine metabolism

In view of well-documented association of hyperhomocysteinaemia with a wide spectrum of diseases and higher incidence of vitamin deficiencies in Indians, we proposed a mathematical model to forecast the role of demographic and genetic variables in influencing homocysteinemetabolism and investigated the influence of life style modulations in controlling homocysteine levels. Total plasma homocysteine levels were measured in fasting samples using reverse phase HPLC. Multiple linear regression (MLR) and neuro-fuzzy models were developed. The MLR model explained 64% variability in homocysteine, while the neurofuzzy model showed higher accuracy in predicting homocysteine with a mean absolute error of 0.00002 μmol/L. Methylene tetrahydrofolate reductase (MTHFR) C677T, 5-methyltetrahydrofolate homocysteine methyltransferase (MTR) A2756G and 5- methyltetrahydrofolate homocysteine methyltransferase reductase (MTRR) A66G were shown to be positively associatiated with homocysteine, while nonvegetarian diet, serine hydroxymethyltransferase 1 (SHMT1) C1420T and TYMS 5'-UTR 28 bp tandem repeat exhibited negative association with homocysteine. The protective role of SHMT1 C1420T was attributed to more H-bonding interactions in the mutant modelled compared to the wild type, as shown through in silico analysis. To conclude, polymorphisms in genes regulating remethylation of homocysteine strongly influence homocysteine levels. The restoration of one-carbon homeostasis by SHMT1 C1420T or increased flux of folate towards remethylation due to TYMS 5'-UTR 28 bp tandem repeat or nonvegetarian diet can lower homocysteine levels.

Essential role of methyl donors in animal productivity

Dietary requirements for the methyl donors, choline, betaine and folate, in livestock species are poorly defined and have not been included in diet formulation software or simulation models for animals. A deficiency of methyl donors may promote an inflammatory state, which is significant for the livestock industry because chronic low-grade inflammation is widespread among livestock under commercial conditions. Furthermore, recent evidence showing that methyl donors activate adenosine monophosphate-activated protein kinase, an anti-inflammatory master switch, indicates that dietary methyl-donor supplementation could be used to prevent or ameliorate chronic inflammation and its sequelae in livestock, which include fatty liver disease in dairy cows, fatty liver and kidney syndrome in broilers, fatty liver haemorrhagic syndrome in layers, gut ulcers in pigs, liver abscesses in feedlot cattle, enteritis in poultry and susceptibility to heat stress in all species. Because of the complexity of interactions among methyl donors, a modelling approach inclusive of a supporting research effort will be required to harness the potential of methyl-donor supplementation in livestock production.

Role of methionine adenosyltransferase 2A in bovine preimplantation development and its associated genomic regions

Methionine adenosyltransferase (MAT) is involved in folate-mediated one-carbon metabolism, which is essential for preimplantation embryos in terms of both short-term periconceptional development and long-term phenotypic programming beyond the periconceptional period. Here, our immunofluorescence analysis of bovine oocytes and preimplantation embryos revealed the consistent expression of MAT2A (the catalytic subunit of the ubiquitously expressed-type of MAT isozyme) during this period. Addition of the MAT2A inhibitor FIDAS to the culture media of bovine preimplantation embryos reduced their blastocyst development, revealing the particular importance of MAT2A in successful blastocyst development. Exploration of MAT2A-associated genomic regions in bovine blastocysts using chromatin immunoprecipitation and sequencing (ChIP-seq) identified candidate MAT2A-associated genes implicated not only in short-term periconceptional embryo development, but also in long-term phenotypic programming during this period in terms of growth, metabolism, and immune functions. These results suggest the critical involvement of MAT2A in the periconceptional period in life-long programming of health and disease as well as successful preimplantation development.

A hybrid stochastic model of folate-mediated one-carbon metabolism: Effect of the common C677T MTHFR variant on de novo thymidylate biosynthesis

Folate-mediated one-carbon metabolism (FOCM) is an interconnected network of metabolic pathways, including those required for the de novo synthesis of dTMP and purine nucleotides and for remethylation of homocysteine to methionine. Mouse models of folate-responsive neural tube defects (NTDs) indicate that impaired de novo thymidylate (dTMP) synthesis through changes in SHMT expression is causative in folate-responsive NTDs. We have created a hybrid computational model comprised of ordinary differential equations and stochastic simulation. We investigated whether the de novo dTMP synthesis pathway was sensitive to perturbations in FOCM that are known to be associated with human NTDs. This computational model shows that de novo dTMP synthesis is highly sensitive to the common MTHFR C677T polymorphism and that the effect of the polymorphism on FOCM is greater in folate deficiency. Computational simulations indicate that the MTHFR C677T polymorphism and folate deficiency interact to increase the stochastic behavior of the FOCM network, with the greatest instability observed for reactions catalyzed by serine hydroxymethyltransferase (SHMT). Furthermore, we show that de novo dTMP synthesis does not occur in the cytosol at rates sufficient for DNA replication, supporting empirical data indicating that impaired nuclear de novo dTMP synthesis results in uracil misincorporation into DNA.

Mathematical modeling of the methionine cycle and transsulfuration pathway in individuals with autism spectrum disorder

Previous research has shown a connection between metabolic abnormalities in the methionine cycle and transsulfuration pathway and autism spectrum disorder. Using clinical data from a case-control study investigating measurements of transmethylation and transsulfuration metabolites, a steady-state model of these metabolites in liver cells was developed and participant-specific parameters were identified. Comparison of mean parameter values and parameter distributions between neurotypical study participants and those on the autism spectrum revealed significant differences for four model parameters. Sensitivity analysis identified the parameter describing the rate of glutamylcysteine synthesis, the rate-limiting step in glutathione production, to be particularly important in determining steady-state metabolite concentrations. These results may provide insight into key reactions to target for potential intervention strategies relating to autism spectrum disorder.

Construction of feasible and accurate kinetic models of metabolism: A Bayesian approach

Kinetic models are essential to quantitatively understand and predict the behaviour of metabolic networks. Detailed and thermodynamically feasible kinetic models of metabolism are inherently difficult to formulate and fit. They have a large number of heterogeneous parameters, are non-linear and have complex interactions. Many powerful fitting strategies are ruled out by the intractability of the likelihood function. Here, we have developed a computational framework capable of fitting feasible and accurate kinetic models using Approximate Bayesian Computation. This framework readily supports advanced modelling features such as model selection and model-based experimental design. We illustrate this approach on the tightly-regulated mammalian methionine cycle. Sampling from the posterior distribution, the proposed framework generated thermodynamically feasible parameter samples that converged on the true values, and displayed remarkable prediction accuracy in several validation tests. Furthermore, a posteriori analysis of the parameter distributions enabled appraisal of the systems properties of the network (e.g., control structure) and key metabolic regulations. Finally, the framework was used to predict missing allosteric interactions.

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.

Using mathematical models to understand metabolism, genes, and disease

Mathematical models are a useful tool for investigating a large number of questions in metabolism, genetics, and gene–environment interactions. A model based on the underlying biology and biochemistry is a platform for in silico biological experimentation that can reveal the causal chain of events that connect variation in one quantity to variation in another. We discuss how we construct such models, how we have used them to investigate homeostatic mechanisms, gene–environment interactions, and genotype–phenotype mapping, and how they can be used in precision and personalized medicine.

Targeted metabolomics and mathematical modeling demonstrate that vitamin B-6 restriction alters one-carbon metabolism in cultured HepG2 cells

Low vitamin B-6 nutritional status is associated with increased risk for cardiovascular disease and certain cancers. Pyridoxal 5'-phosphate (PLP) serves as a coenzyme in many cellular processes, including several reactions in one-carbon (1C) metabolism and the transsulfuration pathway of homocysteine catabolism. To assess the effect of vitamin B-6 deficiency on these processes and associated pathways, we conducted quantitative analysis of 1C metabolites including tetrahydrofolate species in HepG2 cells cultured in various concentrations of pyridoxal. These results were compared with predictions of a mathematical model of 1C metabolism simulating effects of vitamin B-6 deficiency. In cells cultured in vitamin B-6-deficient medium (25 or 35 nmol/l pyridoxal), we observed >200% higher concentrations of betaine (P < 0.05) and creatinine (P < 0.05) and >60% lower concentrations of creatine (P < 0.05) and 5,10-methenyltetrahydrofolate (P < 0.05) compared with cells cultured in medium containing intermediate (65 nmol/l) or the supraphysiological 2,015 nmol/l pyridoxal. Cystathionine, cysteine, glutathione, and cysteinylglycine, which are components of the transsulfuration pathway and subsequent reactions, exhibited greater concentrations at the two lower vitamin B-6 concentrations. Partial least squares discriminant analysis showed differences in overall profiles between cells cultured in 25 and 35 nmol/l pyridoxal vs. those in 65 and 2,015 nmol/l pyridoxal. Mathematical model predictions aligned with analytically derived results. These data reveal pronounced effects of vitamin B-6 deficiency on 1C-related metabolites, including previously unexpected secondary effects on creatine. These results complement metabolomic studies in humans demonstrating extended metabolic effects of vitamin B-6 insufficiency.

A population model of folate-mediated one-carbon metabolism

Background: Previous mathematical models for hepatic and tissue one-carbon metabolism have been combined and extended to include a blood plasma compartment. We use this model to study how the concentrations of metabolites that can be measured in the plasma are related to their respective intracellular concentrations. Methods: The model consists of a set of ordinary differential equations, one for each metabolite in each compartment, and kinetic equations for metabolism and for transport between compartments. The model was validated by comparison to a variety of experimental data such as the methionine load test and variation in folate intake. We further extended this model by introducing random and systematic variation in enzyme activity. Outcomes and Conclusions: A database of 10,000 virtual individuals was generated, each with a quantitatively different one-carbon metabolism. Our population has distributions of folate and homocysteine in the plasma and tissues that are similar to those found in the NHANES data. The model reproduces many other sets of clinical data. We show that tissue and plasma folate is highly correlated, but liver and plasma folate much less so. Oxidative stress increases the plasma S-adenosylmethionine/S-adenosylhomocysteine (SAM/SAH) ratio. We show that many relationships among variables are nonlinear and in many cases we provide explanations. Sampling of subpopulations produces dramatically different apparent associations among variables. The model can be used to simulate populations with polymorphisms in genes for folate metabolism and variations in dietary input.

Mathematical modeling of folate metabolism

Folate metabolism is a complex biological process that is influenced by many variables including transporters, cofactors, and enzymes. Mathematical models provide a useful tool to evaluate this complex system and to elucidate hypotheses that would be otherwise untenable to test in vitro or in vivo. Forty years of model development and refinement along with enhancements in technology have led to systematic improvement in our biological understanding of these models. However, increased complexity does not always lead to increased understanding, and a balanced approach to modeling the system is often advantageous. This approach should address questions about sensitivity of the model to variation and incorporate genomic data. The folate model is a useful platform for investigating the effects of antifolates on the folate pathway. The utility of the model is demonstrated through interrogation of drug resistance, drug-drug interactions, drug selectivity, and drug doses and schedules. Mathematics can be used to create models with the ability to design and improve rationale therapeutic interventions.

Modeling cellular compartmentation in one-carbon metabolism

Folate-mediated one-carbon metabolism (FOCM) is associated with risk for numerous pathological states including birth defects, cancers, and chronic diseases. Although the enzymes that constitute the biological pathways have been well described and their interdependency through the shared use of folate cofactors appreciated, the biological mechanisms underlying disease etiologies remain elusive. The FOCM network is highly sensitive to nutritional status of several B-vitamins and numerous penetrant gene variants that alter network outputs, but current computational approaches do not fully capture the dynamics and stochastic noise of the system. Combining the stochastic approach with a rule-based representation will help model the intrinsic noise displayed by FOCM, address the limited flexibility of standard simulation methods for coarse-graining the FOCM-associated biochemical processes, and manage the combinatorial complexity emerging from reactions within FOCM that would otherwise be intractable.

Biodegradation kinetics and mechanism of 2,4,6-tribromophenol by Bacillus sp. GZT: a phenomenon of xenobiotic methylation during debromination

A strain Bacillus sp. GZT capable of debrominating and mineralizing 2,4,6-tribromophenol (TBP) was isolated and characterized by morphological observation, biochemical and physiological identification as well as 16S rRNA sequence analysis. Biodegradation kinetics experiments demonstrated that initial TBP concentration had a predominant effect on degradation efficiency. Within 120h, the highest TBP degradation and debromination efficiencies were up to 93.2% and 89.3%, respectively, under the optimum condition. Ten metabolic intermediates including five brominated compounds, three oxidative products and two cellular metabolites were all identified by gas chromatography-mass spectrometer, and six key intermediates were doubly validated by authentic standards. The proposed biodegradation mechanism inferred that reductive debromination as a major degradation pathway could simultaneously take place at ortho- and para-positions on TBP, while methylated debromination was also found as a minor degradation pathway during this process. Within 148h degradation, nearly one-third of 3mg/L TBP could be completely mineralized.

A genetic signature of spina bifida risk from pathway-informed comprehensive gene-variant analysis

Despite compelling epidemiological evidence that folic acid supplements reduce the frequency of neural tube defects (NTDs) in newborns, common variant association studies with folate metabolism genes have failed to explain the majority of NTD risk. The contribution of rare alleles as well as genetic interactions within the folate pathway have not been extensively studied in the context of NTDs. Thus, we sequenced the exons in 31 folate-related genes in a 480-member NTD case-control population to identify the full spectrum of allelic variation and determine whether rare alleles or obvious genetic interactions within this pathway affect NTD risk. We constructed a pathway model, predetermined independent of the data, which grouped genes into coherent sets reflecting the distinct metabolic compartments in the folate/one-carbon pathway (purine synthesis, pyrimidine synthesis, and homocysteine recycling to methionine). By integrating multiple variants based on these groupings, we uncovered two provocative, complex genetic risk signatures. Interestingly, these signatures differed by race/ethnicity: a Hispanic risk profile pointed to alterations in purine biosynthesis, whereas that in non-Hispanic whites implicated homocysteine metabolism. In contrast, parallel analyses that focused on individual alleles, or individual genes, as the units by which to assign risk revealed no compelling associations. These results suggest that the ability to layer pathway relationships onto clinical variant data can be uniquely informative for identifying genetic risk as well as for generating mechanistic hypotheses. Furthermore, the identification of ethnic-specific risk signatures for spina bifida resonated with epidemiological data suggesting that the underlying pathogenesis may differ between Hispanic and non-Hispanic groups.

A combined model of hepatic polyamine and sulfur amino acid metabolism to analyze S-adenosyl methionine availability

Many molecular details remain to be uncovered concerning the regulation of polyamine metabolism. A previous model of mammalian polyamine metabolism showed that S-adenosyl methionine availability could play a key role in polyamine homeostasis. To get a deeper insight in this prediction, we have built a combined model by integration of the previously published polyamine model and one-carbon and glutathione metabolism model, published by different research groups. The combined model is robust and it is able to achieve physiological steady-state values, as well as to reproduce the predictions of the individual models. Furthermore, a transition between two versions of our model with new regulatory factors added properly simulates the switch in methionine adenosyl transferase isozymes occurring when the liver enters in proliferative conditions. The combined model is useful to support the previous prediction on the role of S-adenosyl methionine availability in polyamine homeostasis. Furthermore, it could be easily adapted to get deeper insights on the connections of polyamines with energy metabolism.

Genomic and epigenomic instability, fragile sites, schizophrenia and autism

Increasing evidence links genomic and epigenomic instability, including multiple fragile sites regions to neuropsychiatric diseases including schizophrenia and autism. Cancer is the only other disease associated with multiple fragile site regions, and genome and epigenomic instability is a characteristic of cancer. Research on cancer is far more advanced than research on neuropsychiatric disease; hence, insight into neuropsychiatric disease may be derived from cancer research results. Towards this end, this article will review the evidence linking schizophrenia and other neuropsychiatric diseases (especially autism) to genomic and epigenomic instability, and fragile sites. The results of studies on genetic, epigenetic and environmental components of schizophrenia and autism point to the importance of the folate-methionine-transulfuration metabolic hub that is diseases also perturbed in cancer. The idea that the folate-methionine-transulfuration hub is important in neuropsychiatric is exciting because this hub present novel targets for drug development, suggests some drugs used in cancer may be useful in neuropsychiatric disease, and raises the possibility that nutrition interventions may influence the severity, presentation, or dynamics of disease.

A kinetic model describes metabolic response to perturbations and distribution of flux control in the benzenoid network of Petunia hybrida

In recent years there has been much interest in the genetic enhancement of plant metabolism; however, attempts at genetic modification are often unsuccessful due to an incomplete understanding of network dynamics and their regulatory properties. Kinetic modeling of plant metabolic networks can provide predictive information on network control and response to genetic perturbations, which allow estimation of flux at any concentration of intermediate or enzyme in the system. In this research, a kinetic model of the benzenoid network was developed to simulate whole network responses to different concentrations of supplied phenylalanine (Phe) in petunia flowers and capture flux redistributions caused by genetic manipulations. Kinetic parameters were obtained by network decomposition and non-linear least squares optimization of data from petunia flowers supplied with either 75 or 150 mm(2)H(5)-Phe. A single set of kinetic parameters simultaneously accommodated labeling and pool size data obtained for all endogenous and emitted volatiles at the two concentrations of supplied (2)H(5)-Phe. The generated kinetic model was validated using flowers from transgenic petunia plants in which benzyl CoA:benzyl alcohol/phenylethanol benzoyltransferase (BPBT) was down-regulated via RNAi. The determined in vivo kinetic parameters were used for metabolic control analysis, in which flux control coefficients were calculated for fluxes around the key branch point at Phe and revealed that phenylacetaldehyde synthase activity is the primary controlling factor for the phenylacetaldehyde branch of the benzenoid network. In contrast, control of flux through the beta-oxidative and non-beta-oxidative pathways is highly distributed.

Distribution of methionine between cells and incubation medium in suspension of rat hepatocytes

Methionine is an essential amino acid involved in many significant intracellular processes. Aberrations in methionine metabolism are associated with a number of complex pathologies. Liver plays a key role in regulation of blood methionine level. Investigation of methionine distribution between hepatocytes and medium is crucial for understanding the mechanisms of this regulation. For the first time, we analyzed the distribution of methionine between hepatocytes and incubation medium using direct measurements of methionine concentrations. Our results revealed a fast and reversible transport of methionine through the cell membrane that provides almost uniform distribution of methionine between hepatocytes and incubation medium. The steady-state ratio between intracellular and extracellular methionine concentrations was established within a few minutes. This ratio was found to be 1.06±0.38, 0.89±0.26, 0.67±0.16 and 0.82±0.06 at methionine concentrations in the medium of 64±19, 152±39, 413±55, and 1,204±104 μmol/L, respectively. The fast and uniform distribution of methionine between hepatocytes and extracellular compartments provides a possibility for effective regulation of blood methionine levels due to methionine metabolism in hepatocytes.

A theory for the arrangement of sensory organs in Drosophila

We study the arrangements of recurved bristles on the anterior wing margin of wild-type and mutant Drosophila. The epidermal or neural fate of a proneural cell depends on the concentrations of proteins of the achaete-scute complex. At puparium formation, concentrations of proteins are nearly identical in all cells of the anterior wing and each cell has the potential for neural fate. In wild-type flies, the action of regulatory networks drives the initial state to one where a bristle grows out of every fifth cell. Recent experiments have shown that the frequency of recurved bristles can be made to change by adjusting the mean concentrations of the zinc-finger transcription factor Senseless and the micro-RNA miR-9a. Specifically, mutant flies with reduced levels of miR-9a exhibit ectopic bristles, and those with lower levels of both miR-9a and Senseless show regular organization of recurved bristles, but with a lower periodicity of 4. We argue that these characteristics can be explained assuming an underlying Turing-type bifurcation whereby a periodic pattern spontaneously emerges from a uniform background. However, bristle patterns occur in a discrete array of cells, and are not mediated by diffusion. We argue that intracellular actions of transmembrane proteins such as Delta and Notch can play a role of diffusion in destabilizing the homogeneous state. In contrast to diffusion, intercellular actions can be activating or inhibiting; further, there can be lateral cross-species interactions. We introduce a phenomenological model to study bristle arrangements and make several model-independent predictions that can be tested in experiments. In our theory, miRNA-9a is one of the components of the underlying network and has no special regulatory role. The loss of periodicity in its absence is due to the transfer of the system to a bistable state.

The nutrigenetics of hyperhomocysteinemia: quantitative proteomics reveals differences in the methionine cycle enzymes of gene-induced versus diet-induced hyperhomocysteinemia

Hyperhomocysteinemia has long been associated with atherosclerosis and thrombosis and is an independent risk factor for cardiovascular disease. Its causes include both genetic and environmental factors. Although homocysteine is produced in every cell as an intermediate of the methionine cycle, the liver contributes the major portion found in circulation, and fatty liver is a common finding in homocystinuric patients. To understand the spectrum of proteins and associated pathways affected by hyperhomocysteinemia, we analyzed the mouse liver proteome of gene-induced (cystathionine beta-synthase (CBS)) and diet-induced (high methionine) hyperhomocysteinemic mice using two-dimensional difference gel electrophoresis and Ingenuity Pathway Analysis. Nine proteins were identified whose expression was significantly changed by 2-fold (p < or = 0.05) as a result of genotype, 27 proteins were changed as a result of diet, and 14 proteins were changed in response to genotype and diet. Importantly, three enzymes of the methionine cycle were up-regulated. S-Adenosylhomocysteine hydrolase increased in response to genotype and/or diet, whereas glycine N-methyltransferase and betaine-homocysteine methyltransferase only increased in response to diet. The antioxidant proteins peroxiredoxins 1 and 2 increased in wild-type mice fed the high methionine diet but not in the CBS mutants, suggesting a dysregulation in the antioxidant capacity of those animals. Furthermore, thioredoxin 1 decreased in both wild-type and CBS mutants on the diet but not in the mutants fed a control diet. Several urea cycle proteins increased in both diet groups; however, arginase 1 decreased in the CBS(+/-) mice fed the control diet. Pathway analysis identified the retinoid X receptor signaling pathway as the top ranked network associated with the CBS(+/-) genotype, whereas xenobiotic metabolism and the NRF2-mediated oxidative stress response were associated with the high methionine diet. Our results show that hyperhomocysteinemia, whether caused by a genetic mutation or diet, alters the abundance of several liver proteins involved in homocysteine/methionine metabolism, the urea cycle, and antioxidant defense.

The logic of the hepatic methionine metabolic cycle

This review describes our current understanding of the "traffic lights" that regulate sulfur flow through the methionine bionetwork in liver, which supplies two major homeostatic systems governing cellular methylation and antioxidant potential. Theoretical concepts derived from mathematical modeling of this metabolic nexus provide insights into the properties of this system, some of which seem to be paradoxical at first glance. Cellular needs supported by this network are met by use of parallel metabolic tracks that are differentially controlled by intermediates in the pathway. A major task, i.e. providing cellular methylases with the methylating substrate, S-adenosylmethionine, is met by flux through the methionine adenosyltransferase I isoform. On the other hand, a second important function, i.e., stabilization of the blood methionine concentration in the face of high dietary intake of this amino acid, is achieved by switching to an alternative isoform, methionine adenosyltransferase III, and to glycine N-methyl transferase, which together bypass the first two reactions in the methionine cycle. This regulatory strategy leads to two metabolic modes that differ in metabolite concentrations and metabolic rates almost by an order of magnitude. Switching between these modes occurs in a narrow trigger zone of methionine concentration. Complementary experimental and theoretical analyses of hepatic methionine metabolism have been richly informative and have the potential to illuminate its response to oxidative challenge, to methionine restriction and lifespan extension studies and to diseases resulting from deficiencies at specific loci in this pathway.

Individualization of thiopurine therapy: thiopurine S-methyltransferase and beyond

The metabolism of a given drug depends, not solely on a particular enzyme, but rather on a complex metabolic network. Thiopurine S-methyltransferase (TPMT) catalyzes the methylation, and thus deactivation, of 6-mercaptopurine, a thiopurine used in the treatment of acute lymphoblastic leukemia. Low TPMT activity has been associated with severe toxicity of 6-mercaptopurine. Determination of mutations in the TPMT gene before starting 6-mercaptopurine therapy constitutes a quick, simple and cost-effective strategy to individualize thiopurine dosing. However, TPMT phenotype-to-genotype correlation is not complete, indicating a need for identification of novel biomarkers. Based on our recent findings and reviewing seemingly unrelated literature reports we present a synthesis of the current understanding of factors that influence TPMT activity and consequently modulate responsiveness to thiopurine treatment. Identification and understanding of these factors is crucial for improving the efficacy and safety of acute lymphoblastic leukemia treatment.

First steps in computational systems biology: A practical session in metabolic modeling and simulation

A comprehensive understanding of biological functions requires new systemic perspectives, such as those provided by systems biology. Systems biology approaches are hypothesis-driven and involve iterative rounds of model building, prediction, experimentation, model refinement, and development. Developments in computer science are allowing for ever faster numerical simulations of mathematical models. Mathematical modeling plays an essential role in new systems biology approaches. As a complex, integrated system, metabolism is a suitable topic of study for systems biology approaches. However, up until recently, this topic has not been properly covered in biochemistry courses. This communication reports the development and implementation of a practical lesson plan on metabolic modeling and simulation.

Role of sulfur amino acids in controlling nutrient metabolism and cell functions: implications for nutrition

Protein synthesis is affected when an insufficient level of sulfur amino acids is available. This defect may originate from dietary amino acid deficiency and/or excessive amino acid utilisation for other purposes such as the synthesis of glutathione and acute-phase proteins during catabolic stress. Sulfur amino acids are recognised to exert other significant functions since they are precursors of essential molecules, are involved in the methylation process, participate in the control of oxidative status, and may act as mediators affecting metabolism and cell functions. Despite this increased understanding of the role of sulfur amino acids, many questions still remain unanswered due to the complexity of the mechanisms involved. Moreover, surprising effects of dietary sulfur amino acids have been reported, with the development of disorders in cases of both deficiency and excess. These findings indicate the importance of defining adequate levels of intake and providing a rationale for nutritional advice. The aim of the present review is to provide an overview on the roles of sulfur amino acids as regulators of nutrient metabolism and cell functions, with emphasis placed on the implications for nutrition.

A domino effect in antifolate drug action in Escherichia coli

Mass spectrometry technologies for measurement of cellular metabolism are opening new avenues to explore drug activity. Trimethoprim is an antibiotic that inhibits bacterial dihydrofolate reductase (DHFR). Kinetic flux profiling with (15)N-labeled ammonia in Escherichia coli reveals that trimethoprim leads to blockade not only of DHFR but also of another critical enzyme of folate metabolism: folylpoly-gamma-glutamate synthetase (FP-gamma-GS). Inhibition of FP-gamma-GS is not directly due to trimethoprim. Instead, it arises from accumulation of DHFR's substrate dihydrofolate, which we show is a potent FP-gamma-GS inhibitor. Thus, owing to the inherent connectivity of the metabolic network, falling DHFR activity leads to falling FP-gamma-GS activity in a domino-like cascade. This cascade results in complex folate dynamics, and its incorporation in a computational model of folate metabolism recapitulates the dynamics observed experimentally. These results highlight the potential for quantitative analysis of cellular metabolism to reveal mechanisms of drug action.