Molecular and Cellular Microbiology (MCM) | Microbial Genetics, Physiology and Metabolism
Microbiol. Biotechnol. Lett. 2021; 49(3): 458-465
https://doi.org/10.48022/mbl.2104.04014
Juo Choi1,2, Taejoon Park1, Daun Kang2, Jeongju Lee2, Yungpil Kim2, Pilgoo Lee2, Gregory J. Y. Chung2, and Kyungyun Cho1*
1Department of Biotechnology, Hoseo University, Asan 31499, Republic of Korea 2MECOX CureMed Co., Seoul 06744, Republic of Korea
Correspondence to :
Kyungyun Cho, kycho@hoseo.edu
Tubulysins are a group of bioactive secondary metabolites from myxobacteria exhibiting strong anticancer activity against various cancer cell lines. In this study, we describe the identification of putative tubulysin biosynthetic gene clusters (tubA~tubF) in the genome sequences of two tubulysin-producing myxobacterial strains, Archangium gephyra MEHO_002 and MEHO_004. The inactivation of the putative tubulysin biosynthetic genes resulted in a tubulysin-production defect. The DNA sequences of the A. gephyra MEHO_002 and MEHO_004 tubulysin biosynthetic genes were 97% identical, and the amino acid sequences of the encoded proteins shared a similarity of 97−100%. The nucleotide sequences of the tubulysin biosynthetic gene clusters in MEHO_002 and MEHO_004 were 86% identical to that in Cystobacter sp. SBCb004 known as a tubulysin-producing myxobacterium, and the organization of the clusters was identical except for the lack of a tubZ gene in the clusters in MEHO_002 and MEHO_004. The amino acid sequences of the proteins encoded by each gene were 88-97% similar to those encoded by SBCb004, and the domain compositions of the proteins were also identical.
Keywords: Archangium gephyra, tubulysin, myxobacteria, secondary metabolite
포유류 세포에서 미세소관(microtubule)은 세포의 분열,성장, 이동에 결정적인 역할을 한다[1]. 미세소관은 튜불린(tubulin)으로 구성되어 있으며, 튜불린의 중합과 해체를 통해 기능을 수행한다. 이러한 중요한 기능으로 인해 미세소관은 항암제 개발의 주요 타겟이 되어왔다. 미세소관의 기능을저해하는 대표적인 항암물질로는 paclitaxel, vinblastine, epothilone 등이 있다[2].
Tubulysin은 점액세균에 의해 생산되는 이차대사 생리활성물질로 튜불린의 중합을 저해하여 미세소관의 해체를 초래함으로써 세포독성을 보이는데(Fig. 1) [3−5], paclitaxel, vinblastine, epothilone 보다도 훨씬 더 강력한 활성을 나타내는 것으로 보고되었다[6]. Tubulysin은 다양한 암세포주에대해 항암활성을 보이지만[7, 8] 강한 세포독성으로 인해 최근에는 항체약물접합체(antibody-drug conjugate)의 탑재 약물로 관심을 받고 있다[9−12].
Tubulysin은 다섯 개의 아미노산(L-pipecolic acid, L-isoleucine, L-valine, L-cysteine, L-phenylalanine 또는 L-tyrosine)과 두 단위의 acetate로 핵심구조가 이루어져 있는데, 이는 5개의 비리보솜 펩타이드 합성효소(non-ribosomal peptide synthetase, NRPS) 모듈과 2개의 폴리케타이드합성효소(polyketide synthase, PKS) 모듈에 의해 생합성된다[13]. Tubulysin은 점액세균
제조한 플라스미드는 전기천공을 통해 점액세균 균주에 도입하였으며, 50 μg/ml 가나마이신이 들어있는 배지에서 배양함으로써 상동재조합에 의해 플라스미드가
균주들을 1% Amberlite XAD-16 레진(Sigma, USA)을 넣은 CYS 배지에서 7일 동안 진탕배양한 뒤 레진을 회수하였다. 레진을 증류수로 3회 세척한 뒤, 레진 무게의 5배(v/w)에해당하는 100% 메탄올로 3회 추출하고, 회전증발기를 사용하여 건조하였다. 증류수와 에틸아세테이트를 1:1 비율(v/v)로 섞어 넣은 후, 잘 흔들어준 다음, 에틸아세테이트 층을 회수하여 회전증발기로 건조하고 100% 메탄올에 용해하였다.
Tubulysin의 분석을 위해서는 Zorbax SB-C18 칼럼 (4.6 × 150 mm, 5 µm)을 장착한 Agilent 1260 VL Infinity Series HPLC 시스템을 사용하였다(Agilent Technologies, USA). 이동상 A 용매로 0.1% trifluoroacetic acid (FTA)를함유한 증류수, 이동상 B 용매로는 0.1% FTA를 함유한 아세토니트릴을 사용하였다. 용매의 유속은 0.5 ml/min으로 0−5분 간 40% B 용매, 5−25분 간 40에서 45% B 용매 기울기, 25−30분 간 45에서 100% B 용매 기울기, 30−35분 간 100%B 용매, 35−40분 간 40% B 용매로 용리하였다. 용리액은photo diode array (PDA) 검출기를 사용하여 254 nm에서검출하였다.
Tubulysin의 분리를 위해서는 Agilent Zorbax SB-C18 PrepHT 칼럼(21.2 × 250 mm, 7 µm)을 사용하였다. 이동상A는 0.1% TFA를 함유한 증류수, 이동상 B는 0.1% TFA를함유한 아세토니트릴을 사용하였다. 유속은 6 ml/min으로0−60분 간 40−50% B 용매 기울기, 60−70분 간 100% B 용매, 70−71분 간 100−40% B 용매 기울기, 71−80분 간 40%B 용매 방법으로 분리를 하였다.
Tubulysin의 tandem MS 분석은 한국기초과학지원연구원(KBSI)에 의뢰하여 15T FT-ICR 질량분석기(solariX XRTM system, Bruker Daltonics, USA)를 사용하여 이루어졌다.
유전체 DNA의 분리는 일반적으로 알려진 방법을 사용하였다[17]. 유전체 DNA의 염기서열 결정은 마크로젠에 Whole Genome De novo Sequencing을 의뢰하여 실시하였다. 이차대사산물 생합성 유전자의 분석은 antiSMASH (antibiotics & Secondary Metabolite Analysis SHell) 프로그램[18]을이용하였다. 서열 상동성 분석은 미국 국립생물정보센터(NCBI)의 BLAST 프로그램[19]을 사용하였고, 단백질 도메인 분석은 NCBI의 CD-Search 프로그램[20]을 사용하였다.
Tubulysin 생합성 유전자군을 탐색하기 위하여
antiSMASH 프로그램은 생물의 유전체 정보로부터 다양한 종류의 이차대사산물 생합성 유전자군을 밝혀내며, 유전자 각각에 대한 상세 분석도 제공하는 프로그램이다[18]. antiSMASH 프로그램을 이용하여 MEHO_002 균주 유전체에 존재하는 이차대사산물 생합성 유전자군을 분석한 결과, tubulysin 생합성 유전자군으로 추정되는 유전자군이 길이가 가장 긴 contig에 존재하는 것으로 나타났다(Fig. 3A). MEHO_004 균주의 경우에도 마찬가지로 길이가 가장 긴contig에 tubulysin 생합성 유전자군으로 추정되는 유전자군이 존재하는 것으로 분석되었다(Fig. 3B). Tubulysin 생합성유전자군으로 추정되는 유전자군은 약 38 kb 크기로 8개 유전자로 구성되어 있는데(Fig. 3), 각 유전자가 암호화하는 단백질의 아미노산 서열이
사라진 피크가 tubulysin인지 확인하기 위하여 MEHO_ 004 균주의 배양추출물에서 tubulysin A에 해당하는 peak물질을 분리하여 15T FT-ICR 질량분석기를 사용하여 질량을 분석하였다. 1차 이온화한 물질은 [M+H]+의 m/z 값이844.452655로 0.19 ppm에서 분자식이 C43H66N5O10S로 추정되었으며, 이는 피크에 대한 동위원소 미세구조 분석을통해서 확인되었다. 따라서 해당 물질의 분자식은C43H65N5O10S로 tubulysin A일 것으로 예상되었다. 이 피크물질을 2차 이온화하여 분석하였을 때 m/z 값이 742.384353, 504.216237, 239.175383, 211.180471인 조각 피크들이 검출되었는데, 이러한 조각 양상은 기존에 알려진 tubulysin A의 조각 양상(MoNA ID, CCMSLIB00000478582)과 흡사하였다. 이에 더해, 조각피크의 분자량을 기초로 재구성한 물질의 구조도 tubulysin A인 것으로 분석되었다(Fig. 5). 따라서 이러한 결과는 플라스미드 삽입에 의해 사라진 피크물질은 tubulysin이며, 본 연구에서 발견된
Table 1 . Comparison between the tubulysin biosynthetic genes of
Identity/Similarity (%) | ||||
---|---|---|---|---|
Gene | Product size (aa) | Gene | Product size (aa) | |
431 | 96/99 | 430 | ||
219 | 99/100 | 219 | ||
394 | 99/100 | 394 | ||
1,498 | 98/98 | 1,498 | ||
2,622 | 98/98 | 2,622 | ||
3,509 | 98/98 | 3,509 | ||
1,159 | 98/98 | 1,159 | ||
2,831 | 97/97 | 2,827 | ||
367 | 96/98 | 367 | ||
120 | 95/95 | 120 |
Identity/Similarity: Identity and similarity of the amino acid sequences between the proteins encoded by
Table 2 . Comparison between the tubulysin biosynthetic genes of
Gene | Product size (aa) | Predicted function | PKS/NRPS motif predicted by antiSMASH | Identity/ Similarity (%) | Gene | Product size (aa) |
---|---|---|---|---|---|---|
431 | NRPS | C | 81/88 | 432 | ||
219 | hypothetical protein | 76/87 | 219 | |||
- | - | 386 | ||||
394 | hypothetical protein | 91/97 | 394 | |||
1,498 | NRPS | C-A(pip)-nMT-PCP | 84/89 | 1,520 | ||
2,622 | NRPS | C-A(ile)-PCP-C-A(val)-nMT-PCP | 85/90 | 2,626 | ||
3,509 | PKS/NRPS | KS-AT(mal)-KR-DH-ER-ACP-C-A(cys)-PCP | 85/90 | 3,511 | ||
1,159 | NRPS | C-CAL-PCP | 86/90 | 1,161 | ||
2,831 | PKS | KS-AT(mal)-KR-cMT-KR-DH-ER-ACP-TE | 82/88 | 2,843 | ||
- | - | 337 | ||||
367 | patatin-like protein | 80/89 | 368 | |||
120 | hypothetical protein | 60/72 | 124 |
A, adenylation domain; aa, amino acid; ACP, acyl-carrier protein domain; AT, acyltransferase; C, condensation domain; CAL, coenzyme A ligase domain; cMT, carbon methyltransferase; cys, cysteine; DH, dehydrase; ER, enoylreductase domain; ile, isoleucine; KR, ketoreductase domain; KS, ketosynthase domain; mal, malonate; nMT, nitrogen methyltransferase; NRPS, non-ribosomal peptide synthetase; PCP, peptidyl-carrier protein domain; pip, pipecolic acid; PKS, polyketide synthase; TE, thioesterase domain; val, valine.
Identity/Similarity: Identity and similarity to the corresponding proteins encoded by Cystobacter sp. SBCb004.
Table 3 . Comparison between the tubulysin biosynthetic genes of
Gene | Product size (aa) | Predicted function | PKS/NRPS motif predicted by antiSMASH | Identity/Similarity (%) | Gene | Product size (aa) |
---|---|---|---|---|---|---|
430 | NRPS | C | 81/89 | 432 | ||
219 | hypothetical protein | 77/87 | 219 | |||
- | - | 386 | ||||
394 | hypothetical protein | 91/97 | 394 | |||
1,498 | NRPS | C-A(pip)-nMT-PCP | 84/89 | 1,520 | ||
2,622 | NRPS | C-A(ile)-PCP-C-A(val)-nMT-PCP | 85/90 | 2,626 | ||
3,509 | PKS/NRPS | KS-AT(mal)-KR-DH-ER-ACP-C-A(cys)-PCP | 85/90 | 3,511 | ||
1,159 | NRPS | C-CAL-PCP | 86/90 | 1,161 | ||
2,827 | PKS | KS-AT(mal)-KR-cMT-KR-DH-ER-ACP-TE | 82/88 | 2,843 | ||
- | - | 337 | ||||
367 | patatin-like protein | 81/89 | 368 | |||
120 | hypothetical protein | 62/72 | 124 |
A, adenylation domain; aa, amino acid; ACP, acyl-carrier protein domain; AT, acyltransferase; C, condensation domain; CAL, co- enzyme A ligase domain; cMT, carbon methyltransferase; cys, cysteine; DH, dehydrase; ER, enoylreductase domain; ile, isoleucine; KR, ketoreductase domain; KS, ketosynthase domain; mal, malonate; nMT, nitrogen methyltransferase; NRPS, non-ribosomal peptide synthetase; PCP, peptidyl-carrier protein domain; pip, pipecolic acid; PKS, polyketide synthase; TE, thioesterase domain; val, valine.
Identity/Similarity: Identity and similarity to the corresponding proteins encoded by
MEHO_002와 MEHO_004 균주는 tubulysin 생합성 유전자군의 DNA 염기서열이 서로 97% 동일하며 유전체 서열도전반적으로 매우 유사하다. 하지만 MEHO_002 균주는argyrin 생합성 유전자군을 가져 argyrin을 생산하는 반면, MEHO_004 균주는 argyrin 생합성 유전자를 가지지 않아argyrin을 생산하지 않는다는 점이 다르다(Fig. 4). Argyrin은 점액세균과 방선균에서 분리된 cyclic octapeptide 계열의물질로 항암과 항세균 활성을 보인다[24, 25]. MEHO_002와 MEHO_004 균주의 tubulysin 생산수율은 유사한데, HPLC 분석 시 상대적으로 많은 양이 생산되는 argyrin이tubulysin과 유사한 위치에서 분리되기 때문에 argyrin을 생산하지 않는 MEHO_004 균주가 tubulysin을 생산하여 분리하고자 할 때 더 유리할 것으로 보인다.
결론적으로, 본 연구에서는 유전체 분석을 통해
Tubulysin은 다양한 암세포주에 대해 강한 항암활성을 보이는 점액세균 유래 이차대사 생리활성물질이다. 본 연구에서는 tubulysin을 생산하는 두 균주의 점액세균
The authors have no financial conflicts of interest to declare.
Hyesook Hyun, Juo Choi, Daun Kang, Yungpil Kim, Pilgoo Lee , Gregory Chung and Kyungyun Cho
Microbiol. Biotechnol. Lett. 2021; 49(1): 32-38 https://doi.org/10.48022/mbl.2010.10001Juo Choi, Taejoon Park, Daun Kang, Jeongju Lee, Yungpil Kim, Pilgoo Lee, Gregory J. Y. Chung, and Kyungyun Cho
Microbiol. Biotechnol. Lett. 2021; 49(4): 493-500 https://doi.org/10.48022/mbl.2107.07009Yujin Ka, Chaehyeon Park, Gangmin Kim, and Kyungyun Cho
Microbiol. Biotechnol. Lett. 2024; 52(4): 517-519 https://doi.org/10.48022/mbl.2411.11007