GapMind for Amino acid biosynthesis

 

Alignments for a candidate for metH in Acidovorax sp. GW101-3H11

Align methionine synthase; EC 2.1.1.13 (characterized)
to candidate Ac3H11_3274 5-methyltetrahydrofolate--homocysteine methyltransferase (EC 2.1.1.13)

Query= CharProtDB::CH_090726
         (1227 letters)



>FitnessBrowser__acidovorax_3H11:Ac3H11_3274
          Length = 903

 Score = 1052 bits (2721), Expect = 0.0
 Identities = 566/908 (62%), Positives = 666/908 (73%), Gaps = 32/908 (3%)

Query: 341  RLSGLEPLNIGEDSLFVNVGERTNVTGSAKFKRLIKEEKYSEALDVARQQVENGAQIIDI 400
            +LSGLEP++IGE +LFVN+GERTNVTGS  F R+I   +Y EAL VARQQVENGAQ+IDI
Sbjct: 2    KLSGLEPVSIGEGTLFVNIGERTNVTGSKAFARMILNGQYEEALAVARQQVENGAQVIDI 61

Query: 401  NMDEGMLDAEAAMVRFLNLIAGEPDIARVPIMIDSSKWDVIEKGLKCIQGKGIVNSISMK 460
            NMDE MLD++AAMVRFL LIA EPDIARVPIM+DSSKWDVIE GL+CIQGKGIVNSISMK
Sbjct: 62   NMDEAMLDSKAAMVRFLQLIASEPDIARVPIMVDSSKWDVIEAGLRCIQGKGIVNSISMK 121

Query: 461  EGVDAFIHHAKLLRRYGAAVVVMAFDEQGQADTRARKIEICRRAYKILTEEVGFPPEDII 520
            EGV+ F H A+L+RRYGAA VVMAFDE GQADT ARKIEIC RAY++L +EVGFPPEDII
Sbjct: 122  EGVEKFKHEARLVRRYGAAAVVMAFDEVGQADTYARKIEICERAYRVLVDEVGFPPEDII 181

Query: 521  FDPNIFAVATGIEEHNNYAQDFIGACEDIKRELPHALISGGVSNVSFSFRGNDPVREAIH 580
            FDPNIFAVATGIEEHNNYA DFI A   IK+ LP A +SGGVSNVSFSFRGNDPVREAIH
Sbjct: 182  FDPNIFAVATGIEEHNNYAVDFIEATRWIKQNLPGAKVSGGVSNVSFSFRGNDPVREAIH 241

Query: 581  AVFLYYAIRNGMDMGIVNAGQLAIYDDLPAELRDAVEDVILNRRDDGTERLLELAEKYRG 640
             VFLY+AI+ GMDMGIVNAG + +YDDL   LR+ VEDV+LNRR D  ERL+E+AE  + 
Sbjct: 242  TVFLYHAIQAGMDMGIVNAGMVGVYDDLEPVLRERVEDVVLNRRADAGERLVEIAETAKS 301

Query: 641  SKTDDTANAQQAEWRSW-----EVNKRLEYSLVKGITEFIEQDTEEARQ----QATRPIE 691
               D+   +++ EWR        V +RL ++LV GIT+FI +DTEEA Q    +  RP+ 
Sbjct: 302  GAKDE---SKKLEWRGTPEHPKTVGERLSHALVHGITDFIVEDTEEAYQAIVAKGGRPLH 358

Query: 692  VIEGPLMDGMNVVGDLFGEGKMFLPQVVKSARVMKQAVAYLEPFIEASKEQG-------K 744
            VIEGPLMDGMNVVGDLFG GKMFLPQVVKSARVMK AVA+L P+IE  K Q        +
Sbjct: 359  VIEGPLMDGMNVVGDLFGAGKMFLPQVVKSARVMKSAVAHLIPYIEEEKRQDELAGRDVR 418

Query: 745  TNGKMVIATVKGDVHDIGKNIVGVVLQCNNYEIVDLGVMVPAEKILRTAKEVNADLIGLS 804
            + GK++IATVKGDVHDIGKNIV VVLQCNN+E+V++GVMVP  +IL  AK   AD++GLS
Sbjct: 419  SKGKIIIATVKGDVHDIGKNIVTVVLQCNNFEVVNMGVMVPCHEILARAKVEGADIVGLS 478

Query: 805  GLITPSLDEMVNVAKEMERQGF----TIPLLIGGATTSKAHTAVKIEQNYSGPTVYVQNA 860
            GLITPSL+EM  VA EM++        IPLLIGGATTS+ HTAVKI  +Y GP VYV +A
Sbjct: 479  GLITPSLEEMQYVAGEMQKDDHFRIKKIPLLIGGATTSRVHTAVKIAPHYEGPVVYVPDA 538

Query: 861  SRTVGVVAALLSDTQRDDFVARTRKEYETVRIQHGRKKPRTPPVTLEAARDNDFAFDWQA 920
            SR+V V  +LL D  + ++V     +YE VR QH  KK +TP  +LE AR N    ++  
Sbjct: 539  SRSVSVAQSLLGDGVQ-NYVQELNADYEKVRTQHANKK-QTPMWSLEKARANRTPVNFST 596

Query: 921  YTPPVAHRLGVQEVE-ASIETLRNYIDWTPFFMTWSLAGKYPRILEDEVVGVEAQRLFKD 979
            Y PP    LG +  +   +  L  YIDW PFF TW LAG YP IL DEVVGVEA R+F D
Sbjct: 597  YVPPTPRLLGRRVFKNFDLAELAKYIDWGPFFQTWDLAGPYPAILTDEVVGVEATRVFAD 656

Query: 980  ANDMLDKLSAEKTLNPRGVVGLFPANRVGDDIEIYRDETRTHVINVSHHLRQQTEKTGF- 1038
               ML KL   + L+  GV+ L+PAN VGDDIE Y DETRT V    + +RQQTEK    
Sbjct: 657  GQAMLKKLIEGRWLSASGVMALYPANSVGDDIEFYTDETRTQVAMTWYGMRQQTEKQMID 716

Query: 1039 ----ANYCLADFVAPKLSGKADYIGAFAVTGGLEEDALADAFEAQHDDYNKIMVKALADR 1094
                 + CLADFVAPK SG ADY G FAVT GL  +    AF    DDY+ IM+K+LADR
Sbjct: 717  GVMRPSRCLADFVAPKDSGIADYAGLFAVTAGLGVEKKEQAFVDALDDYSAIMLKSLADR 776

Query: 1095 LAEAFAEYLHERVRKVYWGYAPNENLSNEELIRENYQGIRPAPGYPACPEHTEKATIWEL 1154
            LAEAFAE LH+RVR   WGYA  E LSN+E+I E Y GIRPAPGYPACP+H+ K  ++ +
Sbjct: 777  LAEAFAECLHQRVRTDLWGYAAGEQLSNDEMIAEKYAGIRPAPGYPACPDHSAKQELFRV 836

Query: 1155 LEVEKHTGMKLTESFAMWPGASVSGWYFSHPDSKYYAVAQIQRDQVEDYARRKGMSVTEV 1214
            L+ E+  GM LTES AM P ASVSG+Y  HPDS Y+ V +I  DQ+ D A R+GM    +
Sbjct: 837  LQCEE-IGMGLTESLAMTPAASVSGFYIGHPDSVYFNVGKIGEDQLHDMAERRGMDEKVL 895

Query: 1215 ERWLAPNL 1222
            ER LAPNL
Sbjct: 896  ERLLAPNL 903


Lambda     K      H
   0.318    0.134    0.391 

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: 2674
Number of extensions: 118
Number of successful extensions: 11
Number of sequences better than 1.0e-02: 1
Number of HSP's gapped: 2
Number of HSP's successfully gapped: 1
Length of query: 1227
Length of database: 903
Length adjustment: 45
Effective length of query: 1182
Effective length of database: 858
Effective search space:  1014156
Effective search space used:  1014156
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 methionine synthase; EC 2.1.1.13 (characterized)
to candidate Ac3H11_3297 5-methyltetrahydrofolate--homocysteine methyltransferase (EC 2.1.1.13)

Query= CharProtDB::CH_090726
         (1227 letters)



>FitnessBrowser__acidovorax_3H11:Ac3H11_3297
          Length = 354

 Score =  371 bits (953), Expect = e-107
 Identities = 186/344 (54%), Positives = 253/344 (73%), Gaps = 15/344 (4%)

Query: 3   SKVEQLRAQLNERILVLDGGMGTMIQSYRLNEADFRGE----------RFADWPCDLKGN 52
           ++  QL   L +RI++LDG MGTMIQ ++L EA +RGE          RF D+  D+KGN
Sbjct: 7   TRAAQLPDILAQRIVILDGAMGTMIQRFKLGEAQYRGEGYTGPGKAGDRFKDFAYDVKGN 66

Query: 53  NDLLVLSKPEVIAAIHNAYFEAGADIIETNTFNSTTIAMADYQMESLSAEINFAAAKLAR 112
           N+LL L++P+VI+ IH  Y  AGAD+IETNTF +TTIA  DY M +L+ E+N  +A+LAR
Sbjct: 67  NELLSLTRPDVISDIHERYLAAGADLIETNTFGATTIAQEDYHMANLAREMNLRSAELAR 126

Query: 113 ACADEWTARTPEKPRYVAGVLGPTNRTASISPDVNDPAFRNITFDGLVAAYRESTKALVE 172
           A  D+++  TP+KPR+VAG LGPT +TASISPDVNDP  RN+ F+ L AAY E T+ALVE
Sbjct: 127 AACDKYS--TPDKPRFVAGALGPTPKTASISPDVNDPGARNVDFEALRAAYYEQTEALVE 184

Query: 173 GGADLILIETVFDTLNAKAAVFAVKTEFEALGVELPIMISGTITDASGRTLSGQTTEAFY 232
           GG+D++L+ET+FDTLNAKAA+FA+   FE  G  LP++ISGT+TDASGR LSGQT  AF+
Sbjct: 185 GGSDVLLVETIFDTLNAKAALFAIDEYFENSGQRLPLIISGTVTDASGRILSGQTVTAFW 244

Query: 233 NSLRHAEALTFGLNCALGPDELRQYVQELSRIAE-CYVTAHPNAGLPNAFGE--YDLDAD 289
           +S+RH+  L  GLNCALG   +R Y+QEL+R+AE  +++ +PNAGLPN   +  +D   +
Sbjct: 245 HSVRHSRPLAIGLNCALGATLMRPYIQELNRVAEDTFISCYPNAGLPNPMSDTGFDETPE 304

Query: 290 TMAKQIREWAQAGFLNIVGGCCGTTPQHIAAMSRAVEGLAPRKL 333
             ++ + E+A  G +NI+GGCCGTTP HI A+++AVE +  RK+
Sbjct: 305 ITSRLVHEFAAEGLVNILGGCCGTTPDHIGAIAKAVESVPTRKM 348


Lambda     K      H
   0.318    0.134    0.391 

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: 1007
Number of extensions: 32
Number of successful extensions: 5
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: 1227
Length of database: 354
Length adjustment: 38
Effective length of query: 1189
Effective length of database: 316
Effective search space:   375724
Effective search space used:   375724
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: 54 (25.4 bits)

Align candidate Ac3H11_3274 (5-methyltetrahydrofolate--homocysteine methyltransferase (EC 2.1.1.13))
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.26825.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 1191.4   0.0          0 1191.2   0.0    1.0  1  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  5-methyltetrahydrofolate--homocy


Domain annotation for each sequence (and alignments):
>> lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  5-methyltetrahydrofolate--homocysteine methyltransferase (EC 2.1.1.1
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ! 1191.2   0.0         0         0     328    1182 .]       2     873 ..       1     873 [. 0.96

  Alignments for each domain:
  == domain 1  score: 1191.2 bits;  conditional E-value: 0
                                        TIGR02082  328 vlsglealkiaqessfvniGeRtnvaGskkfrklikaedyeealkiakqqveeGaqilDi 387 
                                                        lsgle+++i++ + fvniGeRtnv+Gsk f+++i +++yeeal +a+qqve+Gaq++Di
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274    2 KLSGLEPVSIGEGTLFVNIGERTNVTGSKAFARMILNGQYEEALAVARQQVENGAQVIDI 61  
                                                       68********************************************************** PP

                                        TIGR02082  388 nvDevllDgeadmkkllsllasepdiakvPlmlDssefevleaGLkviqGkaivnsislk 447 
                                                       n+De++lD++a+m+++l+l+asepdia+vP+m+Dss+++v+eaGL++iqGk+ivnsis+k
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274   62 NMDEAMLDSKAAMVRFLQLIASEPDIARVPIMVDSSKWDVIEAGLRCIQGKGIVNSISMK 121 
                                                       ************************************************************ PP

                                        TIGR02082  448 dGeerFlekaklikeyGaavvvmafDeeGqartadkkieiakRayklltekvgfppedii 507 
                                                       +G+e+F ++a+l+++yGaa+vvmafDe Gqa+t+++kiei++Ray++l+++vgfppedii
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  122 EGVEKFKHEARLVRRYGAAAVVMAFDEVGQADTYARKIEICERAYRVLVDEVGFPPEDII 181 
                                                       ************************************************************ PP

                                        TIGR02082  508 fDpniltiatGieehdryaidfieaireikeelPdakisgGvsnvsFslrgndavRealh 567 
                                                       fDpni+++atGieeh++ya+dfiea+r+ik++lP ak+sgGvsnvsFs+rgnd+vRea+h
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  182 FDPNIFAVATGIEEHNNYAVDFIEATRWIKQNLPGAKVSGGVSNVSFSFRGNDPVREAIH 241 
                                                       ************************************************************ PP

                                        TIGR02082  568 svFLyeaikaGlDmgivnagklavyddidkelrevvedlildrrreatekLlelaelykg 627 
                                                       +vFLy+ai+aG+Dmgivnag++ vydd+++ lre ved++l+rr +a e+L+e+ae+ k 
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  242 TVFLYHAIQAGMDMGIVNAGMVGVYDDLEPVLRERVEDVVLNRRADAGERLVEIAETAKS 301 
                                                       *********************************************************776 PP

                                        TIGR02082  628 tkeksskeaqeaewrnlp.....veeRLeralvkGeregieedleear....kklkaple 678 
                                                         +++   +++ ewr++p     v eRL++alv+G++++i ed+eea+     k  +pl+
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  302 GAKDE---SKKLEWRGTPehpktVGERLSHALVHGITDFIVEDTEEAYqaivAKGGRPLH 358 
                                                       66664...44666666543333389***********************444446689*** PP

                                        TIGR02082  679 iiegpLldGmkvvGdLFGsGkmfLPqvvksarvmkkavayLePylekekeed.......k 731 
                                                       +iegpL+dGm+vvGdLFG+GkmfLPqvvksarvmk ava+L+Py+e+ek +d       +
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  359 VIEGPLMDGMNVVGDLFGAGKMFLPQVVKSARVMKSAVAHLIPYIEEEKRQDelagrdvR 418 
                                                       *************************************************9999******* PP

                                        TIGR02082  732 skGkivlatvkGDvhDiGknivdvvLscngyevvdlGvkvPvekileaakkkkaDvigls 791 
                                                       skGki++atvkGDvhDiGkniv vvL+cn++evv++Gv+vP+++il  ak + aD++gls
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  419 SKGKIIIATVKGDVHDIGKNIVTVVLQCNNFEVVNMGVMVPCHEILARAKVEGADIVGLS 478 
                                                       ************************************************************ PP

                                        TIGR02082  792 GLivksldemvevaeemerrg.vk...iPlllGGaalskahvavkiaekYkgevvyvkda 847 
                                                       GLi++sl+em++va em++   ++   iPll+GGa++s+ h+avkia++Y+g+vvyv da
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  479 GLITPSLEEMQYVAGEMQKDDhFRikkIPLLIGGATTSRVHTAVKIAPHYEGPVVYVPDA 538 
                                                       *****************9975333446********************************* PP

                                        TIGR02082  848 seavkvvdkllsekkkaeelekikeeyeeirekfgekkeklialsekaarkevfaldrse 907 
                                                       s++v+v+++ll +   ++++++++++ye++r+++ + k+++   s+++ar ++  ++ s 
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  539 SRSVSVAQSLLGDGV-QNYVQELNADYEKVRTQHAN-KKQTPMWSLEKARANRTPVNFS- 595 
                                                       *************96.89****************98.7789999**********99999. PP

                                        TIGR02082  908 dlevpapkflGtkvleas.ieellkyiDwkalFvqWelrgkypkilkdeleglearklfk 966 
                                                        + +p+p+ lG++v+++  ++el kyiDw ++F +W+l+g yp il+de++g+ea+++f+
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  596 TYVPPTPRLLGRRVFKNFdLAELAKYIDWGPFFQTWDLAGPYPAILTDEVVGVEATRVFA 655 
                                                       9*********************************************************** PP

                                        TIGR02082  967 dakelldklsaekllrargvvGlfPaqsvgddieiytdetvsqetkpiatvrekleqlrq 1026
                                                       d +++l+kl++ + l a+gv+ l+Pa+svgddie ytdet++   ++  t +  ++q+++
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  656 DGQAMLKKLIEGRWLSASGVMALYPANSVGDDIEFYTDETRT---QVAMTWYGMRQQTEK 712 
                                                       ****************************************94...445555555555555 PP

                                        TIGR02082 1027 qsdr.....ylclaDfiaskesGikDylgallvtaglgaeelakkleakeddydsilvka 1081
                                                       q          claDf+a+k+sGi+Dy g+++vtaglg+e+  +++ +  ddy++i++k+
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  713 QMIDgvmrpSRCLADFVAPKDSGIADYAGLFAVTAGLGVEKKEQAFVDALDDYSAIMLKS 772 
                                                       54444667899************************************************* PP

                                        TIGR02082 1082 ladrlaealaellhervRkelwgyaeeenldkedllkerYrGirpafGYpacPdhtekat 1141
                                                       ladrlaea+ae lh+rvR +lwgya+ e+l++++++ e+Y Girpa+GYpacPdh+ k++
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  773 LADRLAEAFAECLHQRVRTDLWGYAAGEQLSNDEMIAEKYAGIRPAPGYPACPDHSAKQE 832 
                                                       ************************************************************ PP

                                        TIGR02082 1142 lleLleaeriGlklteslalaPeasvsglyfahpeakYfav 1182
                                                       l+++l+ e+iG+ ltesla++P+asvsg+y+ hp++ Yf+v
  lcl|FitnessBrowser__acidovorax_3H11:Ac3H11_3274  833 LFRVLQCEEIGMGLTESLAMTPAASVSGFYIGHPDSVYFNV 873 
                                                       ***************************************98 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                            1  (1182 nodes)
Target sequences:                          1  (903 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.08u 0.03s 00:00:00.11 Elapsed: 00:00:00.09
# Mc/sec: 11.27
//
[ok]

This GapMind analysis is from Apr 09 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