Computational sequence analysis of matrix metalloproteinases

Matrix metalloproteinases (MMP) play a cardinal role in the breakdown of extracellular matrix involved in a variety of biological and pathological processes. Research on MMPs has classified and characterized these enzymes according to their matrix substrate specificity, gene and protein domain structure, and regulation of activity and expression. However, the discovery of new MMPs has introduced a need for a more comprehensive and systematic method of classification and quantitative comparison of known and newly discovered members. This study compiles a sequence alignment, constructs a dendrogram, and calculates physical data and homology percentage assignments in order to obtain further insight into MMP structure-function relationships. Thorough analysis of MMP primary sequence domains, physical data patterns, and statistical analysis of sequence homology yields higher resolution in the similarities and differences that group MMP members.     

Structure‐guided approach for detecting large domain inserts in protein sequences as illustrated using the haloacid dehalogenase superfamily

In multi‐domain proteins, the domains typically run end‐to‐end, that is, one domain follows the C‐terminus of another domain. However, approximately 10% of multi‐domain proteins are formed by insertion of one domain sequence into that of another domain. Detecting such insertions within protein sequences is a fundamental challenge in structural biology. The haloacid dehalogenase superfamily (HADSF) serves as a challenging model system wherein a variable cap domain (～5–200 residues in length) accessorizes the ubiquitous Rossmann‐fold core domain, with variations in insertion site and topology corresponding to different classes of cap types. Herein, we describe a comprehensive computational strategy, CapPredictor, for determining large, variable domain insertions in protein sequences. Using a novel sequence‐alignment algorithm in conjunction with a structure‐guided sequence profile from 154 core‐domain‐only structures, more than 40,000 HADSF member sequences were assigned cap types. The resulting data set afforded insight into HADSF evolution. Notably, a similar distribution of cap‐type classes across different phyla was observed, indicating that all cap types existed in the last universal common ancestor. In addition, comparative analyses of the predicted cap‐type and functional assignments showed that different cap types carry out similar chemistries. Thus, while cap domains play a role in substrate recognition and chemical reactivity, cap‐type does not strictly define functional class. Through this example, we have shown that CapPredictor is an effective new tool for the study of form and function in protein families where domain insertion occurs. Proteins 2014; 82:1896–1906. © 2014 Wiley Periodicals, Inc.     

夏家店等古人骨DNA的提取、扩增及序列分析

古DNA从距今 2 0 0 0～ 50 0 0年的古代人骨 (编号为 :M89,ⅢM11,ⅢM12和AM4 )中提取出来 ,利用聚合酶链式反应分别对线粒体DNA中MTND4 基因中 112 10～ 11414(2 0 5bp)以及RegionV中 8196～ 83 16(12 1bP)进行了扩增 ,并被测序 .实验结果表明古代人骨中仍存在有大量的遗传信息 ,通过提取、扩增可以得到真实可靠的序列 .对序列进行比较发现 112 4 8,112 4 9,112 83及 112 93为易发生突变的位点 .结果经同源性分析后发现 ,ⅢM11和ⅢM12有很高同源性 ,而M89和AM4与ⅢM11和ⅢM12的同源性却比较低 .这对人类学、考古学及分子进化的研究有重要的辅助作用 .     

不同栽培条件辣椒疫霉菌ITS序列同源性分析

由疫霉菌引起的辣椒疫病是辣椒生产上的主要病害。明确设施栽培辣椒土传病害和露地栽培条件下辣椒疫病疫霉病菌基因序列和同源性,为制定病害防治措施提供依据。通过病菌形态观察、ITS基因序列测定与同源性分析开展研究。由疫霉菌侵染的辣椒病害在设施和露地不同栽培形式下发生部位略有不同,但在病菌形态上表现一致,ITS基因序列同源性极高,达到了99%以上。设施土传病害和露地栽培条件下辣椒疫病菌属于同一种,均为辣椒疫霉菌(Phytophthora capsici Leonian)。     

Exploration of the active site of Escherichia coli cystathionine γ‐synthase

Cystathionine γ‐synthase (CGS) catalyzes the condensation of O‐succinyl‐L ‐homoserine (L ‐OSHS) and L ‐cysteine (L ‐Cys), to produce L ‐cystathionine (L ‐Cth) and succinate, in the first step of the bacterial transsulfuration pathway. In the absence of L ‐Cys, the enzyme catalyzes the futile α,γ‐elimination of L ‐OSHS, yielding succinate, α‐ketobutyrate, and ammonia. A series of 16 site‐directed variants of Escherichia coli CGS (eCGS) was constructed to probe the roles of active‐site residues D45, Y46, R48, R49, Y101, R106, N227, E325, S326, and R361. The effects of these substitutions on the catalytic efficiency of the α,γ‐elimination reaction range from a reduction of only ～2‐fold for R49K and the E325A,Q variants to 310‐ and 760‐fold for R361K and R48K, respectively. A similar trend is observed for the k_cat/K of the physiological, α,γ‐replacement reaction. The results of this study suggest that the arginine residues at positions 48, 106 and 361 of eCGS, conserved in bacterial CGS sequences, tether the distal and α‐carboxylate moieties, respectively, of the L ‐OSHS substrate. In contrast, with the exception of the 13‐fold increase observed for R106A, the K is not markedly affected by the site‐directed replacement of the residues investigated. The decrease in k_cat observed for the S326A variant reflects the role of this residue in tethering the side chain of K198, the catalytic base. Although no structures exist of eCGS bound to active‐site ligands, the roles of individual residues is consistent with the structures inhibitor complexes of related enzymes. Substitution of D45, E325, or Y101 enables a minor transamination activity for the substrate L ‐Ala.     

我国新分离ECHO30病毒VP1序列分析

测定了引起2003年苏北地区无菌性脑膜炎暴发流行的病因病毒FDJS03分离株的VP1基因序列,并与国外同型流行毒株做比较,以了解本流行株的分子生物学特点及遗传变异规律。随机选取FDJS03分离毒株中的4株,用肠道病毒、VP1序列的特异性引物008／011进行RT-PCR,扩增产物经凝胶纯化后测序。将序列输入GenBank,用BLAgr程序进行核苷酸和氨基酸序列比对;选取32株不同地区不同年代的ECHO30分离毒株,在PHYLIP3．573C和TREE-PUZZLE5．0中构建进化树,比较它们完整VP1序列(876nt)的进化关系。核苷酸和氨基酸同源性比较结果证明：4株分离病毒均为ECHO30。进化树分析显示：本次FDJS03分离株与欧美20世纪70、80年代流行株亲缘关系最近,但自成一簇,与国外毒株仍然有地区差别。ECHO30的VP1基因进化有时间效应,但存在地区差异。本次流行的病原可能是单一基因型的ECHO30病毒。     

In silico structure–function analysis of E. cloacae nitroreductase

Reduction, catalyzed by the bacterial nitroreductases, is the quintessential first step in the biodegradation of a variety of nitroaromatic compounds from contaminated waters and soil. The Enterobacter cloacae nitroreductase (EcNR) enzyme is considered as a prospective biotechnological tool for bioremediation of hazardous nitroaromatic compounds. Using diverse computational methods, we obtain insights into the structural basis of activity and mechanism of its function. We have performed molecular dynamics simulation of EcNR in three different states (free EcNR in oxidized form, fully reduced EcNR with benzoate inhibitor and fully reduced EcNR with nitrobenzene) in explicit solvent and with full electrostatics. Principal Component Analysis (PCA) of the variance‐covariance matrix showed that the complexed nitroreductase becomes more flexible overall upon complexation, particularly helix H6, in the vicinity of the binding site. A multiple sequence alignment was also constructed in order to examine positional constraints on substitution in EcNR. Five regions which are highly conserved within the flavin mononucleotide (FMN) binding site were identified. Obtained results and their implications for EcNR functioning are discussed, and new plausible mechanism has been proposed. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

胡杨Na+/H+反向运输载体(PeNHX2)基因的克隆与序列分析

植物液泡膜上Na /H 反向运输载体(Na /H antiporterNHX)已被证明在盐胁迫中发挥关键作用,该研究旨在从胡杨中克隆获得耐盐相关的的完整cDNA序列,为进一步比较胡杨的耐盐基因的差异,探讨胡杨的耐盐分子机理奠定了基础,并为通过转基因手段提高木本植物的耐盐能力提供侯选基因。根据已发表NHX的同源基因序列,采用特异性引物扩增核心片段,并结合5’与3’RACE技术,成功地从胡杨(populus euphratica)克隆获得PeNHX2,其cDNA全长1638bp。测序和序列分析结果表明该基因与已发表毛白杨(Populus tomentosa)PtNHX基因在核酸水平及蛋白质水平上同源性均达到最高,依次为90.6%和88.25%。利用DNAMAN软件进一步分析得知,该基因与新疆盐生植物中所克隆获得的的犁苞滨藜(Atriplex dimorphostegia)AdNHX1、灰绿藜(Chenopodium glaucum)CgNHX1、盐爪爪(Kalidium foliatum)KfNHX同源性也较高,核酸水平上依次为75.1%、74.8%、74.3%,蛋白质水平依次为78.3%、79.0%、78.3%。PeNHX2与NCBI已注册的部分PeNHX1序列,在核酸或蛋白质水平同源性(72.34%和73.72%)均较低,分析它们可能同属于胡杨NHX基因家族,相关生物学功能有待于进一步研究。     

柞蚕热休克蛋白70基因的克隆和序列分析

HSP70蛋白是受热等因素刺激后而诱导产生的蛋白质,是热休克蛋白家族中最重要的一员。采用RT-PCR方法克隆了柞蚕(Antheraea pernyi)热休克蛋白70基因(HSP70)的ORF序列(GenBank登录号:GU945199),该片段的序列长度为1905bp。生物信息学分析表明,该序列共编码634个氨基酸,预测蛋白的等电点和分子量大小分别为5.62kD和69.5kD。具有HSP70的保守性结构特征,与天蚕(Antheraea yamamai)、家蚕(Bombyx mor)、甘蓝夜蛾(Mamestra brassicae)、棉铃虫(Heliothis viriplaca)、甜菜夜蛾(Spodoptera exigua)、烟草夜蛾(Manduca sexta1)、膜翅目寄生蜂(Cotesia rubecula)的同源性分别为95.7%、78.5%、76.1%、77.3%、76.6%、74.7%、65.9%。根据它们的一级结构构建了系统进化树,进一步确立了它们之间的亲缘关系。     

SPF鸡胚中内源性白血病病毒全基因序列鉴定与分析

摘要:以国内某3家SPF鸡场的SPF鸡胚成纤维细胞提取的基因组DNA为模板,参照已发表的序列,设计合成了4对检测内源性白血病病毒引物,分别检测gag基因、pol基因、env基因和LTR片段,结果显示4者检出阳性率很高(gag,29/46;pol,27/46;env,24/46;LTR,31/46).设计合成了8对引物,选取4片段检测均为阳性的样品之一,经PCR成功扩增出了8段连续的、相互部分重叠的目的DNA片段,分别连接入T载体进行克隆测序.用DNAstar软件对测序结果进行拼接,从一个鸡胚得到了内源性白血病病毒前病毒全基因组序列.比较分析发现,该序列env基因与已知的E亚群内源性病毒代表株env基因的核苷酸序列同源性在98.5%以上,全基因组序列同源性在99.1%以上,而与其他亚群代表株同源性相对较低,env基因同源性仅为56.3%～91.5%.     

Clinanthus microstephium,an Amaryllidaceae Species with Cholinesterase Inhibitor Alkaloids: Structure−Activity Analysis of Haemanthamine Skeleton Derivatives

Plants of the Amaryllidaceae family are well‐known (not only) for their ornamental value but also for the alkaloids that they produce. In this report, the first phytochemical study of Clinanthus genus was carried out. The chemical composition of alkaloid fractions from Clinanthus microstephium was analyzed by GC/MS and NMR. Seven known compounds belonging to three structural types of Amaryllidaceae alkaloids were identified. An epimeric mixture of a haemanthamine‐type compound (6‐hydroxymaritidine) was tested as an inhibitor against acetyl‐ and butyrylcholinesterase enzymes (AChE and BChE, respectively), two enzymes relevant in the treatment of Alzheimer's disease, with good results. Structure–activity relationships through molecular docking studies with this alkaloid and other structurally related compounds were discussed.     

中国狂犬病疫苗生产株CTN-1全基因序列分析

本文首次对我国现行狂犬病疫苗生产用毒株CTN-1进行全长基因组序列测定和分析,为CTN-1疫苗在我国实际应用中的优势提供理论基础。利用RT-PCR方法分段扩增CTN-1全基因组序列,随后将PCR产物克隆到T载体、测序、拼接,用MegAlign软件比较CTN-1全基因组序列与国内外狂犬病疫苗株和街毒株全基因组序列的同源性;再以糖蛋白(G)基因为模板,用ClustalX和MEGA4软件进行系统进化分析。测序结果表明CTN-1全基因组序列的长度为11925nt(GenBank登录号FJ959397),序列分析表明CTN-1为基因Ⅰ型;CTN-1株全基因组序列与国内外狂犬病疫苗株和街毒株全基因组之间的同源性是81.5%～93.4%;与美国蝙蝠分离株SHBRV18的同源性最低,仅为81.5%;与新近从中国国内分离的野毒株HN10株的同源性最高,达93.4%。系统进化分析结果表明,CTN-1与国内不同地区大多数分离的狂犬病街毒株聚类于同一组内;而我国另一疫苗株aG株与国外疫苗株如Flury、PM、PV、ERA、RC-HL和个别中国街毒株分在另一组内。G基因氨基酸对比也显示CTN-1株与国内大多数街毒株的同源性高于其他疫苗株,结果说明CTN-1株较其他疫苗株病毒更适合制备用于预防中国狂犬病的灭活疫苗。     

Structure–Function Relationship between Analogues of the Antibacterial Peptide Indolicidin. I. Synthesis and Biological Activity of Analogues with Increased Amphipathicity and Elevated Net Positive Charge of the Molecule

The antibacterial peptide indolicidin and a number of its analogues were obtained by solid phase synthesis. An optimized method of the synthesis using the Boc strategy was suggested. It was shown that the therapeutic index of indolicidin analogues increased with a decrease in the total positive charge of the molecule and its amphipathicity; i.e., the hemolytic activity of analogues within the range of concentrations examined was practically absent, while the antibacterial activity was preserved.     

Comparison of human solute carriers

Solute carriers are eukaryotic membrane proteins that control the uptake and efflux of solutes, including essential cellular compounds, environmental toxins, and therapeutic drugs. Solute carriers can share similar structural features despite weak sequence similarities. Identification of sequence relationships among solute carriers is needed to enhance our ability to model individual carriers and to elucidate the molecular mechanisms of their substrate specificity and transport. Here, we describe a comprehensive comparison of solute carriers. We link the proteins using sensitive profile–profile alignments and two classification approaches, including similarity networks. The clusters are analyzed in view of substrate type, transport mode, organism conservation, and tissue specificity. Solute carrier families with similar substrates generally cluster together, despite exhibiting relatively weak sequence similarities. In contrast, some families cluster together with no apparent reason, revealing unexplored relationships. We demonstrate computationally and experimentally the functional overlap between representative members of these families. Finally, we identify four putative solute carriers in the human genome. The solute carriers include a biomedically important group of membrane proteins that is diverse in sequence and structure. The proposed classification of solute carriers, combined with experiment, reveals new relationships among the individual families and identifies new solute carriers. The classification scheme will inform future attempts directed at modeling the structures of the solute carriers, a prerequisite for describing the substrate specificities of the individual families.     

Cloning and Sequence Analysis of the Gene Encoding (6-4)photolyase from Dunaliella salina

DNA photolyase can repair UV-induced DNA damage in a light-dependent manner. A cDNA of (6-4)photolyase from Dunaliella salina (GenBank accession number: AY845324) was cloned, sequenced and its amino acid sequence was deduced. The derived amino acid sequence showed high homology with other (6-4)photolyases and a predicted 3D model was constructed by homology modeling. Revisions requested 20 May 2005 and 18 August 2005; Revisions received 2 August 2005 and 28 November 2005     

Conservation analysis of class‐specific positions in cytochrome P450 monooxygenases: Functional and structural relevance

Cytochrome P450 monooxygenases (CYPs) constitute a ubiquitous, highly divergent protein family. Nevertheless, all CYPs share a common fold and conserved catalytic machinery. Based on the electron donor system, 10 classes of CYPs have been described, but most CYPs are members of class I accepting electrons from ferredoxin which is being reduced by FAD‐containing reductase, or class II accepting electrons from FAD‐ and FMN‐containing CPR‐type reductase. Because of the low sequence conservation inside the two classes, the conserved class‐specific positions are expected to be involved in aspects of electron transfer that are specific to the two types of reductases. In this work we present results from a conservation analysis of 16,732 CYP sequences derived from an updated version of the Cytochrome P450 Engineering Database (CYPED), using two class‐specific numbering schemes. While no position was conserved on the distal, substrate‐binding surface of the CYPs, several class‐specific residues were found on the proximal, reductase‐interacting surface; two class I‐specific residues that were negatively charged, and three class II‐specific residues that were aromatic or charged. The class‐specific conservation of glycine and proline residues in the cysteine pocket indicates that there are class‐specific differences in the flexibility of this element. Four heme‐interacting arginines were conserved differently in each class, and a class‐specific substitution of a heme‐interacting tyrosine by histidine was found, pointing to a link between heme stabilization and the reductase type. Proteins 2014; 82:491–504. © 2013 Wiley Periodicals, Inc.     

五步蛇蛇毒金属蛋白酶cDNA的克隆和序列分析

抽提五步蛇毒腺总RNA,通过反转录PCR(RT-PCR)扩增出五步蛇毒腺中一种低分子量金属蛋白酶(aculysinl)的cDNA,克隆到pGMT-vector并测定了全序列.推导其编码的蛋白质序列,发现aculysinl是以酶原形式合成的分泌蛋白,酶原包括信号肽、前肽、金属蛋白酶成熟肽和间隔肽4个部分.金属蛋白酶成熟肽与其它蛇毒金属蛋白酶相比,蛋白质一级结构具有一定的同源性,有一个保守的Zn2+结合位点:HEXXHXXGXXH.Aculysinl含有6个半胱氨酸,推测形成3对链内二硫键.五步蛇低分子量金属蛋白酶cDNA的克隆,为研究蛇毒金属蛋白酶结构与功能的关系,以及开发治疗血栓药物打下了良好的基础     

Probing the Ubiquinone Reduction Site in Bovine Mitochondrial Complex I Using a Series of Synthetic Ubiquinones and Inhibitors

Studies of the structure–activity relationships of ubiquinones and specific inhibitors are helpful to probe the structural and functional features of the ubiquinone reduction site of bovine heart mitochondrial complex I. Bulky exogenous short-chain ubiquinones serve as sufficient electron acceptors from the physiological ubiquinone reduction site of bovine complex I. This feature is in marked contrast to other respiratory enzymes such as mitochondrial complexes II and III. For various complex I inhibitors, including the most potent inhibitors, acetogenins, the essential structural factors that markedly affect the inhibitory potency are not necessarily obvious. Thus, the loose recognition by the enzyme of substrate and inhibitor structures may reflect the large cavitylike structure of the ubiquinone (or inhibitor) binding domain in the enzyme. On the other hand, several phenomena are difficult to explain by a simple one-catalytic site model for ubiquinone.     

Analysis and Organization of Protein Sequence Data: A Retrospective Spanning Four Decades

Protein sequence data are as useful and valuable today as was envisioned by pioneering sequencers and by the organizers of the first sequence database. Sequence analysis was first the province of specialists who developed search, comparison, and tree-building methods. Microcomputers, communication satellites, and the Internet have made these methods accessible to any scientist. The rapid increase in the data has driven a succession of changes in how databases are compiled, distributed, and accessed. Large public databases have become international collaborations. Although they need to develop still more efficient ways to accumulate, organize, annotate, and standardize huge amounts of data, inadequate support is available for such efforts. Thus there will be greater reliance on direct input from the scientific community. The World Wide Web is essential but not sufficient for integrated access to related databases.