Universal similarity measure for comparing protein structures

We introduce a new variant of the root mean square distance (RMSD) for comparing protein structures whose range of values is independent of protein size. This new dimensionless measure (relative RMSD, or RRMSD) is zero between identical structures and one between structures that are as globally dissimilar as an average pair of random polypeptides of respective sizes. The RRMSD probability distribution between random polypeptides converges to a universal curve as the chain length increases. The correlation coefficients between aligned random structures are computed as a function of polypeptide size showing two characteristic lengths of 4.7 and 37 residues. These lengths mark the separation between phases of different structural order between native protein fragments. The implications for threading are discussed.     

A new topological method to measure protein structure similarity

A method for the quantitative evaluation of structural similarity between protein pairs is developed that makes use of a Delaunay-based topological mapping. The result of the mapping is a three-dimensional array which is representative of the global structural topology and whose elements can be used to construe an integral scoring scheme. This scoring scheme was tested for its dependence on the protein length difference in a pairwise comparison, its ability to provide a reasonable means for structural similarity comparison within a family of structural neighbors of similar length, and its sensitivity to the differences in protein conformation. It is shown that such a topological evaluation of similarity is capable of providing insight into these points of interest. Protein structure comparison using the method is computationally efficient and the topological scores, although providing different information about protein similarity, correlate well with the distance root-mean-square deviation values calculated by rigid-body structural alignment.     

Recognizing misfolded and distorted protein structures by the assumption-based similarity score

A new similarity score (sigma-score) is proposed which is able to find the correct protein structure among the very close alternatives and to distinguish between correct and deliberately misfolded structures. This score is based on the general principle 'similar likes similar', and it favors hydrophobic and hydrophilic contacts, and disfavors hydrophobic-to-hydrophilic contacts in proteins. The values of sigma-scores calculated for the high-resolution protein structures from the representative set are compared with those of alternatives: (i) very close alternatives which are only slightly distorted by conformational energy minimization in vacuo; (ii) alternatives with subsequently growing distortions, generated by molecular dynamics simulations in vacuo; (iii) structures derived by molecular dynamics simulation in solvent at 300 K; (iv) deliberately misfolded protein models. In nearly all tested cases the similarity score can successfully distinguish between experimental structure and its alternatives, even if the root mean square displacement of all heavy atoms is less than 1 A. The confidence interval of the similarity score was estimated using the high-resolution X-ray structures of domain pairs related by non-crystallographic symmetry. The similarity score can be used for the evaluation of the general quality of the protein models, choosing the correct structures among the very close alternatives, characterization of models simulating folding/unfolding, etc.     

PSI: indexing protein structures for fast similarity search

MOTIVATION: We consider the problem of finding similarities in protein structure databases. Current techniques sequentially compare the given query protein to all of the proteins in the database to find similarities. Therefore, the cost of similarity queries increases linearly as the volume of the protein databases increase. As the sizes of experimentally determined and theoretically estimated protein structure databases grow, there is a need for scalable searching techniques. RESULTS: Our techniques extract feature vectors on triplets of SSEs (Secondary Structure Elements). Later, these feature vectors are indexed using a multidimensional index structure. For a given query protein, this index structure is used to quickly prune away unpromising proteins in the database. The remaining proteins are then aligned using a popular alignment tool such as VAST. We also develop a novel statistical model to estimate the goodness of a match using the SSEs. Experimental results show that our techniques improve the pruning time of VAST 3 to 3.5 times while maintaining similar sensitivity.     

Evolution and similarity evaluation of protein structures in contact map space

Prediction of fold from amino acid sequence of a protein has been an active area of research in the past few years, but the limited accuracy of existing techniques emphasizes the need to develop newer approaches to tackle this task. In this study, we use contact map prediction as an intermediate step in fold prediction from sequence. Contact map is a reduced graph-theoretic representation of proteins that models the local and global inter-residue contacts in the structure. We start with a population of random contact maps for the protein sequence and "evolve" the population to a "high-feasibility" configuration using a genetic algorithm. A neural network is employed to assess the feasibility of contact maps based on their 4 physically relevant properties. We also introduce 5 parameters, based on algebraic graph theory and physical considerations, that can be used to judge the structural similarity between proteins through contact maps. To predict the fold of a given amino acid sequence, we predict a contact map that will sufficiently approximate the structure of the corresponding protein. Then we assess the similarity of this contact map with the representative contact map of each fold; the fold that corresponds to the closest match is our predicted fold for the input sequence. We have found that our feasibility measure is able to differentiate between feasible and infeasible contact maps. Further, this novel approach is able to predict the folds from sequences significantly better than a random predictor.     

Comparison of gene structures and enzymatic properties between two endoglucanases from Humicola grisea

We have cloned two endoglucanase genes (egl3 and egl4) from a thermophilic fungus, Humicola grisea. The coding region of the egl3 gene was interrupted by an intron of 56-bp, and the deduced amino acid sequence of the egl3 gene was 305 amino acids in length and showed 98.4% identity with Humicola insolens EGV. The coding region of the egl4 gene was also interrupted by an intron of 173-bp, which contains 34 TTC repeated sequence units, and the deduced amino acid sequence of the egl4 gene was 227 amino acids in length and showed 61.5% identity with H. grisea EGL3. The typical hinge and the cellulose-binding domain were observed in the C-terminal region of EGL3, but they were not observed in EGL4. In the 5′ upstream region of both genes, there were a TATA box or its similar sequence, CAAT motifs, and 6-bp sites which are identical or similar to the consensus sequence for binding a catabolite repressor CREA in Aspergillus nidulans. The egl3 and the egl4 genes were expressed in Aspergillus oryzae, and the translation products were purified. The fusion protein, EGL4CBD, which consists of a catalytic domain of EGL4 and the C-terminal region of EGL3, was also constructed and produced by A. oryzae, and purified. These enzymes showed relatively high activity toward carboxymethyl cellulose (CMC) and could not hydrolyze p-nitrophenyl-β-d-glucoside and p-nitrophenyl-β-d-cellobioside. The positive effect of substituting the C-terminal region of EGL4 with that of EGL3 was observed in the hydrolysis of CMC.     

Molecular docking of thiamine reveals similarity in binding properties between the prion protein and other thiamine-binding proteins

Prion-induced diseases are a global health concern. The lack of effective therapy and 100 % mortality rates for such diseases have made the prion protein an important target for drug discovery. Previous NMR experimental work revealed that thiamine and its derivatives bind the prion protein in a pocket near the N-terminal loop of helix 1, and conserved intermolecular interactions were noted between thiamine and other thiamine-binding proteins. Furthermore, water-mediated interactions were observed in all of the X-ray crystallographic structures of thiamine-binding proteins, but were not observed in the thiamine–prion NMR study. To better understand the potential role of water in thiamine–prion binding, a docking study was employed using structural X-ray solvent. Before energy minimization, docked thiamine assumed a “V” shape similar to some of the known thiamine-dependent proteins. Following minimization with NMR-derived restraints, the “F” conformation was observed. Our findings confirmed that water is involved in ligand stabilization and phosphate group interaction. The resulting refined structure of thiamine bound to the prion protein allowed the 4-aminopyrimidine ring of thiamine to π-stack with Tyr150, and facilitated hydrogen bonding between Asp147 and the amino group of 4-aminopyrimidine. Investigation of the π-stacking interaction through mutation of the tyrosine residue further revealed its importance in ligand placement. The resulting refined structure is in good agreement with previous experimental restraints, and is consistent with the pharmacophore model of thiamine-binding proteins.     

A new definition of life

Chirality is often glossed over in theoretical or experimental discussions concerning the origin of life, but the ubiquity of homochiral building blocks in known biological systems demands explanation. Information theory can provide a quantitative framework for understanding the role of chirality in biology. Here I show how conclusions derived from information theory, in particular the concept of equivocation, can explain not only why chiral building blocks are necessary in living systems but also why a homochiral set of building blocks is necessary. These results lead to a new definition of life, and to the conclusion that the simplest form of life exists in the form of self-amplifying, autocatalytic reactions such as the Soai reaction.     

A new definition of genetic distance

Summary A stochastic measure of genetic distance between populations is proposed; unlike currently used measures, it is invariant with respect to union and subdivision of loci. This measure enables a unified quantitative approach to genetic diversity within populations and genetic distance between populations.     

Sequence similarity between protein B and human apolipoprotein A-IV

Sequence comparison of protein B (CAMP-factor) with human apolipoprotein A-IV (apo A-IV) revealed 32% similarity between the N-terminal part of protein B and a part of the putative lipid-binding domain of apo A-IV. The significance of this similarity is discussed with respect to the structure/function relationship of protein B.     

Syntheses, structures and magnetic properties of two new metal complexes based on a pyridyl-diphosphonate ligand

Two new first-raw transition metal diphosphonate complexes, namely, {[Ni₃([hpyedpH)₂(H₂O)₄]·(H₂O)₂}_n (1) and [Mn[hpyedpH₂](H₂O)]_n (2), based on a multidentate ligand 1-hydroxy-2-(3-pyridyl)-ethylidene-1,1-diphosphonic acid (hpyedpH₄) have been synthesized under hydrothermal reaction and structurally characterized by single-crystal X-ray diffraction, IR spectroscopy and element analyses. The data reveals that complex 1 is a 2D layer structure, whereas the complex 2 possesses a 1D motif. The powder X-ray diffraction (PXRD) patterns for complexes 1 and 2 were collected as well, which match well with the ones calculated from their single-crystal structure data. Magnetic measurements show that complex 1 is a ferrimagnet with T_c = 5.0 K. Magnetic studies of complex 2 indicate antiferromagnetic behavior.     

GASH: An improved algorithm for maximizing the number of equivalent residues between two protein structures

Background  

We introduce GASH, a new, publicly accessible program for structural alignment and superposition. Alignments are scored by the Number of Equivalent Residues (NER), a quantitative measure of structural similarity that can be applied to any structural alignment method. Multiple alignments are optimized by conjugate gradient maximization of the NER score within the genetic algorithm framework. Initial alignments are generated by the program Local ASH, and can be supplemented by alignments from any other program.     

A comprehensive analysis of non-sequential alignments between all protein structures

Background  

The majority of relations between proteins can be represented as a conventional sequential alignment. Nevertheless, unusual non-sequential alignments with different connectivity of the aligned fragments in compared proteins have been reported by many researchers. It is interesting to understand those non-sequential alignments; are they unique, sporadic cases or they occur frequently; do they belong to a few specific folds or spread among many different folds, as a common feature of protein structure. We present here a comprehensive large-scale study of non-sequential alignments between available protein structures in Protein Data Bank.