Protein WP_076948622.1 in Rhodococcus qingshengii djl-6-2
Annotation: NCBI__GCF_002893965.1:WP_076948622.1
Length: 478 amino acids
Source: GCF_002893965.1 in NCBI
Candidate for 16 steps in catabolism of small carbon sources
Pathway | Step | Score | Similar to | Id. | Cov. | Bits | Other hit | Other id. | Other bits |
L-arabinose catabolism | xacF | hi | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 (characterized) | 42% | 98% | 373.6 | | | |
D-galacturonate catabolism | dopDH | hi | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 (characterized) | 42% | 98% | 373.6 | | | |
D-glucuronate catabolism | dopDH | hi | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 (characterized) | 42% | 98% | 373.6 | | | |
D-xylose catabolism | dopDH | hi | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 (characterized) | 42% | 98% | 373.6 | | | |
L-fucose catabolism | aldA | lo | lactaldehyde dehydrogenase (EC 1.2.1.22); 2,5-dioxovalerate dehydrogenase (EC 1.2.1.26) (characterized) | 38% | 98% | 313.9 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-rhamnose catabolism | aldA | lo | lactaldehyde dehydrogenase (EC 1.2.1.22); 2,5-dioxovalerate dehydrogenase (EC 1.2.1.26) (characterized) | 38% | 98% | 313.9 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-threonine catabolism | aldA | lo | lactaldehyde dehydrogenase (EC 1.2.1.22); 2,5-dioxovalerate dehydrogenase (EC 1.2.1.26) (characterized) | 38% | 98% | 313.9 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-valine catabolism | mmsA | lo | Methylmalonate-semialdehyde dehydrogenase [inositol] (EC 1.2.1.27) (characterized) | 33% | 97% | 257.7 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-lysine catabolism | amaB | lo | L-piperidine-6-carboxylate dehydrogenase; EC 1.2.1.21 (characterized, see rationale) | 34% | 92% | 248.1 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-isoleucine catabolism | iolA | lo | malonate-semialdehyde dehydrogenase (acetylating) (EC 1.2.1.18) (characterized) | 32% | 96% | 236.9 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
myo-inositol catabolism | mmsA | lo | malonate-semialdehyde dehydrogenase (acetylating) (EC 1.2.1.18) (characterized) | 32% | 96% | 236.9 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
propionate catabolism | iolA | lo | malonate-semialdehyde dehydrogenase (acetylating) (EC 1.2.1.18) (characterized) | 32% | 96% | 236.9 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-threonine catabolism | iolA | lo | malonate-semialdehyde dehydrogenase (acetylating) (EC 1.2.1.18) (characterized) | 32% | 96% | 236.9 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-valine catabolism | iolA | lo | malonate-semialdehyde dehydrogenase (acetylating) (EC 1.2.1.18) (characterized) | 32% | 96% | 236.9 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-arginine catabolism | astD | lo | Succinylglutamic semialdehyde dehydrogenase (EC 1.2.1.71) (characterized) | 32% | 92% | 218.8 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
L-citrulline catabolism | astD | lo | Succinylglutamic semialdehyde dehydrogenase (EC 1.2.1.71) (characterized) | 32% | 92% | 218.8 | Alpha-ketoglutaric semialdehyde dehydrogenase; alphaKGSA dehydrogenase; 2,5-dioxovalerate dehydrogenase; EC 1.2.1.26 | 42% | 373.6 |
Sequence Analysis Tools
View WP_076948622.1 at NCBI
Find papers: PaperBLAST
Find functional residues: SitesBLAST
Search for conserved domains
Find the best match in UniProt
Compare to protein structures
Predict transmenbrane helices: Phobius
Predict protein localization: PSORTb
Find homologs in fast.genomics
Fitness BLAST: loading...
Sequence
MSVAVRSLIGGAWVTGSDGTLTDTNPARPDEVVATGGRAGASEMATGVRVARIAANDWAR
LPMSARGAYLTRAAEIIESNAEAWGEELTREEGKTLAEGIGEVRRAAQILRYYGNAADRE
TGTVYSSPRSSEQILVTRKPLGVVGVVTPFNFPIAIPAWKIAPALVYGNTVVWKPAAAVP
LLAFRLAQALTEAGLPDGVLNMILADSAAGEALVDHPDLDAVTFTGSTGVGRKIAGSAAA
RGVPVQAEMGGKNAAVVLADADLELAVEQVMLGAFRSTGQKCTATSRLVLHESIAEEFLA
RLTVEVKQLRVGDPLDPDVTTGPVVSAYAQQSIIDGINRAVDGGARVVATAEAPQSGYYV
QPTVLELPSSTELWFEELFGPVLAVQRAATTADAFRLANEGEFGLSAALFTQDLTSALAG
IDDLDVGILHVNSESAGADPHVPFGGAKRSGYGPKEQGDAAREFFTHTTTVYLRGGRA
This GapMind analysis is from Sep 24 2021. The underlying query database was built on Sep 17 2021.
Links
Downloads
Related tools
About GapMind
Each pathway is defined by a set of rules based on individual steps or genes. Candidates for each step are identified by using
ublast (a fast alternative to protein BLAST)
against a database of manually-curated proteins (most of which are experimentally characterized) or by using
HMMer with enzyme models (usually from
TIGRFam). Ublast hits may be split across two different proteins.
A candidate for a step is "high confidence" if either:
- ublast finds a hit to a characterized protein at above 40% identity and 80% coverage, and bits >= other bits+10.
- (Hits to curated proteins without experimental data as to their function are never considered high confidence.)
- HMMer finds a hit with 80% coverage of the model, and either other identity < 40 or other coverage < 0.75.
where "other" refers to the best ublast hit to a sequence that is not annotated as performing this step (and is not "ignored").
Otherwise, a candidate is "medium confidence" if either:
- ublast finds a hit at above 40% identity and 70% coverage (ignoring otherBits).
- ublast finds a hit at above 30% identity and 80% coverage, and bits >= other bits.
- HMMer finds a hit (regardless of coverage or other bits).
Other blast hits with at least 50% coverage are "low confidence."
Steps with no high- or medium-confidence candidates may be considered "gaps."
For the typical bacterium that can make all 20 amino acids, there are 1-2 gaps in amino acid biosynthesis pathways.
For diverse bacteria and archaea that can utilize a carbon source, there is a complete
high-confidence catabolic pathway (including a transporter) just 38% of the time, and
there is a complete medium-confidence pathway 63% of the time.
Gaps may be due to:
- our ignorance of proteins' functions,
- omissions in the gene models,
- frame-shift errors in the genome sequence, or
- the organism lacks the pathway.
GapMind relies on the predicted proteins in the genome and does not search the six-frame translation. In most cases, you can search the six-frame translation by clicking on links to Curated BLAST for each step definition (in the per-step page).
For more information, see:
If you notice any errors or omissions in the step descriptions, or any questionable results, please let us know
by Morgan Price, Arkin group, Lawrence Berkeley National Laboratory