Protein WP_043917906.1 in Jannaschia aquimarina GSW-M26
Annotation: NCBI__GCF_000877395.1:WP_043917906.1
Length: 239 amino acids
Source: GCF_000877395.1 in NCBI
Candidate for 4 steps in catabolism of small carbon sources
Pathway | Step | Score | Similar to | Id. | Cov. | Bits | Other hit | Other id. | Other bits |
4-hydroxybenzoate catabolism | mhpE | hi | 4-hydroxy-2-oxo-heptane-1,7-dioate aldolase; 2,4-dihydroxyhept-2-ene-1,7-dioic acid aldolase; HHED aldolase; 4-hydroxy-2-ketoheptane-1,7-dioate aldolase; HKHD aldolase; EC 4.1.2.52 (characterized) | 47% | 93% | 190.3 | 2-keto-3-deoxy-L-rhamnonate aldolase (EC 4.1.2.53) | 44% | 165.2 |
L-tryptophan catabolism | mhpE | hi | 4-hydroxy-2-oxo-heptane-1,7-dioate aldolase; 2,4-dihydroxyhept-2-ene-1,7-dioic acid aldolase; HHED aldolase; 4-hydroxy-2-ketoheptane-1,7-dioate aldolase; HKHD aldolase; EC 4.1.2.52 (characterized) | 47% | 93% | 190.3 | 2-keto-3-deoxy-L-rhamnonate aldolase (EC 4.1.2.53) | 44% | 165.2 |
L-rhamnose catabolism | LRA4 | med | 2-keto-3-deoxy-L-rhamnonate aldolase; KDR aldolase; EC 4.1.2.53; 2-dehydro-3-deoxyrhamnonate aldolase (uncharacterized) | 46% | 87% | 172.2 | 4-hydroxy-2-oxo-heptane-1,7-dioate aldolase; 2,4-dihydroxyhept-2-ene-1,7-dioic acid aldolase; HHED aldolase; 4-hydroxy-2-ketoheptane-1,7-dioate aldolase; HKHD aldolase; EC 4.1.2.52 | 47% | 190.3 |
D-glucuronate catabolism | garL | med | 5-keto-4-deoxy-D-glucarate aldolase; KDGluc aldolase; KDGlucA; 2-dehydro-3-deoxy-D-glucarate aldolase; 2-keto-3-deoxy-D-glucarate aldolase; 5-dehydro-4-deoxy-D-glucarate aldolase; Alpha-keto-beta-deoxy-D-glucarate aldolase; EC 4.1.2.20 (characterized) | 40% | 95% | 164.5 | 4-hydroxy-2-oxo-heptane-1,7-dioate aldolase; 2,4-dihydroxyhept-2-ene-1,7-dioic acid aldolase; HHED aldolase; 4-hydroxy-2-ketoheptane-1,7-dioate aldolase; HKHD aldolase; EC 4.1.2.52 | 47% | 190.3 |
Sequence Analysis Tools
View WP_043917906.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
MSLKSRVLAREPLTGLWQVLPGPVAAEIAARAGFDFLVLDGEHGPWDPSDLRARLIAVPD
AIVRVPANDPVWIKQALDLGAMTVLVPMVHDAEGARAAVAAARYPPDGIRGHGAFVSRAS
AYGQNAGYVTRANEAVGVWVQAESRPALSDLEAICGVEGVDCVFIGPADLAADMGLTPDA
PEVLTAIEDAIGRIVATGTACGAFGDPALYPRWRELGATILSAGVDGSVLAAALARLAP
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