GapMind for Amino acid biosynthesis

 

Alignments for a candidate for metH in Geobacter lovleyi SZ

Align Methionine synthase; EC 2.1.1.13; 5-methyltetrahydrofolate--homocysteine methyltransferase; Methionine synthase, vitamin-B12 dependent; MS (uncharacterized)
to candidate WP_012470188.1 GLOV_RS10575 methionine synthase

Query= curated2:O33259
         (1192 letters)



>NCBI__GCF_000020385.1:WP_012470188.1
          Length = 1169

 Score =  771 bits (1992), Expect = 0.0
 Identities = 459/1193 (38%), Positives = 691/1193 (57%), Gaps = 60/1193 (5%)

Query: 13   LLDVLSQRVMVGDGAMGTQLQAADLTLDDFRG--LEGCNEILNETRPDVLETIHRNYFEA 70
            +L  LSQ +++ DGAMGTQLQA +LT  DF G   EGCNE+L  TRPDV+E +H+ Y EA
Sbjct: 6    VLQALSQCILILDGAMGTQLQARNLTAADFGGEAYEGCNEMLVLTRPDVIEDVHKAYLEA 65

Query: 71   GADAVETNTFGCNLSNLGDYDIADRIRDLSQKGTAIARRVADELGSPDRKRYVLGSMGPG 130
            GAD VET +FG     L +Y +AD++ +L++    +AR+  D   +PD+ R+V GSMGP 
Sbjct: 66   GADIVETCSFGSTDIVLAEYGLADKVFELNKAAALVARKACDAYSTPDKPRFVAGSMGPT 125

Query: 131  TK-LPTLGHTEYAVIRDAYTEAALGMLDGGADAILVETCQDLLQLKAAVLGSRRAMTRAG 189
            T+ +   G   +  +  A+ +  +G+L GG D +L+ET QD L LKA   G R A  + G
Sbjct: 126  TRTISVTGGVTFEQLVTAFRDQTIGLLAGGVDLLLLETAQDTLNLKAGAEGIRLAFEQTG 185

Query: 190  RHIPVFAHVTVETTGTMLLGSEIGAALTAV----EPLGVDMIGLNCATGPAEMSEHLRHL 245
            + +P+    T+E TGTML G  + A   ++    E LG+  IGLNCATGP  M++HLR L
Sbjct: 186  QRVPLMVSGTIEPTGTMLAGQNVEALYASLTHLEENLGLISIGLNCATGPEFMTDHLRAL 245

Query: 246  SRHARIPVSVMPNAGLPVLGAKGAEYPLLPDELAEALAGFIAEFGLSLVGGCCGTTPAHI 305
            S  A   +SV PNAGLP    +   Y   PD LA+ L+ F+ E  L+++GGCCGT+PAHI
Sbjct: 246  SELATCHISVYPNAGLP---DENGNYAESPDSLAQKLSRFVDEGWLNIIGGCCGTSPAHI 302

Query: 306  REVAAAVANIKRPERQVSYEPSVSSLYTAIPFAQDASVLVIGERTNANGSKGFREAMIAE 365
              +A   A  ++P + V+    + S    +   +D   +++GERTN  GS+ F+  ++A 
Sbjct: 303  AAIANMAAG-RQPRKPVTIRRRLVSGIEPLFIEEDNRPILVGERTNVIGSRKFKNMIVAG 361

Query: 366  DYQKCLDIAKDQTRDGAHLLDLCVDYVGRDGVADMKALASRLATSSTLPIMLDSTETAVL 425
             +++  +IA+ Q + GA ++D+CV    RD +ADM+ +   L      PIM+DST+  V 
Sbjct: 362  QFEEAAEIARAQVKTGAQVIDICVANPDRDELADMEQMLVYLPKKVRAPIMIDSTDAKVT 421

Query: 426  QAGLEHLGGRCAINSVNYEDGDGPESRFAKTMALVAEHGAAVVALTIDEE---GQARTAQ 482
            +  L+ L G+C INS+N EDG   E RFAK   L+  +G AVV   IDE+   G A T Q
Sbjct: 422  ELALQRLQGKCVINSINLEDG---EERFAKVAPLLRHYGGAVVVGCIDEDPANGMAVTRQ 478

Query: 483  KKVEIAERLINDITGNWGVDESSILIDTLTFTIATGQEESRRDGIETIEAIRELKKRHPD 542
            +K+E+A+R  + +   +G+    ++ D L F + +G +      +ETIE +R + +  P 
Sbjct: 479  RKLEVAQRSYDLLVNKYGLRPEDLIFDPLVFPVGSGDDNYIGSAVETIEGVRLITETFPQ 538

Query: 543  VQTTLGLSNISFGLNPAARQVLNSVFLHECQEAGLDSAIVHASKILPMNRIPEEQRNVAL 602
              T LG+SN+SFGL  A R+VLN+VFL+ C +AGL  AIV+  K+     IPE +R +A 
Sbjct: 539  CSTILGISNVSFGLPAAGREVLNAVFLYHCVKAGLGYAIVNTEKLERYASIPEAERRLAE 598

Query: 603  DLVYDRRREDYDPLQELMRLFEGVSAASSKEDRLAELAGLPLFERLAQRIVDGERNGLDA 662
            DL++ R     DP+      F      S      A  A LPL ERL   I++G ++GL A
Sbjct: 599  DLIFWR---GADPVAAFAAAFRDKKPVSH-----APAAELPLDERLPLYIIEGSKDGLTA 650

Query: 663  DLDEAMTQ-KPPLQIINEHLLAGMKTVGELFGSGQMQLPFVLQSAEVMKAAVAYLEPHME 721
            DLD A+T+   PL+IIN  L+AGM  VG LF   Q+ +  VLQSAE MKAAV++LEPH+ 
Sbjct: 651  DLDAALTRGDKPLEIINGPLMAGMAEVGRLFNDNQLIVAEVLQSAEAMKAAVSHLEPHLS 710

Query: 722  RSDDDSGKGRIVLATVKGDVHDIGKNLVDIILSNNGYEVVNIGIKQPIATILEVAEDKSA 781
            + D+ S KG+++LATVKGDVHDIGKNLV+IILSNNG+EV+N+GIK     ++  A+ ++ 
Sbjct: 711  K-DETSSKGKMLLATVKGDVHDIGKNLVEIILSNNGFEVINLGIKVGPEELIAAAKKENP 769

Query: 782  DVVGMSGLLVKSTVVMKENLEEMNTRGVAEKFPVLLGGAALTRSYVENDLAEIYQGEVHY 841
            D +G+SGLLVKS + M     ++   G+    P+L+GGAAL+R++ +  +   Y G V Y
Sbjct: 770  DFIGLSGLLVKSALQMIVTAADLKAAGI--DAPLLVGGAALSRAFADTRITPEYNGPVLY 827

Query: 842  ARDAFEGLKLMDTIMS-AKRGEAPDENSPEAIKAREKEAERKARHQRSKRIAAQRKAAEE 900
            A+DA  GL+L + ++  A R +   +       AR++EA        S +IAA + A + 
Sbjct: 828  AKDAMAGLELANQLVDPALRQQLMLD------LARQQEA--------SAKIAAAKAAGQS 873

Query: 901  PVEV-PERSDVAADIEVPAPPFWGSRIVKGLAVADYTGLLDERALFLGQWGLRGQRGGEG 959
             V      S ++++  +   P     I++ + V      L+ + L+    GL G    + 
Sbjct: 874  TVATGSTHSAISSNAPILVAPDLEQHILRDIPVGQIIPYLNRQMLYTKHLGLTGS-VDKL 932

Query: 960  PSYEDLVETEGRPRLRYWLDRLSTDGILAHAAVVYGYFPAVSEGNDIVVLTEPKPDAPVR 1019
             + +D   T+    +   L+R+  +G++   A +Y ++ A  +GND+++  +   +A   
Sbjct: 933  LAGQDEKATKLHLTVEAMLERVLQEGLIKPQA-IYRFYQANGDGNDLILFNQDGSEAT-- 989

Query: 1020 YRFHFPRQQRGRFLCIADFIRSRELAAERGEVDVLPFQLVTMGQPIADFANELFASNAYR 1079
             RF  PRQ+ G  LC+ADF+R         E D +    VT GQ + + + +  A   Y 
Sbjct: 990  -RFTLPRQKSGEQLCVADFVR----PLSGTEKDTMALFAVTCGQGVRELSEQWKAEGDYL 1044

Query: 1080 DYLEVHGIGVQLTEALAEYWHRRIREELKFSGDRAMAAEDPEAKEDYFKLGYRGARFAFG 1139
            +   +  + +++ EA AEY H+RIR       D+ +        +  F   Y+G R +FG
Sbjct: 1045 NSHLLQALALEMAEATAEYLHKRIRTSWGIVDDQTL------TMKQLFNAEYQGIRVSFG 1098

Query: 1140 YGACPDLEDRAKMMALLEPERIGVTLSEELQLHPEQSTDAFVLHHPEAKYFNV 1192
            Y ACP+L+D+ K+ +LL+PE+IG+ L+E   + PE S  A V HHPE KYF++
Sbjct: 1099 YPACPNLDDQKKLFSLLKPEQIGINLTEGDMMDPEASVSALVFHHPEGKYFDL 1151


Lambda     K      H
   0.317    0.135    0.386 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Sequences: 1
Number of Hits to DB: 2978
Number of extensions: 121
Number of successful extensions: 17
Number of sequences better than 1.0e-02: 1
Number of HSP's gapped: 1
Number of HSP's successfully gapped: 1
Length of query: 1192
Length of database: 1169
Length adjustment: 47
Effective length of query: 1145
Effective length of database: 1122
Effective search space:  1284690
Effective search space used:  1284690
Neighboring words threshold: 11
Window for multiple hits: 40
X1: 16 ( 7.3 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.7 bits)
S2: 58 (26.9 bits)

Align candidate WP_012470188.1 GLOV_RS10575 (methionine synthase)
to HMM TIGR02082 (metH: methionine synthase (EC 2.1.1.13))

# hmmsearch :: search profile(s) against a sequence database
# HMMER 3.3.1 (Jul 2020); http://hmmer.org/
# Copyright (C) 2020 Howard Hughes Medical Institute.
# Freely distributed under the BSD open source license.
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query HMM file:                  ../tmp/path.aa/TIGR02082.hmm
# target sequence database:        /tmp/gapView.5930.genome.faa
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query:       TIGR02082  [M=1182]
Accession:   TIGR02082
Description: metH: methionine synthase
Scores for complete sequences (score includes all domains):
   --- full sequence ---   --- best 1 domain ---    -#dom-
    E-value  score  bias    E-value  score  bias    exp  N  Sequence                                 Description
    ------- ------ -----    ------- ------ -----   ---- --  --------                                 -----------
          0 1227.0   0.0          0 1226.8   0.0    1.0  1  lcl|NCBI__GCF_000020385.1:WP_012470188.1  GLOV_RS10575 methionine synthase


Domain annotation for each sequence (and alignments):
>> lcl|NCBI__GCF_000020385.1:WP_012470188.1  GLOV_RS10575 methionine synthase
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ! 1226.8   0.0         0         0       2    1181 ..      11    1150 ..      10    1151 .. 0.95

  Alignments for each domain:
  == domain 1  score: 1226.8 bits;  conditional E-value: 0
                                 TIGR02082    2 nkrilvlDGamGtqlqsanLteadFrgeeadlarelkGnndlLnltkPeviaaihrayfeaGaDive 68  
                                                ++ il+lDGamGtqlq+ nLt+adF ge       ++G+n++L+lt+P+vi+++h+ay+eaGaDive
  lcl|NCBI__GCF_000020385.1:WP_012470188.1   11 SQCILILDGAMGTQLQARNLTAADFGGE------AYEGCNEMLVLTRPDVIEDVHKAYLEAGADIVE 71  
                                                788************************7......9******************************** PP

                                 TIGR02082   69 tntFnsteialadYdledkayelnkkaaklarevadeftltpekkRfvaGslGPtnklatlspdver 135 
                                                t  F+st i+la+Y+l+dk+ elnk+aa +ar+++d ++ tp+k+RfvaGs+GPt+++ +++     
  lcl|NCBI__GCF_000020385.1:WP_012470188.1   72 TCSFGSTDIVLAEYGLADKVFELNKAAALVARKACDAYS-TPDKPRFVAGSMGPTTRTISVT----- 132 
                                                ***************************************.*********************9..... PP

                                 TIGR02082  136 pefrnvtydelvdaYkeqvkglldGGvDllLietvfDtlnakaalfaveevfeekgrelPilisgvi 202 
                                                     vt+++lv a+++q+ gll GGvDllL+et +Dtln+ka   +++ +fe+ g+++P+++sg+i
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  133 ---GGVTFEQLVTAFRDQTIGLLAGGVDLLLLETAQDTLNLKAGAEGIRLAFEQTGQRVPLMVSGTI 196 
                                                ...789**********************************************************876 PP

                                 TIGR02082  203 vdksGrtLsGqtleaflaslehae....ililGLnCalGadelrefvkelsetaealvsviPnaGLP 265 
                                                 + +G++L Gq++ea +asl h e    ++++GLnCa+G++ + +++++lse a++ +sv+PnaGLP
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  197 -EPTGTMLAGQNVEALYASLTHLEenlgLISIGLNCATGPEFMTDHLRALSELATCHISVYPNAGLP 262 
                                                .899***************987542222899************************************ PP

                                 TIGR02082  266 nalgeYdltpeelakalkefaeegllnivGGCCGttPehiraiaeavkdikprkrqeleeksvlsgl 332 
                                                ++ g+Y + p++la++l  f+ eg lni+GGCCGt P+hi+aia+++++ +prk  ++    ++sg+
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  263 DENGNYAESPDSLAQKLSRFVDEGWLNIIGGCCGTSPAHIAAIANMAAGRQPRKPVTIRR-RLVSGI 328 
                                                *****************************************************9887765.5689** PP

                                 TIGR02082  333 ealkiaqessfvniGeRtnvaGskkfrklikaedyeealkiakqqveeGaqilDinvDevllDgead 399 
                                                e+l i++++  + +GeRtnv Gs+kf+++i a+++eea +ia+ qv+ Gaq++Di+v +  +D  ad
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  329 EPLFIEEDNRPILVGERTNVIGSRKFKNMIVAGQFEEAAEIARAQVKTGAQVIDICVANPDRDELAD 395 
                                                ******************************************************************* PP

                                 TIGR02082  400 mkkllsllasepdiakvPlmlDssefevleaGLkviqGkaivnsislkdGeerFlekaklikeyGaa 466 
                                                m+++l  l ++  + + P+m+Ds++ +v e  L+++qGk+++nsi+l+dGeerF++ a l+++yG a
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  396 MEQMLVYLPKK--V-RAPIMIDSTDAKVTELALQRLQGKCVINSINLEDGEERFAKVAPLLRHYGGA 459 
                                                ****9999888..5.9*************************************************** PP

                                 TIGR02082  467 vvvmafDee...GqartadkkieiakRayklltekvgfppediifDpniltiatGieehdryaidfi 530 
                                                vvv  +De+   G+a t+++k+e+a+R y+ll++k+g+ ped+ifDp++++++ G ++    a+++i
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  460 VVVGCIDEDpanGMAVTRQRKLEVAQRSYDLLVNKYGLRPEDLIFDPLVFPVGSGDDNYIGSAVETI 526 
                                                ********6666******************************************************* PP

                                 TIGR02082  531 eaireikeelPdakisgGvsnvsFslrgndavRealhsvFLyeaikaGlDmgivnagklavyddidk 597 
                                                e +r i e++P++ + +G+snvsF+l    a Re+l++vFLy+++kaGl  +ivn+ kl+ y+ i++
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  527 EGVRLITETFPQCSTILGISNVSFGLP--AAGREVLNAVFLYHCVKAGLGYAIVNTEKLERYASIPE 591 
                                                ***************************..999*********************************** PP

                                 TIGR02082  598 elrevvedlildrrreatekLlelaelykgtkeksskeaqeaewrnlpveeRLeralvkGeregiee 664 
                                                  r ++edli+ r     +++  +a  ++++k        +a  ++lp++eRL   +++G ++g+ +
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  592 AERRLAEDLIFWR---GADPVAAFAAAFRDKKPV-----SHAPAAELPLDERLPLYIIEGSKDGLTA 650 
                                                ***********98...57899********99887.....457899********************** PP

                                 TIGR02082  665 dleearkklkapleiiegpLldGmkvvGdLFGsGkmfLPqvvksarvmkkavayLePylekekeedk 731 
                                                dl++a+ +  +pleii+gpL++Gm  vG LF + ++ +  v++sa+ mk+av++LeP+l k++  ++
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  651 DLDAALTRGDKPLEIINGPLMAGMAEVGRLFNDNQLIVAEVLQSAEAMKAAVSHLEPHLSKDE--TS 715 
                                                ***********************************************************9887..99 PP

                                 TIGR02082  732 skGkivlatvkGDvhDiGknivdvvLscngyevvdlGvkvPvekileaakkkkaDviglsGLivksl 798 
                                                skGk++latvkGDvhDiGkn+v+++Ls+ng+ev++lG+kv  e+++ aakk+++D iglsGL+vks 
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  716 SKGKMLLATVKGDVHDIGKNLVEIILSNNGFEVINLGIKVGPEELIAAAKKENPDFIGLSGLLVKSA 782 
                                                ******************************************************************* PP

                                 TIGR02082  799 demvevaeemerrgvkiPlllGGaalskahvavkiaekYkgevvyvkdaseavkvvdkllsekkkae 865 
                                                 +m+  a +++  g++ Pll+GGaals+a+ ++ i ++Y+g+v+y+kda+++++++++l++   +++
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  783 LQMIVTAADLKAAGIDAPLLVGGAALSRAFADTRITPEYNGPVLYAKDAMAGLELANQLVDPALRQQ 849 
                                                *************************************************************998877 PP

                                 TIGR02082  866 elekikeeyee.irekfgekkeklialsekaarkevfaldrsedlevpapkflGtkvleas.ieell 930 
                                                 +  + +++e  ++  + +    + ++        + +   s+++ + ++++l +++l+++ + +++
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  850 LMLDLARQQEAsAKIAAAK-AAGQSTV-----ATGSTHSAISSNAPILVAPDLEQHILRDIpVGQII 910 
                                                7655544433303333322.2222222.....2233344555578899999**************** PP

                                 TIGR02082  931 kyiDwkalFv.qWelrgkypkilkdeleglearklfkdakelldklsaekllrargvvGlfPaqsvg 996 
                                                +y++ ++l++ + +l g+ +k l    ++++a+kl+ +++++l+++ +e+l+++++++ ++ a+  g
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  911 PYLNRQMLYTkHLGLTGSVDKLLAG--QDEKATKLHLTVEAMLERVLQEGLIKPQAIYRFYQANGDG 975 
                                                **********899*******88875..67899*********************************** PP

                                 TIGR02082  997 ddieiytdetvsqetkpiatvrekleqlrqqsdrylclaDfiaskesGikDylgallvtaglgaeel 1063
                                                +d+ ++++++        a   +++  +rq+s+++lc+aDf+++ +  +kD +++++vt g+g++el
  lcl|NCBI__GCF_000020385.1:WP_012470188.1  976 NDLILFNQDGS------EA---TRFTLPRQKSGEQLCVADFVRPLSGTEKDTMALFAVTCGQGVREL 1033
                                                ******96665......33...35899*****************9999******************* PP

                                 TIGR02082 1064 akkleakeddydsilvkaladrlaealaellhervRkelwgyaeeenldkedllkerYrGirpafGY 1130
                                                +++++a+ d+ +s l++ala  +aea ae+lh+r+R   wg  +++ l +++l+  +Y+Gir +fGY
  lcl|NCBI__GCF_000020385.1:WP_012470188.1 1034 SEQWKAEGDYLNSHLLQALALEMAEATAEYLHKRIRTS-WGIVDDQTLTMKQLFNAEYQGIRVSFGY 1099
                                                ************************************87.666699********************** PP

                                 TIGR02082 1131 pacPdhtekatlleLleaeriGlklteslalaPeasvsglyfahpeakYfa 1181
                                                pacP+++++++l++Ll++e+iG++lte  ++ Peasvs+l+f+hpe kYf 
  lcl|NCBI__GCF_000020385.1:WP_012470188.1 1100 PACPNLDDQKKLFSLLKPEQIGINLTEGDMMDPEASVSALVFHHPEGKYFD 1150
                                                **************************************************6 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                            1  (1182 nodes)
Target sequences:                          1  (1169 residues searched)
Passed MSV filter:                         1  (1); expected 0.0 (0.02)
Passed bias filter:                        1  (1); expected 0.0 (0.02)
Passed Vit filter:                         1  (1); expected 0.0 (0.001)
Passed Fwd filter:                         1  (1); expected 0.0 (1e-05)
Initial search space (Z):                  1  [actual number of targets]
Domain search space  (domZ):               1  [number of targets reported over threshold]
# CPU time: 0.07u 0.03s 00:00:00.10 Elapsed: 00:00:00.10
# Mc/sec: 13.57
//
[ok]

This GapMind analysis is from Apr 10 2024. The underlying query database was built on Apr 09 2024.

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:

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:

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:

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