Quantifying the propagation of parametric uncertainty on flux balance analysis.

Flux balance analysis (FBA) and associated techniques operating on stoichiometric genome-scale metabolic models play a central role in quantifying metabolic flows and constraining feasible phenotypes. At the heart of these methods lie two important assumptions: (i) the biomass precursors and energy requirements neither change in response to growth conditions nor environmental/genetic perturbations, and (ii) metabolite production and consumption rates are equal at all times (i.e., steady-state). Despite the stringency of these two assumptions, FBA has been shown to be surprisingly robust at predicting cellular phenotypes. In this paper, we formally assess the impact of these two assumptions on FBA results by quantifying how uncertainty in biomass reaction coefficients, and departures from steady-state due to temporal fluctuations could propagate to FBA results. In the first case, conditional sampling of parameter space is required to re-weigh the biomass reaction so as the molecular weight remains equal to 1 g mmol-1, and in the second case, metabolite (and elemental) pool conservation must be imposed under temporally varying conditions. Results confirm the importance of enforcing the aforementioned constraints and explain the robustness of FBA biomass yield predictions.

Dissecting the metabolic reprogramming of maize root under nitrogen-deficient stress conditions.

The growth and development of maize (Zea mays L.) largely depends on its nutrient uptake through the root. Hence, studying its growth, response, and associated metabolic reprogramming to stress conditions is becoming an important research direction. A genome-scale metabolic model (GSM) for the maize root was developed to study its metabolic reprogramming under nitrogen stress conditions. The model was reconstructed based on the available information from KEGG, UniProt, and MaizeCyc. Transcriptomics data derived from the roots of hydroponically grown maize plants were used to incorporate regulatory constraints in the model and simulate nitrogen-non-limiting (N+) and nitrogen-deficient (N-) condition. Model-predicted flux-sum variability analysis achieved 70% accuracy compared with the experimental change of metabolite levels. In addition to predicting important metabolic reprogramming in central carbon, fatty acid, amino acid, and other secondary metabolism, maize root GSM predicted several metabolites (l-methionine, l-asparagine, l-lysine, cholesterol, and l-pipecolate) playing a regulatory role in the root biomass growth. Furthermore, this study revealed eight phosphatidylcholine and phosphatidylglycerol metabolites which, even though not coupled with biomass production, played a key role in the increased biomass production under N-deficient conditions. Overall, the omics-integrated GSM provides a promising tool to facilitate stress condition analysis for maize root and engineer better stress-tolerant maize genotypes.

Elucidation of trophic interactions in an unusual single-cell nitrogen-fixing symbiosis using metabolic modeling.

Marine nitrogen-fixing microorganisms are an important source of fixed nitrogen in oceanic ecosystems. The colonial cyanobacterium Trichodesmium and diatom symbionts were thought to be the primary contributors to oceanic N2 fixation until the discovery of the unusual uncultivated symbiotic cyanobacterium UCYN-A (Candidatus Atelocyanobacterium thalassa). UCYN-A has atypical metabolic characteristics lacking the oxygen-evolving photosystem II, the tricarboxylic acid cycle, the carbon-fixation enzyme RuBisCo and de novo biosynthetic pathways for a number of amino acids and nucleotides. Therefore, it is obligately symbiotic with its single-celled haptophyte algal host. UCYN-A receives fixed carbon from its host and returns fixed nitrogen, but further insights into this symbiosis are precluded by both UCYN-A and its host being uncultured. In order to investigate how this syntrophy is coordinated, we reconstructed bottom-up genome-scale metabolic models of UCYN-A and its algal partner to explore possible trophic scenarios, focusing on nitrogen fixation and biomass synthesis. Since both partners are uncultivated and only the genome sequence of UCYN-A is available, we used the phylogenetically related Chrysochromulina tobin as a proxy for the host. Through the use of flux balance analysis (FBA), we determined the minimal set of metabolites and biochemical functions that must be shared between the two organisms to ensure viability and growth. We quantitatively investigated the metabolic characteristics that facilitate daytime N2 fixation in UCYN-A and possible oxygen-scavenging mechanisms needed to create an anaerobic environment to allow nitrogenase to function. This is the first application of an FBA framework to examine the tight metabolic coupling between uncultivated microbes in marine symbiotic communities and provides a roadmap for future efforts focusing on such specialized systems.

SNPeffect: identifying functional roles of SNPs using metabolic networks.

Genetic sources of phenotypic variation have been a focus of plant studies aimed at improving agricultural yield and understanding adaptive processes. Genome-wide association studies identify the genetic background behind a trait by examining associations between phenotypes and single-nucleotide polymorphisms (SNPs). Although such studies are common, biological interpretation of the results remains a challenge; especially due to the confounding nature of population structure and the systematic biases thus introduced. Here, we propose a complementary analysis (SNPeffect) that offers putative genotype-to-phenotype mechanistic interpretations by integrating biochemical knowledge encoded in metabolic models. SNPeffect is used to explain differential growth rate and metabolite accumulation in A. thaliana and P. trichocarpa accessions as the outcome of SNPs in enzyme-coding genes. To this end, we also constructed a genome-scale metabolic model for Populus trichocarpa, the first for a perennial woody tree. As expected, our results indicate that growth is a complex polygenic trait governed by carbon and energy partitioning. The predicted set of functional SNPs in both species are associated with experimentally characterized growth-determining genes and also suggest putative ones. Functional SNPs were found in pathways such as amino acid metabolism, nucleotide biosynthesis, and cellulose and lignin biosynthesis, in line with breeding strategies that target pathways governing carbon and energy partition.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms.

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

A diurnal flux balance model of Synechocystis sp. PCC 6803 metabolism.

Phototrophic organisms such as cyanobacteria utilize the sun's energy to convert atmospheric carbon dioxide into organic carbon, resulting in diurnal variations in the cell's metabolism. Flux balance analysis is a widely accepted constraint-based optimization tool for analyzing growth and metabolism, but it is generally used in a time-invariant manner with no provisions for sequestering different biomass components at different time periods. Here we present CycleSyn, a periodic model of Synechocystis sp. PCC 6803 metabolism that spans a 12-hr light/12-hr dark cycle by segmenting it into 12 Time Point Models (TPMs) with a uniform duration of two hours. The developed framework allows for the flow of metabolites across TPMs while inventorying metabolite levels and only allowing for the utilization of currently or previously produced compounds. The 12 TPMs allow for the incorporation of time-dependent constraints that capture the cyclic nature of cellular processes. Imposing bounds on reactions informed by temporally-segmented transcriptomic data enables simulation of phototrophic growth as a single linear programming (LP) problem. The solution provides the time varying reaction fluxes over a 24-hour cycle and the accumulation/consumption of metabolites. The diurnal rhythm of metabolic gene expression driven by the circadian clock and its metabolic consequences is explored. Predicted flux and metabolite pools are in line with published studies regarding the temporal organization of phototrophic growth in Synechocystis PCC 6803 paving the way for constructing time-resolved genome-scale models (GSMs) for organisms with a circadian clock. In addition, the metabolic reorganization that would be required to enable Synechocystis PCC 6803 to temporally separate photosynthesis from oxygen-sensitive nitrogen fixation is also explored using the developed model formalism.

Chlorin e6 Conjugated Methoxy-Poly(Ethylene Glycol)-Poly(D,L-Lactide) Glutathione Sensitive Micelles for Photodynamic Therapy.

PURPOSE: In this study, we developed a polymeric micellar system for glutathione-mediated intracellular delivery of a photosensitizer, chlorin e6 (Ce6) by synthesizing an amphiphilic polymer, methoxy-poly(ethylene glycol)-poly(D,L-lactide)-disulfide-Ce6 (mPEG-PLA-S-S-Ce6), which self-assembled in aqueous environment to form micelles. METHODS: The polymer-drug conjugate was characterized by NMR. The singlet oxygen (2O1) generation and in vitro release of Ce6 micelles were evaluated. Further, glutathione-mediated intracellular drug delivery was assessed in human alveolar adenocarcinoma cells (A549), mouse mammary carcinoma cells (4 T1) and 3D A549 spheroids. RESULTS: The micellar system protected Ce6 from aggregation leading to improved 2O1 generation compared to free Ce6. Due to the availability of glutathione, the disulfide bonds in the micelles were cleaved resulting in rapid release of Ce6 evident from the in vitro study. The Ce6 micelles displayed quicker drug release in presence of glutathione monoester (GSH-OEt) pre-treated A549 and 4 T1 cells compared to without pre-treated cells. In vitro phototoxicity of micelles displayed enhanced toxicity in 10 mM GSH-OEt pre-treated A549 and 4 T1 cells compared to untreated cells. As anticipated, Ce6 micelles showed lower drug release in presence of 0.1 mM of buthionine sulfoximine (BSO) pretreated A549 and 4 T1 cells exhibiting lower phototoxicity. Further, A549 3D spheroids treated with Ce6 micelles showed significant inhibition in growth, enhanced phototoxicity, and cellular apoptosis in comparison to free Ce6. CONCLUSION: The above results showed that the developed strategy could be effective in improving the PDT efficacy of Ce6, and the developed polymeric micellar system could be utilized as a PDT regimen for cancer.

Blood pressure variation with gestational age and birth weight in Indian newborn.

Prospective observational study was conducted in a tertiary care hospital of India over 8 months to measure blood pressure (BP) in healthy term and preterm neonates using oscillometric method and explore the associations with gestational age and birth weight. Consecutive BP measurements were taken by standard oscillometric method on 1617 neonates on day 4, 7 and 14 of life. Mean birth weight was 2.7 ± 0.46 kg, and mean gestational age was 38.2 ± 2.12 weeks. The mean arterial pressure (MAP) on day 4, 7 and 14 were 59.3 ± 7.33, 63.2 ± 6.55 and 66.4 ± 6.13 mmHg, respectively. Larger and mature newborns had significantly higher BP than those who were smaller and premature. Birth weight more strongly correlated with MAP than gestational age. Predictive equations linking MAP with gestational age and birth weight were deduced, which can be used for judicious fluid inotrope management.

Genome streamlining to improve performance of a fast-growing cyanobacterium Synechococcus elongatus UTEX 2973.

Cyanobacteria are photosynthetic organisms that have garnered significant recognition as potential hosts for sustainable bioproduction. However, their complex regulatory networks pose significant challenges to major metabolic engineering efforts, thereby limiting their feasibility as production hosts. Genome streamlining has been demonstrated to be a successful approach for improving productivity and fitness in heterotrophs but is yet to be explored to its full potential in phototrophs. Here, we present the systematic reduction of the genome of the cyanobacterium exhibiting the fastest exponential growth, Synechococcus elongatus UTEX 2973. This work, the first of its kind in a photoautotroph, involved an iterative process using state-of-the-art genome-editing technology guided by experimental analysis and computational tools. CRISPR-Cas3 enabled large, progressive deletions of predicted dispensable regions and aided in the identification of essential genes. The large deletions were combined to obtain a strain with 55-kb genome reduction. The strains with streamlined genome showed improvement in growth (up to 23%) and productivity (by 22.7%) as compared to the wild type (WT). This streamlining strategy not only has the potential to develop cyanobacterial strains with improved growth and productivity traits but can also facilitate a better understanding of their genome-to-phenome relationships.IMPORTANCEGenome streamlining is an evolutionary strategy used by natural living systems to dispense unnecessary genes from their genome as a mechanism to adapt and evolve. While this strategy has been successfully borrowed to develop synthetic heterotrophic microbial systems with desired phenotype, it has not been extensively explored in photoautotrophs. Genome streamlining strategy incorporates both computational predictions to identify the dispensable regions and experimental validation using genome-editing tool, and in this study, we have employed a modified strategy with the goal to minimize the genome size to an extent that allows optimal cellular fitness under specified conditions. Our strategy has explored a novel genome-editing tool in photoautotrophs, which, unlike other existing tools, enables large, spontaneous optimal deletions from the genome. Our findings demonstrate the effectiveness of this modified strategy in obtaining strains with streamlined genome, exhibiting improved fitness and productivity.

Platforms for High-Throughput Screening and Force Measurements on Fungi and Oomycetes.

Pathogenic fungi and oomycetes give rise to a significant number of animal and plant diseases. While the spread of these pathogenic microorganisms is increasing globally, emerging resistance to antifungal drugs is making associated diseases more difficult to treat. High-throughput screening (HTS) and new developments in lab-on-a-chip (LOC) platforms promise to aid the discovery of urgently required new control strategies and anti-fungal/oomycete drugs. In this review, we summarize existing HTS and emergent LOC approaches in the context of infection strategies and invasive growth exhibited by these microorganisms. To aid this, we introduce key biological aspects and review existing HTS platforms based on both conventional and LOC techniques. We then provide an in-depth discussion of more specialized LOC platforms for force measurements on hyphae and to study electro- and chemotaxis in spores, approaches which have the potential to aid the discovery of alternative drug targets on future HTS platforms. Finally, we conclude with a brief discussion of the technical developments required to improve the uptake of these platforms into the general laboratory environment.

A Genome-Scale Metabolic Model of Anabaena 33047 to Guide Genetic Modifications to Overproduce Nylon Monomers.

Nitrogen fixing-cyanobacteria can significantly improve the economic feasibility of cyanobacterial production processes by eliminating the requirement for reduced nitrogen. Anabaena sp. ATCC 33047 is a marine, heterocyst forming, nitrogen fixing cyanobacteria with a very short doubling time of 3.8 h. We developed a comprehensive genome-scale metabolic (GSM) model, iAnC892, for this organism using annotations and content obtained from multiple databases. iAnC892 describes both the vegetative and heterocyst cell types found in the filaments of Anabaena sp. ATCC 33047. iAnC892 includes 953 unique reactions and accounts for the annotation of 892 genes. Comparison of iAnC892 reaction content with the GSM of Anabaena sp. PCC 7120 revealed that there are 109 reactions including uptake hydrogenase, pyruvate decarboxylase, and pyruvate-formate lyase unique to iAnC892. iAnC892 enabled the analysis of energy production pathways in the heterocyst by allowing the cell specific deactivation of light dependent electron transport chain and glucose-6-phosphate metabolizing pathways. The analysis revealed the importance of light dependent electron transport in generating ATP and NADPH at the required ratio for optimal N2 fixation. When used alongside the strain design algorithm, OptForce, iAnC892 recapitulated several of the experimentally successful genetic intervention strategies that over produced valerolactam and caprolactam precursors.

The within-host fitness of HIV-1 increases with age in ART-naïve HIV-1 subtype C infected children.

As the immune system develops with age, children combat infections better. HIV-1, however, targets an activated immune system, potentially rendering children increasingly permissive to HIV-1 infection as they grow. How HIV-1 fitness changes with age in children is unknown. Here, we estimated the within-host basic reproductive ratio, R0, a marker of viral fitness, in HIV-1 subtype C-infected children in India, aged between 84 days and 17 years. We measured serial viral load and CD4 T cell counts in 171 children who initiated first-line ART. For 25 children, regular and frequent measurements provided adequate data points for analysis using a mathematical model of viral dynamics to estimate R0. For the rest, we used CD4 counts for approximate estimation of R0. The viral load decline during therapy was biphasic. The mean lifespans of productively and long-lived infected cells were 1.4 and 27.8 days, respectively. The mean R0 was 1.5 in children aged < 5 years, increased with age, and approached 6.0 at 18 years, close to 5.8 estimated previously for adults. The tolerogenic immune environment thus compromises HIV-1 fitness in young children. Early treatment initiation, when the R0 is small, will likely improve viral control, in addition to suppressing the latent reservoir.

Bacterial colonization reprograms the neonatal gut metabolome.

Initial microbial colonization and later succession in the gut of human infants are linked to health and disease later in life. The timing of the appearance of the first gut microbiome, and the consequences for the early life metabolome, are just starting to be defined. Here, we evaluated the gut microbiome, proteome and metabolome in 88 African-American newborns using faecal samples collected in the first few days of life. Gut bacteria became detectable using molecular methods by 16 h after birth. Detailed analysis of the three most common species, Escherichia coli, Enterococcus faecalis and Bacteroides vulgatus, did not suggest a genomic signature for neonatal gut colonization. The appearance of bacteria was associated with reduced abundance of approximately 50 human proteins, decreased levels of free amino acids and an increase in products of bacterial fermentation, including acetate and succinate. Using flux balance modelling and in vitro experiments, we provide evidence that fermentation of amino acids provides a mechanism for the initial growth of E. coli, the most common early colonizer, under anaerobic conditions. These results provide a deep characterization of the first microbes in the human gut and show how the biochemical environment is altered by their appearance.

Enhanced clearance of HIV-1-infected cells by broadly neutralizing antibodies against HIV-1 in vivo.

Antiretroviral drugs and antibodies limit HIV-1 infection by interfering with the viral life cycle. In addition, antibodies also have the potential to guide host immune effector cells to kill HIV-1-infected cells. Examination of the kinetics of HIV-1 suppression in infected individuals by passively administered 3BNC117, a broadly neutralizing antibody, suggested that the effects of the antibody are not limited to free viral clearance and blocking new infection but also include acceleration of infected cell clearance. Consistent with these observations, we find that broadly neutralizing antibodies can target CD4(+) T cells infected with patient viruses and can decrease their in vivo half-lives by a mechanism that requires Fcγ receptor engagement in a humanized mouse model. The results indicate that passive immunotherapy can accelerate elimination of HIV-1-infected cells.

Normative blood pressure data for Indian neonates.

OBJECTIVE: To establish the normative blood pressure (BP) values in healthy Indian neonates using oscillometric method, and to develop BP percentile charts. DESIGN: Prospective observational study. SETTING: Neonatal unit of a teaching hospital in Eastern India. PARTICIPANTS: 1617 hemodynamically stable inborn neonates without birth asphyxia, major congenital anomaly, maternal complications (e.g. preeclampsia, hypertension, diabetes) or critical neonatal illness. PROCEDURE: Quite state measurements of systolic BP (SBP), diastolic BP (DBP) and mean arterial pressure (MAP) were recorded by oscillometric method on day 4, 7 and 14 of postnatal life. The averages of three readings at 2-minute intervals were used. RESULTS: Percentile charts (providing 5th, 10th, 25th, 50th, 75th, 95th, and 99th percentile values) have been developed. SBP, DBP and MAP showed a steady rise from day 4 to day 14, and were comparable between males and females, but were significantly lower in preterms than in term neonates. CONCLUSIONS: Normative neonatal BP data along with gestational age-wise percentile charts shall be of help for decision-making and planning for sick newborns.