PE16 Family assigned · medium auto-curated
H37Rv Rv1430 · MTBC0 - ·
528 aa · 1606386–1607972 (+) ·
RefSeq YP_177810.1
Annotation: from legacy to revised
| Legacy (H37Rv / Mycobrowser) | PE family protein PE16 |
|---|---|
| MTBC0 PGAP re-annotation | — |
| Revised (this work) | PE family protein PE16. Pfam: PE (PF00934.26), PE-PPE (PF08237.18). |
Auto-curated: this verdict and function were generated by rules from PGAP + Pfam + Foldseek and have not been hand-reviewed.
Annotated on the H37Rv protein: this gene has no 1:1 ancestral MTBC0 anchor (PE/PPE, paralogue, IS element, or otherwise unanchored CDS).
Curated reference (UniProt)
| UniProt |
L7N697
SwissProt · reviewed
· Evidence at protein level
|
|---|---|
| UniProt name | Esterase PE16 |
| EC (curated) |
EC 3.1.1.-
|
| Curated function | Esterase that hydrolyzes short to medium chain fatty acid esters with the highest specific activity for p-nitrophenyl caproate (pNPC6). Has lower activity with p-nitrophenyl caprylate (pNPC8) and p-nitrophenyl butyrate (pNPC4). Has weak activity with p-nitrophenyl caprate (pNPC10) and p-nitrophenyl laurate (pNPC12). Does not possess lipolytic activity and cutinase activity. |
Functional vocabulary (eggNOG-mapper, orthology transfer)
| COG category |
N Cell motility
|
|---|---|
| Preferred name | PE16 |
| eggNOG description | PE-PPE domain |
| Orthologous group | COG3391 |
Orthology-based transfer (eggNOG 5.0.2, diamond). EC/KO/GO/CAZy are computed annotations, not manual curation; cross-check against the primary literature before treating a specific reaction as established.
Conservation & selection (intra-MTBC, 145 209 strains)
| pN/pS | 0.718 · relaxed/neutral |
|---|---|
| Polymorphic sites (≥ 0.1% of strains) | 2 synonymous, 4 missense, 0 nonsense, 2 frameshift |
| Disruption | 2 distinct premature-stop/frameshift site(s); most common in 0.11% of strains (154) · clonal |
pN/pS from segregating SNPs (singletons removed) normalised by possible sites. Low pN/pS = purifying selection (a strong signal that a "hypothetical" is a real, constrained gene). A high pN/pS is ambiguous: relaxed constraint or positive selection (drug resistance, antigenic variation) inflate it; e.g. rpoB/katG/pncA score high here for resistance, not loss of function. A clonal disruption (one allele over a clade) suggests lineage pseudogenisation; a convergent one (many independent alleles) is typical of resistance loss-of-function.
Domains (Pfam, hmmscan --cut_ga)
| Pfam | Accession | i-Evalue | Residues | Description |
|---|---|---|---|---|
PE | PF00934.26 | 3.6e-33 | 4–94 | PE family |
PE-PPE | PF08237.18 | 1.0e-88 | 146–369 | PE-PPE domain |
Functional interaction network (STRING v12, guilt-by-association)
Closest characterised functional partner: Rv1432 (dehydrogenase), medium confidence from genomic context alone (score 551 excluding text-mining).
| Partner | Product | Score | No text-mining | Channels (≥400) |
|---|---|---|---|---|
Rv1431 hyp |
hypothetical protein | 551 | 551 ctx | neighborhood:551 |
Rv1432 |
dehydrogenase | 551 | 551 ctx | neighborhood:551 |
Rv1364c |
sigma factor regulatory protein | 520 | 515 | coexpression:500 |
Rv1429 hyp |
hypothetical protein | 457 | 457 ctx | neighborhood:457 |
Rv3005c hyp |
hypothetical protein | 473 | 100 | textmining:439 |
Rv0024 |
NLP/P60 family protein | 534 | 63 | textmining:523 |
Rv1337 |
integral membrane protein | 425 | 63 | textmining:412 |
Rv2431c PE25 |
PE family protein PE25 | 775 | 41 | textmining:775 |
Rv3893c PE36 |
PE family protein PE36 | 626 | 41 | textmining:626 |
Rv2430c PPE41 |
PPE family protein PPE41 | 411 | 41 | textmining:411 |
STRING combines evidence channels (neighborhood, fusion, cooccurrence, coexpression, experimental, database, text-mining) into a 0–1000 score. The ctx badge marks edges carried by the genomic-context channels (conserved neighborhood, fusion, phylogenetic co-occurrence), which are independent of orthology and structure and the strongest signal for an unknown gene. The no text-mining column recomputes the score from data alone, so a link that does not depend on the literature is visible. Association is a function hypothesis, not proof: corroborate with the operon context and the primary literature before assigning a function.
Evidence
- Annotation from H37Rv (no MTBC0 1:1 anchor; H37Rv protein used): PE family protein PE16
- Pfam (hmmscan --cut_ga): PE PF00934.26 (E=4e-33), PE-PPE PF08237.18 (E=1e-88)
- (auto-curated by rules from PGAP + Pfam + Foldseek; not hand-reviewed)
Sources
- Ancestral sequence & coordinates: Harrison LB et al. (2024), An imputed ancestral reference genome for the MTBC, doi:10.1101/2023.09.07.556366
- Product annotation: NCBI PGAP on MTBC0; legacy from H37Rv NC_000962.3 (RefSeq YP_177810.1)
- Domains: Pfam-A via hmmscan --cut_ga — PE (PF00934.26), PE-PPE (PF08237.18)
- Sequence-level signal: ESM Atlas (EvolutionaryScale × BioHub) — exploratory
- Controlled vocabulary: eggNOG-mapper 2.1.12 (Cantalapiedra et al. 2021,
doi:10.1093/molbev/msab293), eggNOG 5.0 DB
(Huerta-Cepas et al. 2019) — OG
COG3391 - Curated reference: UniProt L7N697 (SwissProt, reviewed; Evidence at protein level)
- Intra-MTBC selection: pN/pS and disruption from SPDI variants of 145 209 MTBC strains (this work, local collection vs H37Rv NC_000962.3)
- Interaction network: STRING v12.0 (Szklarczyk et al. 2023,
doi:10.1093/nar/gkac1000), taxon 83332, CC-BY 4.0 —
10 functional partner(s); context anchor
Rv1432 - Primary literature: none located yet; annotation rests on the domain/homology sources above.
Ancestral MTBC0 protein sequence
>H37Rv|Rv1430|PE16 MSFVFAVPEMVAATASDLASLGAALSEATAAAAIPTTQVLAAAADEVSAAIAELFGAHGQEFQALSAQASAFHDRFVRALSAAAGWYVDAEAANAALVDTAATGASELGSGGRTALILGSTGTPRPPFDYMQQVYDRYIAPHYLGYAFSGLYTPAQFQPWTGIPSLTYDQSVAEGAGYLHTAIMQQVAAGNDVVVLGFSQGASVATLEMRHLASLPAGVAPSPDQLSFVLLGNPNNPNGGILARFPGLYLQSLGLTFNGATPDTDYATTIYTTQYDGFADFPKYPLNILADVNALLGIYYSHSLYYGLTPEQVASGIVLPVSSPDTNTTYILLPNEDLPLLQPLRGIVPEPLLDLIEPDLRAIIELGYDRTGYADVPTPAALFPVHIDPIAVPPQIGAAIGGPLTALDGLLDTVINDQLNPVVTSGIYQAGAELSVAAAGYGAPAGVTNAIFIGQQVLPILVEGPGALVTADTHYLVDAIQDLAAGDLSGFNQNLQLIPATNIALLVFAAGIPAVAAVAILTGQDFPV