利用序列保守模体和局部构象信息预测转录因子结合位点

转录因子结合位点的计算预测是研究基因转录调控的重要环节,但常用的位置特异得分矩阵方法预测特异性偏低.通过深入分析结合位点的生物特征,提出了一种综合利用序列保守模体和局部构象信息的结合位点预测方法,以极大相关得分矩阵作为保守模体的描述模型,并根据二苷参数模型计算位点序列的局部构象,将两类信息得分组合为多维特征向量,在二次判别分析的框架下进行训练和滑动预测.预测过程中还引入了位置信息量以优化似然得分和过滤备选结果.针对大肠杆菌CRP和Fis结合位点数据的留一法测试结果表明,描述模型的改进和多种信息的融合能有效地改善预测方法的性能,大幅度提高特异性.     

Entropy capacity determines protein folding

Search and study of the general principles that govern kinetics and thermodynamics of protein folding generate a new insight into the factors controlling this process. Here, based on the known experimental data and using theoretical modeling of protein folding, we demonstrate that there exists an optimal relationship between the average conformational entropy and the average energy of contacts per residue-that is, an entropy capacity-for fast protein folding. Statistical analysis of conformational entropy and number of contacts per residue for 5829 protein structures from four general structural classes (all-alpha, all-beta, alpha/beta, alpha+beta) demonstrates that each class of proteins has its own class-specific average number of contacts (class alpha/beta has the largest number of contacts) and average conformational entropy per residue (class all-alpha has the largest number of rotatable angles phi, psi, and chi per residue). These class-specific features determine the folding rates: alpha proteins are the fastest folding proteins, then follow beta and alpha+beta proteins, and finally alpha/beta proteins are the slowest ones. Our result is in agreement with the experimental folding rates for 60 proteins. This suggests that structural and sequence properties are important determinants of protein folding rates.     

NMR investigation of DNA primer-template models: structural insights into dislocation mutagenesis in DNA replication

Slipped frameshift intermediates can occur when DNA polymerase slows or stalls at sites of DNA lesions. However, this phenomenon is much less common when unmodified DNA is replicated. In order to study the effect of templating bases on the alignment of primer-templates, NMR structural investigation has been performed on primer-template oligonucleotide models which mimic the situation that dNTP has just been incorporated opposite template. NMR evidence reveals the occurrence of misalignment when dGTP is incorporated opposite template T with a downstream nucleotide C. Depending on the template sequence, further extension of the primer can lead to realignment.     

STITCHER: Dynamic assembly of likely amyloid and prion β‐structures from secondary structure predictions

The supersecondary structure of amyloids and prions, proteins of intense clinical and biological interest, are difficult to determine by standard experimental or computational means. In addition, significant conformational heterogeneity is known or suspected to exist in many amyloid fibrils. Previous work has demonstrated that probability‐based prediction of discrete β‐strand pairs can offer insight into these structures. Here, we devise a system of energetic rules that can be used to dynamically assemble these discrete β‐strand pairs into complete amyloid β‐structures. The STITCHER algorithm progressively ‘stitches’ strand‐pairs into full β‐sheets based on a novel free‐energy model, incorporating experimentally observed amino‐acid side‐chain stacking contributions, entropic estimates, and steric restrictions for amyloidal parallel β‐sheet construction. A dynamic program computes the top 50 structures and returns both the highest scoring structure and a consensus structure taken by polling this list for common discrete elements. Putative structural heterogeneity can be inferred from sequence regions that compose poorly. Predictions show agreement with experimental models of Alzheimer's amyloid beta peptide and the Podospora anserina Het‐s prion. Predictions of the HET‐s homolog HET‐S also reflect experimental observations of poor amyloid formation. We put forward predicted structures for the yeast prion Sup35, suggesting N‐terminal structural stability enabled by tyrosine ladders, and C‐terminal heterogeneity. Predictions for the Rnq1 prion and alpha‐synuclein are also given, identifying a similar mix of homogenous and heterogeneous secondary structure elements. STITCHER provides novel insight into the energetic basis of amyloid structure, provides accurate structure predictions, and can help guide future experimental studies. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization

Water plays an essential role in determining the structure and function of all biological systems. Recent methodological advances allow for an accurate and efficient estimation of the thermodynamic properties of water molecules at the surface of proteins. In this work, we characterize these thermodynamic properties and relate them to various structural and functional characteristics of the protein. We find that high-energy hydration sites often exist near protein motifs typically characterized as hydrophilic, such as backbone amide groups. We also find that waters around alpha helices and beta sheets tend to be less stable than waters around loops. Furthermore, we find no significant correlation between the hydration site-free energy and the solvent accessible surface area of the site. In addition, we find that the distribution of high-energy hydration sites on the protein surface can be used to identify the location of binding sites and that binding sites of druggable targets tend to have a greater density of thermodynamically unstable hydration sites. Using this information, we characterize the FKBP12 protein and show good agreement between fragment screening hit rates from NMR spectroscopy and hydration site energetics. Finally, we show that water molecules observed in crystal structures are less stable on average than bulk water as a consequence of the high degree of spatial localization, thereby resulting in a significant loss in entropy. These findings should help to better understand the characteristics of waters at the surface of proteins and are expected to lead to insights that can guide structure-based drug design efforts.     

The dynamical response of hen egg white lysozyme to the binding of a carbohydrate ligand

It has become clear that the binding of small and large ligands to proteins can invoke significant changes in side chain and main chain motion in the fast picosecond to nanosecond timescale. Recently, the use of a "dynamical proxy" has indicated that changes in these motions often reflect significant changes in conformational entropy. These entropic contributions are sometimes of the same order as the total entropy of binding. Thus, it is important to understand the connections amongst motion between the manifold of states accessible to the native state of proteins, the corresponding entropy, and how this impacts the overall energetics of protein function. The interaction of proteins with carbohydrate ligands is central to a range of biological functions. Here, we examine a classic carbohydrate interaction with an enzyme: the binding of wild-type hen egg white lysozyme (HEWL) to the natural, competitive inhibitor chitotriose. Using NMR relaxation experiments, backbone amide and side chain methyl axial order parameters were obtained across apo and chitotriose-bound HEWL. Upon binding, changes in the apparent amplitude of picosecond to nanosecond main chain and side chain motions are seen across the protein. Indeed, binding of chitotriose renders a large contiguous fraction of HEWL effectively completely rigid. Changes in methyl flexibility are most pronounced closest to the binding site, but average to only a small overall change in the dynamics across the protein. The corresponding change in conformational entropy is unfavorable and estimated to be a significant fraction of the total binding entropy.     

一种新的蛋白质结构类预测方法

基于氨基酸的16种分类模型,给出蛋白质序列的派生序列,进而结合加权拟熵和LZ复杂度构造出34维特征向量来表示蛋白质序列。借助于贝叶斯分类器对同源性不超过25%的640数据集进行蛋白质结构类预测,准确度达到71.28%。     

Structural insights into the recognition of peroxisomal targeting signal 1 by Trypanosoma brucei peroxin 5

Glycosomes are peroxisome-like organelles essential for trypanosomatid parasites. Glycosome biogenesis is mediated by proteins called “peroxins,” which are considered to be promising drug targets in pathogenic Trypanosomatidae. The first step during protein translocation across the glycosomal membrane of peroxisomal targeting signal 1 (PTS1)-harboring proteins is signal recognition by the cytosolic receptor peroxin 5 (PEX5). The C-terminal PTS1 motifs interact with the PTS1 binding domain (P1BD) of PEX5, which is made up of seven tetratricopeptide repeats. Obtaining diffraction-quality crystals of the P1BD of Trypanosoma brucei PEX5 (TbPEX5) required surface entropy reduction mutagenesis. Each of the seven tetratricopeptide repeats appears to have a residue in the α_L conformation in the loop connecting helices A and B. Five crystal structures of the P1BD of TbPEX5 were determined, each in complex with a hepta- or decapeptide corresponding to a natural or nonnatural PTS1 sequence. The PTS1 peptides are bound between the two subdomains of the P1BD. These structures indicate precise recognition of the C-terminal Leu of the PTS1 motif and important interactions between the PTS1 peptide main chain and up to five invariant Asn side chains of PEX5. The TbPEX5 structures reported here reveal a unique hydrophobic pocket in the subdomain interface that might be explored to obtain compounds that prevent relative motions of the subdomains and interfere selectively with PTS1 motif binding or release in trypanosomatids, and would therefore disrupt glycosome biogenesis and prevent parasite growth.     

The application of the entropy-based statistic for genomic association study of QTL

An entropy-based statistic TPE has been proposed for genomic association study for disease-susceptibility locus.The statistic TPE may be directly adopted and/or extended to quantitative-trait locus (QTL)mapping for quantitative traits.In this article,the statistic TPE was extended and applied to quantitative trait for association analysis of QTL by means of selective genotyping.The statistical properties (the type I error rate and the power) were examined under a range of parameters and population-sampling strategies (e.g.,various genetic models,various heritabilities,and various sample-selection threshold values) by simulation studies.The results indicated that the statistic Tee is robust and powerful for genomic association study of QTL.A simulation study based on the haplotype frequencies of 10 single nucleotide polymorphisms (SNPs) of angiotensin-I converting enzyme genes was conducted to evaluate the performance of the statistic TPE for genetic association study.     

An empirical EEG analysis in brain death diagnosis for adults

Electroencephalogram (EEG) is often used in the confirmatory test for brain death diagnosis in clinical practice. Because EEG recording and monitoring is relatively safe for the patients in deep coma, it is believed to be valuable for either reducing the risk of brain death diagnosis (while comparing other tests such as the apnea) or preventing mistaken diagnosis. The objective of this paper is to study several statistical methods for quantitative EEG analysis in order to help bedside or ambulatory monitoring or diagnosis. We apply signal processing and quantitative statistical analysis for the EEG recordings of 32 adult patients. For EEG signal processing, independent component analysis (ICA) was applied to separate the independent source components, followed by Fourier and time-frequency analysis. For quantitative EEG analysis, we apply several statistical complexity measures to the EEG signals and evaluate the differences between two groups of patients: the subjects in deep coma, and the subjects who were categorized as brain death. We report statistically significant differences of quantitative statistics with real-life EEG recordings in such a clinical study, and we also present interpretation and discussions on the preliminary experimental results.