动物组织蛋白酶L-like基因家族的进化分析

目的:组织蛋白酶L-like家族是在溶酶体中发现的一类非常重要的半胱氨酸组织蛋白酶。其主要功能为催化各种蛋白质的水解,并通过水解蛋白质参与到许多的生理调节过程当中。根据序列比对分析和传统的功能分类,在动物中,组织蛋白酶L-like家族成员包括组织蛋白酶L、V、S、K、H和F。但是这些家族成员之间的进化关系仍然没有详细研究分析清楚。本课题主要研究组织蛋白酶L-like家族成员之间的进化关系。方法:本研究通过搜集整理22个物种的177条组织蛋白酶L-like家族蛋白的序列,并构建系统发育进化树来分析组织蛋白酶L-like家族各成员之间的进化关系。结果:序列数据结果显示,串联重复在组织蛋白酶L-like家族的进化过程中发生。斑马鱼的组织蛋白酶L,爪蟾的组织蛋白酶S和K,大鼠和小鼠的组织蛋白酶L都发生过明显的串联重复事件。进化树结果显示了组织蛋白酶H、S和K、L和V之间的进化关系,组织蛋白酶S和K在脊椎动物出现的进化过程中,从组织蛋白酶L中分化出来,与他们在脊椎动物体内的特异性功能,以及脊椎动物在进化过程中分化产生的特异性功能相对应。结论:在物种进化的过程中,组织蛋白酶L-like家族成员F、H、S和K、L和V按时间顺序分化,这表明组织蛋白酶L-like基因家族结构和功能的分化与新的物种和新的功能出现密切相关。     

半胱氨酸组织蛋白酶与免疫调节

半胱氨酸组织蛋白酶(CC)属于蛋白酶clanCA中的木瓜蛋白酶家族(C1),由组织蛋白酶B、L、H、S、K、F、V、X、W、O和C等组成。这些蛋白酶参与了各种生理和病理过程。     

脊椎动物Sox基因家族的系统发生分析

脊椎动物Sox基因是一类高度保守的基因家族,它们编码一类转录因子,参与多种发育过程的调控,这个家族的共有特征是Sox蛋白具有一个约79个氨基酸,可与DNA特异结合的HMG盒区,该基因家族成员众多复杂,对它们的基因结构,功能及进化关系的研究有助于对该基因家族的全面认识,利用已克隆的全部脊椎动物全长核苷酸／蛋白质数据,对这些数据进行序列比对及进化树的构建,分析Sox基因家族成员的分类及分子进化模式。     

RHR转录因子家族起源、功能以及进化机制的研究进展

转录因子是一类能够通过与基因特异性序列进行结合,从而调控基因转录与表达的蛋白质,对细胞的生物学活性具有重要的调节作用。RHR(Rel-homology region,RHR)转录因子家族属于IF(immunoglobulin fold)转录因子超家族最主要的成员,其成员含有保守的Rel结构域和IPT(immunoglobulin-like fold)结构域。作为古老的转录因子家族,RHR家族成员随着物种演化,通过基因的复制、突变和沉默,不断分化出新型同源基因的同时也伴随着基因的丢失。自然选择导致了各家族成员不同的进化速率,并且在一些功能结构域上展现出了特殊的进化机制。然而,目前有关RHR家族起源和分化的综述比较少见。本文综述了RHR家族各成员的分布、分类、功能及家族进化等方面的研究成果,以期为研究整个转录因子家族的演化机制和物种之间的进化关系提供参考和新的思路。     

七鳃鳗KCTD基因家族成员全基因组鉴定及进化分析

钾离子通道四聚化结构域(KCTD)蛋白基因家族是一个保守的基因家族,该家族成员的共同特征是具有一个含有BTB保守结构域的N-末端和一个可变的C-末端。KCTD基因的突变或不正常调控与人类多种疾病相关。七鳃鳗是现存最原始的脊椎动物,作为联系无脊椎动物和脊椎动物之间的桥梁,在生物进化研究中占有重要地位。本研究通过对海七鳃鳗(Petromyzon marinus)和日本七鳃鳗(Lethenteron japonicum)基因组和转录组数据分析,全面系统地鉴定了海七鳃鳗和日本七鳃鳗KCTD基因家族成员,并对其基因结构特征、蛋白保守基序和基因表达模式进行了分析。在海七鳃鳗和日本七鳃鳗中分别鉴定出13个和14个KCTD基因,基因长度和外显子数目在不同KCTD基因间变化很大,KCTD蛋白中4个基序保守性显著,大多数KCTD基因呈泛表达模式,并且在胚胎发育时期明显高表达。除七鳃鳗外,对12个无脊椎动物和脊椎动物代表物种KCTD基因家族成员进行了鉴定,并对KCTD基因家族成员的进化关系进行了分析。根据进化树聚类情况,将KCTD基因家族成员分为11个亚家族。进化分析结果显示,KCTD基因家族从低等的无脊椎动物线虫和果蝇到高等的人类都存在;线虫中仅有5个成员,果蝇中有8个成员,随着物种进化程度由低到高,KCTD家族成员数目呈现增加的趋势;从爬行类开始,脊椎动物KCTD基因数目稳定在24个左右。硬骨鱼类特有的全基因组复制事件影响鱼类KCTD基因数目。本研究结果不仅丰富了七鳃鳗KCTD基因家族信息,同时也对KCTD家族基因间的进化关系进行了探究,为深入研究该家族基因功能提供了一定的依据。     

烟草蛋白酶抑制剂基因的电子克隆及生物信息学分析

蛋白酶抑制剂可以增强植物对病虫害的抵抗能力,为了深入研究蛋白酶抑制剂在烟草中的作用机制,利用生物信息学的方法,成功获得了烟草品种K326中的4种蛋白酶抑制剂cDNA序列(NtPI-1、NtPI-2、NtPI-3和NtPI-4)并对其进行了序列分析.它们编码的氨基酸序列都具有典型的马铃薯蛋白酶抑制剂Ⅰ家族功能结构域,属于马铃薯蛋白酶抑制剂Ⅰ家族成员.序列分析表明,4种蛋白酶抑制剂基因cDNA序列均具有完整的开放读码框,依次编码128、95、94和72个氨基酸残基,它们的氨基酸序列一致性在31％-38％之间.系统树分析表明,烟草品种K326中的4种蛋白酶抑制剂分散在进化树的不同位置,形成不同亚群.     

水稻Trihelix转录因子家族全基因组分析及功能预测

Trihelix转录因子家族在植物生长发育以及响应逆境胁迫等方面发挥着重要作用,但目前基于水稻全基因组水平鉴定和分析该基因家族的研究尚未见相关报道。本文利用生物信息学方法在水稻基因组数据库中鉴定到Trihelix家族成员31个,序列聚类和功能结构域分析发现该家族均含有高度保守的、特征性的Trihelix结构域;根据亲缘关系远近和结构域特点,将其分为5个亚家族(Ⅰ~Ⅰ)。通过与拟南芥、二穗短炳草和高粱中Trihelix家族的聚类分析发现,这4个物种中Trihelix家族的分类相一致,但每个物种均含有不同亚家族的成员,表明该基因家族的分化早于物种的分化。基于MEME程序分析水稻Trihelix转录因子家族的保守基序与聚类分析结果具有较高的一致性。染色体区段复制分析表明,部分Trihelix家族成员在水稻以及水稻与其他物种之间存在种内和种间的染色体区段复制;生物芯片数据分析发现,Trihelix基因家族在水稻不同组织中、以及对6种不同植物激素的响应呈现多样化的表达谱。采用Rice FREND在线数据库分析发现,水稻Trihelix转录因子家族的20个成员与其他蛋白存在互作关系。本研究结果初步明确了水稻Trihelix转录因子家族的进化特点、染色体分布、染色体区段复制关系、组织表达、激素应答,以及该家族蛋白与其他蛋白质的互作情况,为进一步揭示Trihelix转录因子家族的分子进化规律和生物学功能奠定了基础。     

家蚕组织蛋白酶基因家族的鉴定及表达特征分析

家蚕是鳞翅目完全变态昆虫,在其变态过程中伴随着巨大的形态变化,包括旧组织的解离和新组织的形成,在这过程中有多种组织蛋白酶参与。组织蛋白酶是一类细胞内蛋白酶,广泛存在于各个物种中,包括组织蛋白酶B、H、L等几个亚家族。对家蚕组织蛋白酶的研究将有利于阐明家蚕变态发育的详细过程。通过对家蚕基因组数据库进行筛选,共在家蚕中鉴定到13种组织蛋白酶,并对这13种组织蛋白酶的基本信息和表达模式进行了分析。另外,利用家蚕基因芯片数据和荧光定量PCR分析,鉴定编号为BGIBMGA004622的基因为卵巢特异表达的组织蛋白酶L亚家族基因。该基因全长1 209 bp,编码402个氨基酸。经过序列分析,该酶与其他物种的组织蛋白酶L具有较高的同源性,其活性位点高度保守,且与鳞翅目的组织蛋白酶L在进化上聚为一支。同时,对该基因进行克隆并原核表达,结果显示重组蛋白以包涵体的形式表达。定量PCR结果显示,该酶在蛹发育初期表达量逐渐升高,至蛹3 d达到最高值,推测其可能参与卵巢与卵母细胞的发育过程。     

植物miR159家族成员分子特性及其进化规律研究

为了解植物miR159家族成员的分子特性及其进化规律,该研究对miRBase数据库中登录的miR159家族成员进行分类统计、进化树构建、科间比较、二级结构预测及靶基因分析。结果表明:miR159家族在植物界分布非常广泛,蕨类植物可能是miR159家族的进化祖先;系统发育进化树分析显示,植物miR159家族成员间存在多个进化分支,且进化关系与植物属性有关,即植物亲缘关系越近的成员更易成枝,且具有相同进化方向的成员序列高度同源;Mfold预测显示,pre-miR159均会自发形成典型、稳定的茎环二级结构,并包含19~21个碱基为单位的功能片段,每个单位均有可能形成miR159成熟体;靶基因分析发现,miR159家族成员主要作用于MYB转录因子、转座因子和假定蛋白等,但在不同物种间及相同物种的不同成员间作用的靶标种类与靶基因ID数量均存在差异,尤其是miR159-3p与miR159-5p间的差异最为显著。     

趋化因子及其受体基因家族的系统进化分析

通过分析现有的趋化因子和趋化因子受体的氨基酸序列,用距离法和最简约法构建了聚类图,探讨了趋化因子和趋化因子受体基因家族的系统演化特征。可见基因家族成员的分化早于脊椎动物的分化。不同物种的同一种基因的聚类关系能较好地反映物种经因子受体的进化速度不同,其中ＣＸＣＲ４的进化速率最低。趋化因子和趋化因子受体可能都起源于少数几个原始的基因,病毒编码与寄主相似的趋化因子或受体是进化过程中分子模拟的结果。     

Evolutionary History of Cathepsin L (L-like) Family Genes in Vertebrates

Cathepsin L family, an important cysteine protease found in lysosomes, is categorized into cathepsins B, F, H, K, L, S, and W in vertebrates. This categorization is based on their sequence alignment and traditional functional classification, but the evolutionary relationship of family members is unclear. This study determined the evolutionary relationship of cathepsin L family genes in vertebrates through phylogenetic construction. Results showed that cathepsins F, H, S and K, and L and V were chronologically diverged. Tandem-repeat duplication was found to occur in the evolutionary history of cathepsin L family. Cathepsin L in zebrafish, cathepsins S and K in xenopus, and cathepsin L in mice and rats underwent evident tandem-repeat events. Positive selection was detected in cathepsin L-like members in mice and rats, and amino acid sites under positive selection pressure were calculated. Most of these sites appeared at the connection of secondary structures, suggesting that the sites may slightly change spatial structure. Severe positive selection was also observed in cathepsin V (L2) of primates, indicating that this enzyme had some special functions. Our work provided a brief evolutionary history of cathepsin L family and differentiated cathepsins S and K from cathepsin L based on vertebrate appearance. Positive selection was the specific cause of differentiation of cathepsin L family genes, confirming that gene function variation after expansion events was related to interactions with the environment and adaptability.     

Ovarian cysteine proteinases in the teleost Fundulus heteroclitus: molecular cloning and gene expression during vitellogenesis and oocyte maturation

The cysteine proteinases cathepsins B and L are members of the multigene family of lysosomal proteases that have been implicated in the processing of yolk proteins (YPs) in teleost oocytes. However, the full identification of the type of cathepsins expressed in fish ovarian follicles and embryos, as well as their regulatory mechanisms and specific function(s), are not yet elucidated. In this study, cDNAs encoding cathepsins B, L, F, K, S, Z, C, and H have been isolated from the teleost Fundulus heteroclitus, and the analysis of their deduced amino acid sequences revealed highly similar structural features to vertebrate orthologs, and confirmed in this species the existence of cathepsin L-like, cathepsin B-like, and cathepsin F-like subfamilies of cysteine proteinases. While all identified cathepsins were expressed in ovarian follicles, the corresponding mRNAs showed different temporal expression patterns. Thus, similar mRNA levels of cathepsins L, F, S, B, C, and Z were found throughout the oocyte growth or vitellogenesis period, whereas those for cathepsin H and K appeared to decrease as vitellogenesis advanced. During oocyte maturation, a transient accumulation of cathepsins L, S, H, and F mRNAs, approximately a 3-, 1.5-, 1.6-, and 6-fold increase, respectively, was detected in ovarian follicles within the 20-25 hr after hormone stimulation, coincident with the maximum proteolysis of the oocyte major YPs. The specific temporal pattern of expression of these genes may indicate a potential role of cathepsin L-like and cathepsin F proteases in the YP processing events occurring during fish oocyte maturation and/or early embryogenesis.     

Human cathepsin V functional expression, tissue distribution, electrostatic surface potential, enzymatic characterization, and chromosomal localization

Cathepsin V, a thymus and testis-specific human cysteine protease, was expressed in Pichia pastoris, and its physicokinetic properties were determined. Recombinant procathepsin V is autocatalytically activated at acidic pH and is effectively inhibited by various cysteine protease class-specific inhibitors. The S2P2 subsite specificity of cathepsin V was found to be intermediate between those of cathepsins S and L. The substrate binding pocket, S2, accepted both aromatic and nonaromatic hydrophobic residues, whereas cathepsins L and S preferred either an aromatic or nonaromatic hydrophobic residue, respectively. In contrast to cathepsin L, but similar to cathepsin S, cathepsin V exhibited only a very weak collagenolytic activity. Furthermore, cathepsin V was determined to be significantly more stable at mildly acidic and neutral pH than cathepsin L, but distinctly less stable than cathepsin S. A homology structure model of cathepsin V revealed completely different electrostatic potentials on the molecular surface when compared with human cathepsin L. The model-based electrostatic potential of human cathepsin V was neutral to weakly positive at and in the vicinity of the active site cleft, whereas that of cathepsin L was negative over extended regions of the surface. Surprisingly, the electrostatic potential of the human cathepsin V model structure resembled that of the model structure of mouse cathepsin L. These differences in the electrostatic potential at the molecular surfaces provide a reactivity determinant that may be the source of differences in substrate selectivity and pH stability. Cathepsin V was mapped to the chromosomal region 9q22.2, a site adjacent to the cathepsin L locus. The high sequence identity and the overlapping chromosomal gene loci suggest that both proteases evolved from an ancestral cathepsin L-like precursor by gene duplication.     

Selective inhibition of the collagenolytic activity of human cathepsin K by altering its S2 subsite specificity

The primary specificity of papain-like cysteine proteases (family C1, clan CA) is determined by S2-P2 interactions. Despite the high amino acid sequence identities and structural similarities between cathepsins K and L, only cathepsin K is capable of cleaving interstitial collagens in their triple helical domains. To investigate this specificity, we have engineered the S2 pocket of human cathepsin K into a cathepsin L-like subsite. Using combinatorial fluorogenic substrate libraries, the P1-P4 substrate specificity of the cathepsin K variant, Tyr67Leu/Leu205Ala, was determined and compared with those of cathepsins K and L. The introduction of the double mutation into the S2 subsite of cathepsin K rendered the unique S2 binding preference of the protease for proline and leucine residues into a cathepsin L-like preference for bulky aromatic residues. Homology modeling and docking calculations supported the experimental findings. The cathepsin L-like S2 specificity of the mutant protein and the integrity of its catalytic site were confirmed by kinetic analysis of synthetic di- and tripeptide substrates as well as pH stability and pH activity profile studies. The loss of the ability to accept proline in the S2 binding pocket by the mutant protease completely abolished the collagenolytic activity of cathepsin K whereas its overall gelatinolytic activity remained unaffected. These results indicate that Tyr67 and Leu205 play a key role in the binding of proline residues in the S2 pocket of cathepsin K and are required for its unique collagenase activity.     

An evolutionarily conserved tripartite tryptophan motif stabilizes the prodomains of cathepsin L-like cysteine proteases.

Cathepsin L-like cysteine proteinases contain an evolutionarily highly conserved alpha-helical motif in the proregion. This is called the ER(F/W)N(I/V)N motif according to the conserved amino acids along one side of the helix. We studied the function of this motif using site-directed mutagenesis experiments of human procathepsin S. We replaced each of these amino acids with alanine and constructed deletion mutants lacking parts of the helix. All mutants were expressed in HEK 293 cells, but only one, W52A, was not processed to mature cathepsin S, nor was it phosphorylated or secreted into the culture medium. W52 is part of the hydrophobic core in the propeptide region of cathepsin S comprising two additional tryptophan residues, W28 and W31, also conserved among cathepsin L-like cysteine peptidases. Replacement of the latter with alanine led to consequences similar to those with the W52A mutation. Recombinant propeptides containing mutations of one of the three tryptophan residues were three orders of magnitude less effective as inhibitors of mature cathepsin S than the wild-type propeptide. The results point to a dominant role of the respective hydrophobic stack in the proper folding, transport and maturation of procathepsin S and related cathepsin L-like cysteine proteinases.     

Structural and functional conservation profiles of novel cathepsin L-like proteins identified in the Drosophila melanogaster genome

Cathepsin L is a cysteine protease which degrades connective tissue proteins including collagen, elastin, and fibronectin. In this study, five well-characterized cathepsin L proteins from different arthropods were used as query sequences for the Drosophila genome database. The search yielded 10 cathepsin L-like sequences, of which eight putatively represent novel cathepsin L-like proteins. To understand the phylogenetic relationship among these cathepsin L-like proteins, a phylogenetic tree was constructed based on their sequences. In addition, models of the tertiary structures of cathepsin L were constructed using homology modeling methods and subjected to molecular dynamics simulations to obtain reasonable structure to understand its dynamical behavior. Our findings demonstrate that all of the potential Drosophila cathepsin L-like proteins contain at least one cathepsin propeptide inhibitor domain. Multiple sequence alignment and homology models clearly highlight the conservation of active site residues, disulfide bonds, and amino acid residues critical for inhibitor binding. Furthermore, comparative modeling indicates that the sequence/structure/function profiles and active site architectures are conserved.     

Manipulating substrate and pH in zymography protocols selectively distinguishes cathepsins K, L, S, and V activity in cells and tissues

Cathepsins K, L, S, and V are cysteine proteases that have been implicated in tissue-destructive diseases such as atherosclerosis, tumor metastasis, and osteoporosis. Among these four cathepsins are the most powerful human collagenases and elastases, and they share 60% sequence homology. Proper quantification of mature, active cathepsins has been confounded by inhibitor and reporter substrate cross-reactivity, but is necessary to develop properly dosed therapeutic applications. Here, we detail a method of multiplex cathepsin zymography to detect and distinguish the activity of mature cathepsins K, L, S, and V by exploiting differences in individual cathepsin substrate preferences, pH effects, and electrophoretic mobility under non-reducing conditions. Specific identification of cathepsins K, L, S, and V in one cell/tissue extract was obtained with cathepsin K (37 kDa), V (35 kDa), S (25 kDa), and L (20 kDa) under non-reducing conditions. Cathepsin K activity disappeared and V remained when incubated at pH 4 instead of 6. Application of this antibody free, species independent, and medium-throughput method was demonstrated with primary human monocyte-derived macrophages and osteoclasts, endothelial cells stimulated with inflammatory cytokines, and normal and cancer lung tissues, which identified elevated cathepsin V in lung cancer.     

Species differences between human and rat in the substrate specificity of cathepsin K

Cathepsin K is known to play an important role in bone resorption, and it has the P2 specificity for proline. Rat cathepsin K has 88% identity with the human enzyme. However, it has been reported that its enzymatic activity for a Cbz-Leu-Arg-MCA substrate is lower than that of human cathepsin K, and that the rat enzyme is not well inhibited by human cathepsin K inhibitors. For this study, we prepared recombinant enzyme to investigate the substrate specificity of rat cathepsin K. Cleavage experiments using the fragment of type I collagen and peptidic libraries demonstrated that rat cathepsin K preferentially hydrolyses the substrates at the P2 Hyp position. Comparison of the S2 site between rat and human cathepsin K sequences indicated that two S2 residues at Ser134 and Val160 in rat are varied to Ala and Leu, respectively, in the human enzyme. Cleavage experiments using two single mutants, S134A and V160L, and one double mutant, S134A/V160L, of rat cathepsin K showed that all the rat mutants lost the P2 Hyp specificity. The information obtained from our comparative studies on rat and human cathepsin K should make a significant impact on developing specific inhibitors of human cathepsin K since rat is usually used as test species.     

The alpha1/2 helical backbone of the prodomains defines the intrinsic inhibitory specificity in the cathepsin L-like cysteine protease subfamily

Proregions of papain-like cysteine proteases are potent and often highly selective inhibitors of their parental enzymes. The molecular basis of their selectivity is poorly understood. For two closely related members of the cathepsin L-like subfamily we established strong selectivity differences. The propeptide of cathepsin S was observed to inhibit cathepsin L with a K(i) of 0.08 nM, yet cathepsin L propeptide inhibited cathepsin S only poorly. To identify the respective structural correlates we engineered chimeric propeptides and compared their inhibitory specificity with the wild-types. Specificity resided in the N-terminal part, strongly suggesting that the backbone of the prodomain was the underlying structure.