Skip to content

Major mass-rebalance of HumanGEM - #131

Merged
haowang-bioinfo merged 59 commits into
develfrom
rebalanceModel
Nov 20, 2019
Merged

Major mass-rebalance of HumanGEM#131
haowang-bioinfo merged 59 commits into
develfrom
rebalanceModel

Conversation

@JonathanRob

@JonathanRob JonathanRob commented Sep 19, 2019

Copy link
Copy Markdown
Collaborator

Main improvements in this PR:

This PR implements a substantial re-balancing of reactions in HumanGEM, resulting in a model that is 93.4% charge-balanced, 99.4% mass-balanced, and 100% stoichiometrically consistent. The changes also address issues #43, #58, and #60.

The manual curation process identified and implemented the following changes:

Curation description Number
Updated metabolite formulas and/or charges 465
Updated reaction equations 724
Identification and removal of duplicated reactions 203
Inactivation (LB=UB=0) of reactions due to mass inconsistency 196
Repair and re-activation of previously inactivated reactions 10
Identification and removal of duplicated metabolites 34
Addition of new metabolites 6
Addition of existing metabolites to new compartments 20
Addition of new reactions 21
Updated metabolite annotations (external IDs) 26
Updated reaction annotations (external IDs) 15
Updated reaction grRules 28


The changes are implemented in the rebalanceHumanGEM script, which reads change information organized into 12 different tsv files located in the /ComplementaryData/modelCuration/fullRebalance directory.

Other changes in this PR:

  • A new field rxnHMRID is added to the humanGEMRxnAssoc_lip.JSON reaction annotation file, which contains HMR1 and HMR2 reaction IDs associated with each reaction (if there are any). This provides more explicit association to these models, and was also necessary to document the removal of some duplicated HMR reactions.

I hereby confirm that I have:

  • Tested my code on my own computer for running the model
  • Selected devel as a target branch

@JonathanRob

Copy link
Copy Markdown
Collaborator Author

This PR should not be merged until after #130 has been merged.

@JonathanRob
JonathanRob changed the base branch from master to devel September 19, 2019 15:40
@JonathanRob JonathanRob reopened this Oct 3, 2019
@JonathanRob
JonathanRob requested a review from avlant October 3, 2019 06:48

@avlant avlant left a comment

Copy link
Copy Markdown

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Very nice

haowang-bioinfo and others added 6 commits November 19, 2019 14:10
- These files had been dealt with in #141, now only keep the xlsx version.
- Newly assigned KEGG rxn id to HMR_4284 through manual curation
- adjust and use this function for updating grRules and other related fields
- add notes to the model change log files
Comment on lines +1 to +7
metNames newFormula newCharge
alpha-D-mannosyl-beta-D-mannosyl-diacetylchitobiosyldiphosphodolichol C128H210N2O27P2 -2
beta-D-mannosyldiacetylchitobiosyldiphosphodolichol C122H200N2O22P2 -2
dehydrodolichol C100H162O 0
dehydrodolichol-diphosphate C100H161O7P2 -3
dehydrodolichol-phosphate C100H161O4P -2
dolichol C100H164O 0

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This file contains significant updates to metabolite formulas and charges. It would very nice to add two more columns: oldFormula and oldCharge (which actually can be found in the model change log files).

Comment on lines +2 to +12
HMR_5164 ENSG00000072274
HMR_5154 ENSG00000084674
HMR_5159 ENSG00000110245
HMR_5162 ENSG00000122194
HMR_5165 ENSG00000130203
HMR_5157 ENSG00000130208
HMR_5160 ENSG00000171560
HMR_5163 ENSG00000180210
HMR_5152 ENSG00000196136
HMR_5153 ENSG00000197249
HMR_5158 ENSG00000234906

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If there are references that support these grRule changes, the pmid and suggested confidenceScores can be added as additional columns in this file.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants