Giter VIP home page Giter VIP logo

rbc-gem's Introduction

RBC-GEM: Genome-scale metabolic model for the erythrocyte of Homo sapiens

GitHub version Zenodo Gitter chat memote tested

Description

TODO

Citation

If you use RBC-GEM in your research, please cite the following:

RBC-GEM: a Knowledge Base for Systems Biology of Human Red Blood Cell Metabolism; Zachary Brandon Haiman, Angelo D'Alessandro, Bernhard Palsson bioRxiv 2024.04.26.591249; doi: https://doi.org/10.1101/2024.04.26.591249

All the releases are also archived in [Zenodo] from which specific version can be cited if used.

Keywords

Utilisation: experimental data reconstruction; multi-omics integrative analysis; model template Field: metabolic-network reconstruction Type of model: reconstruction; curated Model source: iAB-RBC-283, Human-GEM Omic source: genomics; proteomics; metabolomics; lipidomics Taxonomic name: Homo sapiens Taxonomy ID: taxonomy:9606 Metabolic system: general metabolism Tissue: Blood Cell type: Red Blood Cell (erythrocyte) Condition: generic metabolism

Installation

With Python

The recommended method is to install rbc_gem_utils is to use pip to install the software. It is recommended to do this inside a virtual environment. If you have conda, follow the instructions for managing environments

  1. Clone the main branch of this repository, or download the latest release.

  2. Navigate to the code directory containing the pyproject.toml file and install the package:

    cd "/my/path/RBC-GEM/code"
    pip install "."
    
  3. Test your install:

     python -c "from rbc_gem_utils import show_versions; show_versions()"
    

Usage

TODO

Contributing

Contributions are always welcome! Please read the contributing guideline to get started.

Contributors

Code contributors are reported automatically by GitHub under Contributors, while other contributions come in as Issues.

rbc-gem's People

Contributors

z-haiman avatar

Watchers

 avatar

rbc-gem's Issues

Ammonia/Ammonium Reactions

Brief description:

Similar to Human-GEM issues #665 and PR#677, all reactions should depend on one form of NH3/NH4+ and have the deprotonation reaction. Human-GEM resolved issue by removing all of the NH4+ versions and keep all of the NH3 version.

Unlike Human-GEM though, sticking with NH4 over NH3 would be preferable.
Reasons for NH4 over NH3:

From Mohiuddin SS, Khattar D. Biochemistry, Ammonia. [Updated 2023 Feb 20]. In: StatPearls [Internet]. Treasure Island (FL): StatPearls Publishing; 2024 Jan-. Available from: https://www.ncbi.nlm.nih.gov/books/NBK541039/

In nature, ammonia exists as both NH3 and in the ionic form of ammonia as ammonium ion (NH4+). A buffering reaction: NH3+ H+ --> NH4+ is used to maintain the relative amount of each form. Under biological conditions, the pKa of this reaction is about 9.15 and this reaction occurs almost instantaneously. As a result, the majority of ammonia under physiological conditions exists as NH4+, and only about 1.7% of total ammonia presents as NH3 at pH 7.4.

Most physiological scenarios for RBCs will have NH4 prevalent over NH3 due to the pH.

Current status:

Identifier Reaction equation Genes
13DAMPPOX 13dampp_c + h2o_c + o2_c --> bamppald_c + h2o2_c + h_c + nh3_c GAPFILLING
2IMBUTDA 2ibut_c + h2o_c --> 2obut_c + h_c + nh3_c RIDA
3MOXTYROX 3moxtyr_c + h2o_c + o2_c --> 3mox4hpac_c + h2o2_c + h_c + nh3_c GAPFILLING
42A12BOOX dpam__L_c + h2o_c + o2_c --> 34dhpac_c + h2o2_c + h_c + nh3_c GAPFILLING
CSND csn_c + h2o_c --> nh3_c + ura_c GAPFILLING
CTPS1 atp_c + nh3_c + utp_c --> adp_c + ctp_c + h_c + pi_c CTPS1 or CTPS2
CYTD cytd_c + h2o_c --> nh3_c + uri_c

Proposed status/to-do list:

Identifier Reaction equation Genes
13DAMPPOX 13dampp_c + h2o_c + o2_c --> bamppald_c + h2o2_c + nh4_c GAPFILLING
2IMBUTDA 2ibut_c + h2o_c --> 2obut_c + nh4_c RIDA
3MOXTYROX 3moxtyr_c + h2o_c + o2_c --> 3mox4hpac_c + h2o2_c + nh4_c GAPFILLING
42A12BOOX dpam__L_c + h2o_c + o2_c --> 34dhpac_c + h2o2_c + nh4_c GAPFILLING
CSND csn_c + h2o_c + h_c --> nh4_c + ura_c GAPFILLING
CTPS1 atp_c + nh4_c + utp_c --> adp_c + ctp_c + 2 h_c + pi_c CTPS1 or CTPS2
CYTD cytd_c + h2o_c + h_c --> nh4_c + uri_c
NH4t nh4_c <=> nh4_e AQP3 or AQP1 or (RHAG and RHCE and RHD) or (RHAG and RHCE) or RHAG
  • Update stoichiometry for 13DAMPPOX.
  • Update stoichiometry for 2IMBUTDA.
  • Update stoichiometry for 3MOXTYROX.
  • Update stoichiometry for 42A12BOOX.
  • Update stoichiometry for CSND.
  • Update stoichiometry for CTPS1.
  • Update stoichiometry for CYTD.
  • Add NH4t as reversible transport by RHAG

Reaction reversibility updates

Brief description:

From Thuillier L. Purification and kinetic properties of human erythrocyte Mg2+-dependent inorganic pyrophosphatase. Biochim Biophys Acta. 1978;524(1):198-206. doi:10.1016/0005-2744(78)90118-3:

Biosynthetic reactions which produce pyrophosphate (polysaccharides, proteins, DNA and RNA synthesis} are subject to a "thermodynamic pull" [1--3] due to the subsequent complete and irreversible hydrolysis of pyrophosphate by widely distributed pyrophosphatases (standard free energy of hydrolysis: AG ° = --23.56 kJ • mo1-1 of PPi in the absence of divalent cations[4]).

The following enzyme catalyzed reactions would be considered to be reversible:

  • Orotate phosphoribosyltransferase, catalyzed by UMPS
  • Nicotinamide phosphoribosyltransferase, catalyzed by NAMPT
  • Hypoxanthine phosphoribosyltransferase, catalyzed by HPRT1

Current status:

Identifier Reaction equation Genes
ORPT orot5p_c + ppi_c <=> orot_c + prpp_c UMPS
NMNS h_c + ncam_c + prpp_c --> nmn_c + ppi_c NAMPT
HXPRT hxan_c + prpp_c --> imp_c + ppi_c HPRT1

Proposed status/to-do list:

Identifier Reaction equation Genes
ORPT orot_c + prpp_c <=> orot5p_c + ppi_c UMPS
NMNS h_c + ncam_c + prpp_c <=> nmn_c + ppi_c NAMPT
HXPRT hxan_c + prpp_c <=> imp_c + ppi_c HPRT1
  • Reverse stoichiometry for ORPT.
  • Set NMNS as reversible.
  • Set HXPRT as reversible.

Hydrogen sulfide reactions

Brief description:

All reactions should depend on one form of H2S/HS(-) and have the deprotonation reaction. Issue can be resolved by removing all of the H2S versions and keep all of the HS(-) versions.

Reasons for HS- over H2S:

From Wang R. Physiological implications of hydrogen sulfide: a whiff exploration that blossomed. Physiol Rev. 2012 Apr;92(2):791-896. doi: 10.1152/physrev.00017.2011. PMID: 22535897.

The definition of what constitutes free sulfide, as opposed to bound sulfide, should also be noted. Free sulfide is dissolved H2S gas, which is a weak acid and in solution exists in the equilibrium H2S ↔ HS− ↔ S2−. With the pKa1 ∼7.0 and the pKa2 >17, there is essentially no S2− in biological tissues, nearly equal amounts of H2S and HS− within the cell, and approximately a 20% H2S/80% HS− ratio in extracellular fluid and plasma at 37°C and pH 7.4. Due to the temperature sensitivity of pKa1, the H2S/HS− ratio remains nearly constant in blood and tissues of ectothermic vertebrates over a wide range of body temperatures. For example, in trout blood at 10°C, the ratio is 15% H2S/85% HS−.

From Jennings ML. Transport of H2S and HS(-) across the human red blood cell membrane: rapid H2S diffusion and AE1-mediated Cl(-)/HS(-) exchange. Am J Physiol Cell Physiol. 2013 Nov 1;305(9):C941-50. doi: 10.1152/ajpcell.00178.2013. Epub 2013 Jul 17. PMID: 23864610; PMCID: PMC4042536.

The data indicate that HS(-) is a very good substrate for AE1; the Cl(-)/HS(-) exchange rate is about one-third as rapid as Cl(-)/HCO3(-) exchange. The H2S permeability coefficient must also be high (>10(-2) cm/s, half time <0.003 s) to account for the pH equilibration data.

Most physiological scenarios for RBCs will have HS- prevalent over H2S due to the pH, and Band 3 can exchange HS(-).

Current status:

Identifier Reaction equation Genes
H2SMT amet__L_c + h2s_c --> ahcys__L_c + ch4s_c + h_c TMT1A
HTSULDHLA dhlipoate_c + htsul_c --> h2s_c + h_c + lipoate_c + so3_c MPST
MTOX ch4s_c + h2o_c + o2_c --> fald_c + h2o2_c + h2s_c SELENBP1
SPODM2 h2s_c + 2.0 o2_c --> 2.0 h_c + so4_c SOD1

Proposed status/to-do list:

Identifier Reaction equation Genes
H2SMT amet__L_c + sh1_c --> ahcys__L_c + ch4s_c TMT1A
HTSULDHLA dhlipoate_c + htsul_c --> sh1_c + 2 h_c + lipoate_c + so3_c MPST
MTOX ch4s_c + h2o_c + o2_c --> fald_c + h2o2_c + sh1_c + h_c SELENBP1
SPODM2 sh1_c + 2.0 o2_c --> h_c + so4_c SOD1
  • Update stoichiometry for H2SMT.
  • Update stoichiometry for HTSULDHLA.
  • Update stoichiometry for MTOX.
  • Update stoichiometry for SPODM2.

Studies exploring human G6PD-deficient individuals/RBCs in acute exercise conditions

Brief explanation of proposed change

Inclusion of additional references for G6PD. Studies explore G6PD-deficient RBCs

PMID or DOI Explanation Title
PMID:16888458 Important (possibly first) human study highlighting "G6PD-deficient individuals are able to perform short- and long-duration acute exercise without experiencing greater oxidative stress than non-G6PD-deficient individuals. Findings demonstrate "that despite the theoretically lower capacity of G6PD-deficient individuals to resist perturbations in their redox status, they are not more susceptible to oxidative stress" Exercise-induced oxidative stress in G6PD-deficient individuals
PMID:19997026 Follow-up study with additional results complimenting PMID:16888458 Comparison between glucose-6-phosphate dehydrogenase-deficient and normal individuals after eccentric exercise

Affected items

Entry Type PubMed ID or DOI
G6PD Gene pubmed:16888458;pubmed:19997026

Anything else:

Additional evidence for CD36

Brief explanation of proposed change

Inclusion of additional references for CD36.

PMID or DOI Explanation Title
PMID: 38326223 "[...] detected CD36 expression on RBCs and reticulocytes from 20 blood donors. By mining reticulocyte and RBC datasets, we found evidence for CD36-derived peptides enriched in the membrane fractions." Evidence that CD36 is expressed on red blood cells and constitutes a novel blood group system of clinical importance

Affected items

Add the following publications as evidence for existing genes:

Entry Type PubMed ID or DOI
CD36 Gene pubmed:38326223

Anything else:

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.