Predicting drug-protein interactions by preserving the graph information of multi source data.

Examining potential drug-target interactions (DTIs) is a pivotal component of drug discovery and repurposing. Recently, there has been a significant rise in the use of computational techniques to predict DTIs. Nevertheless, previous investigations have predominantly concentrated on assessing either the connections between nodes or the consistency of the network's topological structure in isolation. Such one-sided approaches could severely hinder the accuracy of DTI predictions. In this study, we propose a novel method called TTGCN, which combines heterogeneous graph convolutional neural networks (GCN) and graph attention networks (GAT) to address the task of DTI prediction. TTGCN employs a two-tiered feature learning strategy, utilizing GAT and residual GCN (R-GCN) to extract drug and target embeddings from the diverse network, respectively. These drug and target embeddings are then fused through a mean-pooling layer. Finally, we employ an inductive matrix completion technique to forecast DTIs while preserving the network's node connectivity and topological structure. Our approach demonstrates superior performance in terms of area under the curve and area under the precision-recall curve in experimental comparisons, highlighting its significant advantages in predicting DTIs. Furthermore, case studies provide additional evidence of its ability to identify potential DTIs.

CeCaFLUX: the first web server for standardized and visual instationary 13C metabolic flux analysis.

SUMMARY: The number of instationary 13C-metabolic flux (INST-MFA) studies grows every year, making it more important than ever to ensure the clarity, standardization and reproducibility of each study. We proposed CeCaFLUX, the first user-friendly web server that derives metabolic flux distribution from instationary 13C-labeled data. Flux optimization and statistical analysis are achieved through an evolutionary optimization in a parallel manner. It can visualize the flux optimizing process in real-time and the ultimate flux outcome. It will also function as a database to enhance the consistency and to facilitate sharing of flux studies. AVAILABILITY AND IMPLEMENTATION: CeCaFLUX is freely available at https://www.cecaflux.net, the source code can be downloaded at https://github.com/zhzhd82/CeCaFLUX. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CeCaFDB: a curated database for the documentation, visualization and comparative analysis of central carbon metabolic flux distributions explored by 13C-fluxomics.

The Central Carbon Metabolic Flux Database (CeCaFDB, available at http://www.cecafdb.org) is a manually curated, multipurpose and open-access database for the documentation, visualization and comparative analysis of the quantitative flux results of central carbon metabolism among microbes and animal cells. It encompasses records for more than 500 flux distributions among 36 organisms and includes information regarding the genotype, culture medium, growth conditions and other specific information gathered from hundreds of journal articles. In addition to its comprehensive literature-derived data, the CeCaFDB supports a common text search function among the data and interactive visualization of the curated flux distributions with compartmentation information based on the Cytoscape Web API, which facilitates data interpretation. The CeCaFDB offers four modules to calculate a similarity score or to perform an alignment between the flux distributions. One of the modules was built using an inter programming algorithm for flux distribution alignment that was specifically designed for this study. Based on these modules, the CeCaFDB also supports an extensive flux distribution comparison function among the curated data. The CeCaFDB is strenuously designed to address the broad demands of biochemists, metabolic engineers, systems biologists and members of the -omics community.

Infection rate of Eperythrozoon spp. in Chinese population: a systematic review and meta-analysis since the first Chinese case reported in 1991.

BACKGROUND: Eperythrozoonosis is an important animal health problem worldwide, it not only has a major impact on the economic viability, but also makes a significant impact on public health issues. The present systemic review intends to collate all relevant published data to assess the burden of Eperythrozoon infection in Chinese population and discuss the implications of these findings for public health policy. METHODS: A meta-analysis was conducted to review the published studies that reported Eperythrozoon spp. in Chinese population. Inclusion criteria comprised of the use of microscopic venous blood smear examination for Eperythrozoon detection and a detailed description of sampling techniques. RESULTS: Twenty-four cross-sectional studies with 52,433 participants and 14,951 positive cases, within the range of China mainland, were included in the present analysis. The infection rate of Eperythrozoon varied from 0 to 97.29% with geographical and seasonal variations, people with mild infection intensity contributed the major part (68.93%). The infection rates were highest in the children and adolescents group, significantly increased risk of Eperythrozoon infection was found among herdsmen. CONCLUSIONS: The current study raises awareness about the human eperythrozoonosis in China, which is a newly emerging zoonosis. The majority of Eperythrozoon infection intensity was asymptomatic mild infection. The infection rate of Eperythrozoon in Chinese population varied by geographical region, season, age and occupation. These factors need to be considered when conducting health education campaigns and comparing the surveillance results from different studies.

13C metabolic flux analysis: Classification and characterization from the perspective of mathematical modeling and application in physiological research of neural cell.

13C metabolic flux analysis (13C-MFA) has emerged as a forceful tool for quantifying in vivo metabolic pathway activity of different biological systems. This technology plays an important role in understanding intracellular metabolism and revealing patho-physiology mechanism. Recently, it has evolved into a method family with great diversity in experiments, analytics, and mathematics. In this review, we classify and characterize the various branch of 13C-MFA from a unified perspective of mathematical modeling. By linking different parts in the model to each step of its workflow, the specific technologies of 13C-MFA are put into discussion, including the isotope labeling model (ILM), isotope pattern measuring technique, optimization algorithm and statistical method. Its application in physiological research in neural cell has also been reviewed.

Parallel isotope differential modeling for instationary 13C fluxomics at the genome scale.

BACKGROUND: A precise map of the metabolic fluxome, the closest surrogate to the physiological phenotype, is becoming progressively more important in the metabolic engineering of photosynthetic organisms for biofuel and biomass production. For photosynthetic organisms, the state-of-the-art method for this purpose is instationary 13C fluxomics, which has arisen as a sibling of transcriptomics or proteomics. Instationary 13C data processing requires solving high-dimensional nonlinear differential equations and leads to large computational and time costs when its scope is expanded to a genome-scale metabolic network. RESULT: Here, we present a parallelized method to model instationary 13C labeling data. The elementary metabolite unit (EMU) framework is reorganized to allow treating individual mass isotopomers and breaking up of their networks into strongly connected components (SCCs). A variable domain parallel algorithm is introduced to process ordinary differential equations in a parallel way. 15-fold acceleration is achieved for constant-step-size modeling and ~ fivefold acceleration for adaptive-step-size modeling. CONCLUSION: This algorithm is universally applicable to isotope granules such as EMUs and cumomers and can substantially accelerate instationary 13C fluxomics modeling. It thus has great potential to be widely adopted in any instationary 13C fluxomics modeling.

Three-Dimensional Graphene Enhances Neural Stem Cell Proliferation Through Metabolic Regulation.

Graphene consists of two-dimensional sp2-bonded carbon sheets, a single or a few layers thick, which has attracted considerable interest in recent years due to its good conductivity and biocompatibility. Three-dimensional graphene foam (3DG) has been demonstrated to be a robust scaffold for culturing neural stem cells (NSCs) in vitro that not only supports NSCs growth, but also maintains cells in a more active proliferative state than 2D graphene films and ordinary glass. In addition, 3DG can enhance NSCs differentiation into astrocytes and especially neurons. However, the underlying mechanisms behind 3DG's effects are still poorly understood. Metabolism is the fundamental characteristic of life and provides substances for building and powering the cell. Metabolic activity is tightly tied with the proliferation, differentiation, and self-renewal of stem cells. This study focused on the metabolic reconfiguration of stem cells induced by culturing on 3DG. This study established the correlation between metabolic reconfiguration metabolomics with NSCs cell proliferation rate on different scaffold. Several metabolic processes have been uncovered in association with the proliferation change of NSCs. Especially, culturing on 3DG triggered pathways that increased amino acid incorporation and enhanced glucose metabolism. These data suggested a potential association between graphene and pathways involved in Parkinson's disease. Our work provides a very useful starting point for further studies of NSC fate determination on 3DG.

Improving metabolic flux estimation via evolutionary optimization for convex solution space.

MOTIVATION: Flux estimation by using (13) C-labeling pattern information of metabolites is currently the only method that can give accurate, detailed quantification of all intracellular fluxes in the central metabolism of a microorganism. In essence, it corresponds to a constrained optimization problem which minimizes a weighted distance between measured and simulated results. Characteristics, such as existence of multiple local minima, non-linear and non-differentiable make this problem a special difficulty. RESULTS: In the present work, we propose an evolutionary-based global optimization algorithm taking advantage of the convex feature of the problem's solution space. Based on the characteristics of convex spaces, specialized initial population and evolutionary operators are designed to solve (13)C-based metabolic flux estimation problem robustly and efficiently. The algorithm was applied to estimate the central metabolic fluxes in Escherichia coli and compared with conventional optimization technique. Experimental results illustrated that our algorithm is capable of achieving fast convergence to good near-optima and maintaining the robust nature of evolutionary algorithms at the same time. AVAILABILITY: Available from the authors upon request.

A genome-scale metabolic network alignment method within a hypergraph-based framework using a rotational tensor-vector product.

Biological network alignment aims to discover important similarities and differences and thus find a mapping between topological and/or functional components of different biological molecular networks. Then, the mapped components can be considered to correspond to both their places in the network topology and their biological attributes. Development and evolution of biological network alignment methods has been accelerated by the rapidly increasing availability of such biological networks, yielding a repertoire of tens of methods based upon graph theory. However, most biological processes, especially the metabolic reactions, are more sophisticated than simple pairwise interactions and contain three or more participating components. Such multi-lateral relations are not captured by graphs, and computational methods to overcome this limitation are currently lacking. This paper introduces hypergraphs and association hypergraphs to describe metabolic networks and their potential alignments, respectively. Within this framework, metabolic networks are aligned by identifying the maximal Z-eigenvalue of a symmetric tensor. A shifted higher-order power method was utilized to identify a solution. A rotational strategy has been introduced to accelerate the tensor-vector product by 250-fold on average and reduce the storage cost by up to 1,000-fold. The algorithm was implemented on a spark-based distributed computation cluster to significantly increase the convergence rate further by 50- to 80-fold. The parameters have been explored to understand their impact on alignment accuracy and speed. In particular, the influence of initial value selection on the stationary point has been simulated to ensure an accurate approximation of the global optimum. This framework was demonstrated by alignments among the genome-wide metabolic networks of Escherichia coli MG-1655 and Halophilic archaeon DL31. To our knowledge, this is the first genome-wide metabolic network alignment at both the metabolite level and the enzyme level. These results demonstrate that it can supply quite a few valuable insights into metabolic networks. First, this method can access the driving force of organic reactions through the chemical evolution of metabolic network. Second, this method can incorporate the chemical information of enzymes and structural changes of compounds to offer new way defining reaction class and module, such as those in KEGG. Third, as a vertex-focused treatment, this method can supply novel structural and functional annotation for ill-defined molecules. The related source code is available on request.

The Relationship between Functional and Evolutionary Correlations of Enzyme Reactions in Metabolic Network Evolution

Abstract The evolution of species is inevitably accompanied by the evolution of metabolic networks to adapt to different environments. The metabolic networks of different species were collected from the Kyoto Encyclopedia of Genes and Genomes (KEGG) website, and some enzyme reactions with the highest occurrence frequency in all species were found and are reported in this paper. The correlation coefficients of whether the enzyme reactions appear in all species were calculated, and the corresponding evolutionary correlation connection networks were calculated according to different correlation coefficient thresholds. These studies show that, as the evolutionary correlation of enzyme reactions increases, the weighted average of the mean functional concentration ratios of the enzyme reactions also increases, indicating that the functional concentration ratio of enzyme reactions has a certain correlation with the evolutionary correlation. The work presented in this paper enhances our understanding of the characteristics and general rules of metabolic network evolution.

Elementary Flux Mode Analysis Revealed Cyclization Pathway as a Powerful Way for NADPH Regeneration of Central Carbon Metabolism.

NADPH regeneration capacity is attracting growing research attention due to its important role in resisting oxidative stress. Besides, NADPH availability has been regarded as a limiting factor in production of industrially valuable compounds. The central carbon metabolism carries the carbon skeleton flux supporting the operation of NADPH-regenerating enzyme and offers flexibility in coping with NADPH demand for varied intracellular environment. To acquire an insightful understanding of its NADPH regeneration capacity, the elementary mode method was employed to compute all elementary flux modes (EFMs) of a network representative of central carbon metabolism. Based on the metabolic flux distributions of these modes, a cluster analysis of EFMs with high NADPH regeneration rate was conducted using the self-organizing map clustering algorithm. The clustering results were used to study the relationship between the flux of total NADPH regeneration and the flux in each NADPH producing enzyme. The results identified several reaction combinations supporting high NADPH regeneration, which are proven to be feasible in cells via thermodynamic analysis and coincident with a great deal of previous experimental report. Meanwhile, the reaction combinations showed some common characteristics: there were one or two decarboxylation oxidation reactions in the combinations that produced NADPH and the combination constitution included certain gluconeogenesis pathways. These findings suggested cyclization pathways as a powerful way for NADPH regeneration capacity of bacterial central carbon metabolism.

Metabolic flux ratio analysis and multi-objective optimization revealed a globally conserved and coordinated metabolic response of E. coli to paraquat-induced oxidative stress.

The ability of a microorganism to adapt to changes in the environment, such as in nutrient or oxygen availability, is essential for its competitive fitness and survival. The cellular objective and the strategy of the metabolic response to an extreme environment are therefore of tremendous interest and, thus, have been increasingly explored. However, the cellular objective of the complex regulatory structure of the metabolic changes has not yet been fully elucidated and more details regarding the quantitative behaviour of the metabolic flux redistribution are required to understand the systems-wide biological significance of this response. In this study, the intracellular metabolic flux ratios involved in the central carbon metabolism were determined by fractional (13)C-labeling and metabolic flux ratio analysis (MetaFoR) of the wild-type E. coli strain JM101 at an oxidative environment in a chemostat. We observed a significant increase in the flux through phosphoenolpyruvate carboxykinase (PEPCK), phosphoenolpyruvate carboxylase (PEPC), malic enzyme (MEZ) and serine hydroxymethyltransferase (SHMT). We applied an ε-constraint based multi-objective optimization to investigate the trade-off relationships between the biomass yield and the generation of reductive power using the in silico iJR904 genome-scale model of E. coli K-12. The theoretical metabolic redistribution supports that the trans-hydrogenase pathway should not play a direct role in the defence mounted by E. coli against oxidative stress. The agreement between the measured ratio and the theoretical redistribution established the significance of NADPH synthesis as the goal of the metabolic reprogramming that occurs in response to oxidative stress. Our work presents a framework that combines metabolic flux ratio analysis and multi-objective optimization to investigate the metabolic trade-offs that occur under varied environmental conditions. Our results led to the proposal that the metabolic response of E. coli to paraquat-induced oxidative stress is globally conserved and coordinated.

A systematic investigation of Escherichia coli central carbon metabolism in response to superoxide stress.

BACKGROUND: The cellular responses of bacteria to superoxide stress can be used to model adaptation to severe environmental changes. Superoxide stress promotes the excessive production of reactive oxygen species (ROS) that have detrimental effects on cell metabolic and other physiological activities. To antagonize such effects, the cell needs to regulate a range of metabolic reactions in a coordinated way, so that coherent metabolic responses are generated by the cellular metabolic reaction network as a whole. In the present study, we have used a quantitative metabolic flux analysis approach, together with measurement of gene expression and activity of key enzymes, to investigate changes in central carbon metabolism that occur in Escherichia coli in response to paraquat-induced superoxide stress. The cellular regulatory mechanisms involved in the observed global flux changes are discussed. RESULTS: Flux analysis based on nuclear magnetic resonance (NMR) and mass spectroscopy (MS) measurements and computation provided quantitative results on the metabolic fluxes redistribution of the E. coli central carbon network under paraquat-induced oxidative stress. The metabolic fluxes of the glycolytic pathway were redirected to the pentose phosphate pathway (PP pathway). The production of acetate increased significantly, the fluxes associated with the TCA cycle decreased, and the fluxes in the glyoxylate shunt increased in response to oxidative stress. These global flux changes resulted in an increased ratio of NADPH:NADH and in the accumulation of α-ketoglutarate. CONCLUSIONS: Metabolic flux analysis provided a quantitative and global picture of responses of the E. coli central carbon metabolic network to oxidative stress. Systematic adjustments of cellular physiological state clearly occurred in response to changes in metabolic fluxes induced by oxidative stress. Quantitative flux analysis therefore could reveal the physiological state of the cell at the systems level and is a useful complement to molecular systems approaches, such as proteomics and transcription analyses.

Increasing the accuracy of mass isotopomer analysis through calibration curves constructed using biologically synthesized compounds.

Mass isotopomer analysis is an important technique to measure the production and flow of metabolites in living cells, tissues, and organisms. This technique depends on accurate quantifications of different mass isotopomers using mass spectrometry. Constructing calibration curves using standard samples is the most universal approach to convert raw mass spectrometry measurements into quantitative distributions of mass isotopomers. Calibration curve approach has been, however, of very limited use in comprehensive analyses of biological systems, mainly suffering from the lack of extensive range of standard samples with accurately known isotopic enrichment. Here, we present a biological method capable of synthesizing specifically labeled amino acids. These amino acids have well-determined and estimable mass isotopomer distributions and thus can serve as standard samples. In this method, the bacterium strain Methylobacterium salsuginis sp. nov. was cultivated with partially 13C-labeled methanol as the only carbon source to produce 13C-enriched compounds. We show that the mass isotopomer distributions of the various biosynthesized amino acids are well determined and can be reasonably estimated based on proposed binomial approximation if the labeling state of the biomass reached an isotopic steady state. The interference of intramolecular inhomogeneity of 13C isotope abundances caused by biological isotope fractionation was eliminated by estimating average 13C isotope abundance. Further, the predictions are tested experimentally by mass spectrometry (MS) spectra of the labeled glycine, alanine, and aspartic acid. Most of the error in mass spectrometry measurements was less than 0.74 mol% in the test case, significantly reduced as compared with uncalibrated results, and this error is expected to be less than 0.4 mol% in real experiment as revealed by theoretical analysis.