人P2X_1蛋白生物信息预测

利用生物信息学方法对配体门控离子通道受体人P2X_1蛋白的理化性质、信号肽、跨膜区、亲疏水性、二级结构及三级结构、蛋白质间的相互作用、GO注释等进行预测分析。结果显示:P2X_1蛋白是两次跨膜的亲水蛋白质,由399个氨基酸组成,等电点为8.75,有一段核定位序列;二级结构存在6个α螺旋区和15个β折叠区,三级结构预测结果的可靠性为45.35%,拉曼图分析表明预测结构较稳定;与P2X_1相互作用的蛋白质主要是传导信号的嘌呤受体和离子通道蛋白质,而且P2X_1蛋白可能参与血管凝集及炎症反应等生理过程。上述对人P2X_1蛋白结构及功能的预测为研究人P2X_1蛋白在生物过程及疾病治疗方面的作用提供了重要的理论依据。     

斑马鱼TATA结合蛋白的生物信息学分析

斑马鱼TATA结合蛋白(TBP)是转录过程中的重要起始因子。利用生物信息学方法对斑马鱼TBP的理化性质、物种间同源性、保守结构域、跨膜区、亲水性/疏水性、蛋白质二级结构、蛋白质三级结构、蛋白质相互作用进行预测分析。分析表明,斑马鱼TBP全长302个氨基酸,等电点9.8,属于TATA结合蛋白超家族,不含跨膜区,属于亲水蛋白;二级结构以无规则卷曲为主,含5个α螺旋区和8个β折叠区,三维建模空间结构可信度98.9%,进一步分析建模结果可靠;与斑马鱼TBP相互作用的蛋白质均为转录因子或TFⅡD复合物组分。分析结果对于深入研究斑马鱼TBP在基因转录中的作用具有一定的理论指导意义。     

牛NAD(+)异柠檬酸脱氢酶生物信息学分析

异柠檬酸脱氢酶(IDH)是三羧酸循环中的关键酶。为了进一步探索IDH的结构与功能,利用生物信息学方法对牛NAD(+)IDH进行分析。结果表明,该蛋白为亲水性蛋白,呈碱性,无跨膜区,在β亚基N端存在长度为16个氨基酸的信号肽,亚细胞定位在线粒体。二级结构预测结果显示该蛋白的三个亚基均以α螺旋为主,并且其上的螺旋和折叠均紧密有序排列,这有助于形成特定的结构域并发挥特定的生物学功能。对α和γ两个亚基采用同源建模法预测其三级结构模型,β亚基采用折叠识别法预测其三级结构模型,预测结果质量评估均较好。此外,分析了不同亚基的表面电荷分布,并预测了最可能的异柠檬酸和ADP的结合位点。通过生物信息学方法进行合理预测,为进一步研究IDH家族尤其是哺乳动物NAD(+)IDH提供了重要理论依据。     

赤链蛇Sox基因的克隆及序列分析

参照人SRY基因HMG-box保守区的序列,设计一对兼并引物,采用PCR技术扩增了赤链蛇的Sox基因,并对扩增产物进行了克隆和测序,结果在雌雄个体中共筛选出三个Sox基因,其中有一个为雌雄共有,显示出性别差异性;三个Sox基因编码的氨基酸序列与人相应SOX3,SOX4,SOX22基因的相似性分别为96%,98%,96%。显示出Sox基因在进化上的高度保守性,本文为赤链蛇的性别决定机制研究提供了分子资料。     

基于生物信息学方法预测人线粒体转录终止因子3蛋白的结构与功能

目的:基于生物信息学预测人线粒体转录终止因子3(hMTERF3)蛋白的结构与功能。方法:利用GenBank、Uniprot、ExPASy、SWISS-PROT数据库资源和不同的生物信息学软件对hMTERF3蛋白进行系统研究,包括hMTERF3的理化性质、跨膜区和信号肽、二级结构功能域、亚细胞定位、蛋白质的功能分类预测、同源蛋白质多重序列比对、系统发育树构建、三级结构同源建模。结果:软件预测hMTERF3蛋白的相对分子质量为47.97×103,等电点为8.60,不具信号肽和跨膜区;二级结构分析显示主要为螺旋和无规则卷曲,包含6个MTERF基序,三级结构预测结果与二级结构预测结果相符;亚细胞定位分析结果显示该蛋白定位于人线粒体;功能分类预测其为转运和结合蛋白,参与基因转录调控;同源蛋白质多重序列比对和进化分析显示,hMTERF3蛋白与大鼠、小鼠等哺乳动物的MTERF3蛋白具有高度同源性,在系统发育树上聚为一类。结论:hMTERF3蛋白的生物信息学分析为进一步开展对该蛋白的结构和功能的实验研究提供了理论依据。     

人STC2蛋白的生物学分析

运用生物信息学方法预测分析不同物种间STC2蛋白的理化性质、同源性以及人STC2蛋白亲水性、核定位序列,跨膜区域、信号肽结构、亚细胞定位,二级结构、三级结构、互作蛋白、GO注释分析。STC2蛋白由302个氨基酸组成,理论等电点为6.93,具有较强亲水性,在哺乳动物中较为保守,不存在核定位序列和跨膜结构,含有信号肽,主要集中在细胞内质网或分泌到胞外;STC2蛋白二级结构预测有11个α螺旋区和1个β折叠区,拉曼图表明三级预测结构可靠,互作蛋白及GO注释提示STC2可参与多种细胞生物学过程。通过对人STC2蛋白结构和功能的预测分析,为STC2蛋白的进一步研究提供一定的理论依据,也为STC2相关疾病的诊治提供新的思路。     

2019新型冠状病毒S蛋白的结构和功能分析

运用生物信息学,预测急性呼吸系统综合症冠状病毒2(SARS-CoV-2/2019-nCoV)的基本理化性质、结构、功能和抗原表位等,为新型冠状病毒肺炎(COVID-19)的防治提供思路。应用ExPASy分析S蛋白的消光系数、不稳定系数和半衰期等理化性质;利用SignaIP v5.0分析S蛋白的信号肽;应用TMHMM分析S蛋白的跨膜区;利用NetPhos3.1在线工具预测S蛋白的磷酸化位点;应用Pfam预测S蛋白的结构域;应用PSIPRED分析S蛋白的二级结构特征;利用SWISS-MODEL构建S蛋白的三级结构;利用BLAST分析SARS-CoV-2的S蛋白与其他物种的相似性;利用MEGA软件分析2019-nCoV的S蛋白与其他物种的进化关系。S蛋白由1 273个氨基酸组成,其相对分子质量为141 178.47,等电点为6.24,含有一个跨膜区,是低亲水性分泌蛋白;S蛋白的基本组成单位为纤突蛋白,其二级结构中以无规则卷曲和螺旋结构为主,三级结构中纤突糖蛋白和ACE2复合体具有重要的意义;2019-nCoV与蝙蝠冠状病毒和SARS-CoV同源;S蛋白存在多个潜在的线性T细胞和B细胞表位,1 202～1 210位氨基酸区域的抗原性和应答频率最高。生物信息学技术有利于了解S蛋白的理化性质、结构、功能和潜在的线性T细胞表位等,可为新型冠状肺炎的研究和防治提供参考依据。     

家蝇小热休克蛋白(sHsp20.6)的生物信息学分析

研究家蝇小热休克蛋白sHsp20.6的生物学功能。方法应用生物信息学的方法和工具对家蝇sHsp20.6的理化性质、疏水性、跨膜区和信号肽、膜体分析、二级结构功能域、蛋白质的功能分类预测、多重序列比对与系统发育树构建、三级结构建模进行分析。结果表明:家蝇sHsp20.6是一个亲水蛋白,分子量为20.64kD,等电点为5.66,不具有跨膜区和信号肽,包含有一个HSP20的结构域,主要构成原件为α螺旋和无规则卷曲,三维结构预测显示该蛋白为棒状结构,C端结构域具有7个片层结构。聚类分析显示,家蝇sHsp20.6蛋白与昆虫中的直系同源小热休克蛋白(orthologoussmallheatshockprotein)聚为一类。     

肺癌中SOX4基因的调控网络分析验证和蛋白质预测

目的:分析SOX4基因的调控网络和蛋白产物特征,并进行验证。方法:利用生物信息学方法预测人SOX4基因的调控网络及其蛋白产物与其他蛋白的相互作用;在肺癌细胞系和肺癌临床样本中通过实时荧光定量PCR对mi RNA-30s进行定量分析。结果:人类SOX4蛋白可能与TP53、DICER1和SDCBP2等3种蛋白相互作用;SOX4基因可能受多种lnc RNA和mi RNA的调控;hsa-mi R30s在临床样本中表达模式不完全相同。与癌旁组织相比,hsami R30c-5p和hsa-mi R30d-5p在癌组织中显示低表达,hsa-mi R30d-5p则相反;与正常人胚肺成纤维细胞KMB17相比,肺癌细胞系A549和H1299中3个mi RNA-30成员的表达模式均为高表达。结论:SOX4基因及其蛋白产物的生物信息学分析为相关研究提供了重要信息基础,但尚须相关实验验证,以确定相关性。     

Tritium planigraphy: Differences in the spatial structures of the influenza virus M1 protein in crystal,solution, and virion

The structure of the M1 protein of the influenza virus A/Puerto Rico/8/34 (PR8, subtype H1N1) in solution at acidic pH and in the composition of the virion has been studied by the tritium planigraphy method. A model of the spatial structure was constructed using a special algorithm simulating the experiment and a set of algorithms for predicting the secondary structure and disordered regions in proteins. The tertiary structure was refined using the Rosetta program. For a comparison of the structures in solution and inside the virion, the data of X-ray diffraction analysis for the NM domain were also used. The main difference in the structures of the protein in solution and the crystalline state is observed in the region of contact of N and M domains, which in the crystalline state is packed more densely. The regions of the maximum label incorporation almost completely coincide with unstructured regions in the protein that were predicted by the bioinformatics analysis. These regions are concentrated in the C domain and in loop regions between M, N, and C domains. The data were confirmed by analytical centrifugation and dynamic light scattering. Anomalous hydrodynamic dimensions and a low structuration of the M1 protein in solution were found. The polyfunctionality of the protein in the cell is probably related to its flexible tertiary structure, which, owing to unstructured regions, provides contact with various partner molecules.     

Differences in spatial structures of the influenza virus M1 protein in crystal, solution and virion

Spatial structure of the influenza virus A/Puerto Rico/8/34 (PR8, subtype H1N1) M1 protein in a solution and composition of the virion was studied by tritium planigraphy technique. The special algorithm for modeling of the spatial structure is used to simulate the experiment, as well as a set of algorithms predicting secondary structure and disordered regions in proteins. Tertiary structures were refined using the program Rosetta. To compare the structures in solution and in virion, also used the X-ray diffraction data for NM-domain. The main difference between protein structure in solution and crystal is observed in the contact region of N- and M-domains, which are more densely packed in the crystalline state. Locations include the maximum label is almost identical to the unstructured regions of proteins predicted by bioinformatics analysis. These areas are concentrated in the C-domain and in the loop regions between the M-, N-, and C-domains. Analytical centrifugation and dynamic laser light scattering confirm data of tritium planigraphy. Anomalous hydrodynamic size, and low structuring of the M1 protein in solution were found. The multifunctionality of protein in the cell appears to be associated with its plastic tertiary structure, which provides at the expense of unstructured regions of contact with various molecules-partners.     

A novel approach to segregate and identify functional loop regions in protein structures using their Ramachandran maps

The loops which connect or flank helices/sheets in protein structures are known to be functionally important. However, ironically they also belong to the part of protein whose structure is least accurately predicted. Here, a new method to isolate and analyze loop regions in protein structure is proposed using the spatial coordinates of the solved three‐dimensional structure. The extent of dispersion among points of successive amino acid residues in the Ramachandran map of protein region is utilized to calculate the Mean Separation between these points in the Ramachandran Plot (MSRP). Based on analysis of 2935 protein secondary structure regions obtained using DSSP software, spanning a range from 2 to 64 residues, taken from a set of 170 proteins, it is shown that helices (MSRP < 17) and strands (MSRP < 64) stand effectively demarcated from the loop regions (MSRP > 130). Analysis of 43 DNA binding and 98 ligand binding proteins revealed several loop regions with clear change in MSRP subsequent to binding. The population of such loops correlated with the magnitude of backbone displacement in the protein subsequent to binding. Can changes in MSRP quantify the temporal oscillations in dihedral angles among structured/unstructured regions in proteins? Molecular dynamics simulations (10 ns) revealed that deviations in MSRP among different snapshots in the trajectory were at least twofold higher for unstructured proteins in comparison with ordered proteins. The above results validate the use of MSRP parameter as a tool to identify and investigate functionally active loops and unstructured regions in protein structures. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

东北虎中五个Sox基因的保守区序列分析(英文)

在哺乳动物中,Y染色体决定着雄性性别,这是由于在其短臂上存在一个编码睾丸决定因子 (TDF) 的基因。1990年,人们克隆了睾丸决定因子基因并命名为SRY。序列分析表明SRY基因中存在一个保守区,与染色体高迁移率组 (HMG) 蛋白质上的DNA结合结构域具有一定的相似性。基于HMG基序的保守性人们发现了一个新的基因家族Sox基因家族。凡是在HMG区域与SRY基因具有50%以上相似性的基因被称为Sox基因。Sox基因在早期胚胎发育过程中参与多种发育途径,具有重要的作用。参与诸如性别决定、骨组织的发育、血细胞生成过程、神经系统的发育、晶状体的发育等多种早期胚胎发育过程。 虎(Panthera tigris)作为世界上最濒危的兽类之一,东北虎(Panthera tigris altaica Temminck)是其中的一个亚种,被列为一类珍稀动物。本文对其发育基因家族—SOX基因进行了研究。 利用肌肉组织为材料制备基因组DNA, 应用特异于HMG-box区域的兼并引物, 扩增了东北虎的SOX基因。在扩增产物中发现一条220bp的扩增带。经过克隆与序列测定和同源性检索,发现5个基因片段(Fig.1&2)。其与人SRY基因的相似性分别为75%、56%、51%、67% 和48%;与小鼠Sry基因的相似性分别为73%、54%、57%、66% 和 48% (Table 1)。因此这5个DNA片段为东北虎的5个Sox基因片     

A data bank merging related protein structures and sequences.

A data collection which merges protein structural and sequence information is described. Structural superpositions amongst proteins with similar main-chain fold were performed or collected from the literature. Sequences taken from the protein primary structure databases were associated with the multiple structural alignments providing they were at least 50% homologous in residue identity to one of the structural sequences and at least 50% of the structural sequence residues were alignable. Such restrictions allow reasonable confidence that the primary sequences share the conformation of the tertiary structural templates, except in the less conserved loop regions. Multiple structural superpositions were collected for 38 familial groups containing a total of 209 tertiary structures; 45 structures had no superposable mates and were used individually. Other information is also provided as main-chain and side-chain conformational angles, secondary structural assignments and the like. Wedding the primary and tertiary structural data resulted in an 8-fold increase of data bank sequence entries over those associated with the known three-dimensional architectures alone.     

Identifying the Tertiary Fold of Small Proteins with Different Topologies from Sequence and Secondary Structure using the Genetic Algorithm and Extended Criteria Specific for Strand Regions

Grid-free protein folding simulations based on sequence and secondary structure knowledge (using mostly experimentally determined secondary structure information but also analysing results from secondary structure predictions) were investigated using the genetic algorithm, a backbone representation, and standard dihedral angular conformations. Optimal structures are selected according to basic protein building principles. Having previously applied this approach to proteins with helical topology, we have now developed additional criteria and weights for β-strand- containing proteins, validated them on four small β-strand-rich proteins with different topologies, and tested the general performance of the method on many further examples from known protein structures with mixed secondary structural type and less than 100 amino acid residues.Topology predictions close to the observed experimental structures were obtained in four test cases together with fitness values that correlated with the similarity of the predicted topology to the observed structures. Root-mean-square deviation values of C^αatoms in the superposed predicted and observed structures, the latter of which had different topologies, were between 4.5 and 5.5 Å (2.9 to 5.1 Å without loops). Including 15 further protein examples with unique folds, root-mean-square deviation values ranged between 1.8 and 6.9 Å with loop regions and averaged 5.3 Å and 4.3 Å, including and excluding loop regions, respectively.     

Molecular modeling and characterization of the B. thuringiensis and B. thuringiensis LDC-9 cytolytic proteins

The Cyt toxins are able to lyse a wide range of cell types in vitro, unlike the Cry delta-endotoxins. It exerts its activity by the formation of pores within target cell membranes. The structural information available for Cyt2Aa (PDB id: 1CBY) consists of a single domain in which two outer layers of alpha-helix wrap around a mixed beta-sheet. Beta-barrel was suggested as a possible structure of the pores. Hence, this study seeks to investigate the structural properties of other Cytolytic proteins by predicting the three-dimensional (3D) model using Cyt2Aa as template. The predicted models are expected to be significantly more accurate as all the Cyt proteins showed significant similarity with the template (PDB id: 1CBY). The refined homology models revealed similar secondary structures (alpha-helices and beta-sheets) and tertiary features as Cyt2Aa. The variation in the loop regions of the tertiary structure accounts for the differential toxicity.     

Knowledge based modelling of homologous proteins, Part I: Three-dimensional frameworks derived from the simultaneous superposition of multiple structures

An approach is described for modelling the three-dimensional structure of a protein from the tertiary structures of several homologous proteins that have been determined by X-ray analysis. A method is developed for the simultaneous superposition of several protein molecules and for the calculation of an 'average structure' or 'framework'. Investigation of the convergence properties of this method, in the case of both weighted and unweighted least squares, demonstrates that both give a unique answer and the latter is robust for an homologous family of proteins. Multi-dimensional scaling is used to subgroup of the proteins with respect to structural homology. The framework calculated on the basis of the family of homologous proteins, or of an appropriate subgroup, is used to align fragments of the known protein structures of high sequence homology with the unknown. This alignment provides a basis for model building the tertiary structure. Different techniques for using the framework to model the mainchain of various globins and an immunoglobulin domain in the structurally conserved regions are investigated.     

Three-dimensional model of Escherichia coli ribosomal 5 S RNA as deduced from structure probing in solution and computer modeling.

The conformation of Escherichia coli 5 S rRNA was investigated using chemical and enzymatic probes. The four bases were monitored at one of their Watson-Crick positions with dimethylsulfate (at C(N-3) and A(N-1], with a carbodiimide derivative (at G(N-1) and U(N-3] and with kethoxal (at G(N-1, N-2]. Position N-7 of purine was probed with diethylpyrocarbonate (at A(N-7] and dimethylsulfate (at G(N-7]. Double-stranded or stacked regions were tested with RNase V1 and unpaired guanine residues with RNase T1. We also used lead(II) that has a preferential affinity for interhelical and loop regions and a high sensitivity for flexible regions. Particular care was taken to use uniform conditions of salt, magnesium, pH and temperature for the different enzymatic chemical probes. Derived from these experimental data, a three dimensional model of the 5 S rRNA was built using computer modeling which integrates stereochemical constraints and phylogenetic data. The three domains of 5 S rRNA secondary structure fold into a Y-shaped structure that does not accommodate long-range tertiary interactions between domains. The three domains have distinct structural and dynamic features as revealed by the chemical reactivity and the lead(II)-induced hydrolysis: domain 2 (loop B/helix III/loop C) displays a rather weak structure and possesses dynamic properties while domain 3 (helix V/region E/helix IV/loop D) adopts a highly structured and overall helical conformation. Conserved nucleotides are not crucial for the tertiary folding but maintain an intrinsic structure in the loop regions, especially via non-canonical pairing (A.G, G.U, G.G, A.C, C.C), which can close the loops in a highly specific fashion. In particular, nucleotides in the large external loop C fold into an organized conformation leading to the formation of a five-membered loop motif. Finally, nucleotides at the hinge region of the Y-shape are involved in a precise array of hydrogen bonds based on a triple interaction between U14, G69 and G107 stabilizing the quasi-colinearity of helices II and V. The proposed tertiary model is consistent with the localization of the ribosomal protein binding sites and possesses strong analogy with the model proposed for Xenopus laevis 5 S rRNA, indicating that the Y-shape model can be generalized to all 5 S rRNAs.