Relationship between multiple sequence alignments and quality of protein comparative models

Comparative modeling is the method of choice, whenever applicable, for protein structure prediction, not only because of its higher accuracy compared to alternative methods, but also because it is possible to estimate a priori the quality of the models that it can produce, thereby allowing the usefulness of a model for a given application to be assessed beforehand. By and large, the quality of a comparative model depends on two factors: the extent of structural divergence between the target and the template and the quality of the sequence alignment between the two protein sequences. The latter is usually derived from a multiple sequence alignment (MSA) of as many proteins of the family as possible, and its accuracy depends on the number and similarity distribution of the sequences of the protein family. Here we describe a method to evaluate the expected difficulty, and by extension accuracy, of a comparative model on the basis of the MSA used to build it. The parameter that we derive is used to compare the results obtained in the last two editions of the Critical Assessment of Methods for Structure Prediction (CASP) experiment as a function of the difficulty of the modeling exercise. Our analysis demonstrates that the improvement in the scope and quality of comparative models between the two experiments is largely due to the increased number of available protein sequences and to the consequent increased chance that a large and appropriately spaced set of protein sequences homologous to the proteins of interest is available.     

BAliBASE 3.0: latest developments of the multiple sequence alignment benchmark

Multiple sequence alignment is one of the cornerstones of modern molecular biology. It is used to identify conserved motifs, to determine protein domains, in 2D/3D structure prediction by homology and in evolutionary studies. Recently, high-throughput technologies such as genome sequencing and structural proteomics have lead to an explosion in the amount of sequence and structure information available. In response, several new multiple alignment methods have been developed that improve both the efficiency and the quality of protein alignments. Consequently, the benchmarks used to evaluate and compare these methods must also evolve. We present here the latest release of the most widely used multiple alignment benchmark, BAliBASE, which provides high quality, manually refined, reference alignments based on 3D structural superpositions. Version 3.0 of BAliBASE includes new, more challenging test cases, representing the real problems encountered when aligning large sets of complex sequences. Using a novel, semiautomatic update protocol, the number of protein families in the benchmark has been increased and representative test cases are now available that cover most of the protein fold space. The total number of proteins in BAliBASE has also been significantly increased from 1444 to 6255 sequences. In addition, full-length sequences are now provided for all test cases, which represent difficult cases for both global and local alignment programs. Finally, the BAliBASE Web site (http://www-bio3d-igbmc.u-strasbg.fr/balibase) has been completely redesigned to provide a more user-friendly, interactive interface for the visualization of the BAliBASE reference alignments and the associated annotations.     

Tissue arrays as fiducial markers for section alignment in 3-D reconstruction technology

Combination of conventional histology and the three-dimensional spatial view of tissue structures offers new prospects for understanding and diagnosing nature and development of human diseases. The essential technical problem related to three-dimensional reconstruction in histopathology is represented by the correct alignment of serial sections. During the past years several methods have been proposed but failed to become popular because of their limits in terms of time consume and restricted applicability. We aimed to overcome this problem by applying the technology of Tissue Array, thus by positioning adequate fiducial markers from specific "donor" blocks into the "recipient" paraffin block of interest. Digitized pictures of serially cut sections were aligned according to the tissue markers embedded by Tissue Array, and then processed with specific softwares for three-dimensional reconstruction. Thirteen models, including fetal hearts, breast and thyroid carcinomas, were elaborated. We found the procedure to be easy, fast and reproducible. Moreover, by selectively embedding the fiducial markers according to specific angles, the Tissue Arrays can be exploited in order to establish the distance between sections. This original methodology of incorporating Tissue Arrays into paraffin blocks as fiducial markers for three-dimensional reconstruction has a potential impact on histology for research purposes and diagnostic applications.     

Analysis of FMRFamide-like peptide (FLP) diversity in phylum Nematoda

This study reports a series of systematic BLAST searches of nematode ESTs on the Genbank database, using search strings derived from known nematode FLPs (those encoded by Caenorhabditis elegans flp genes as well as those isolated from other nematodes including Ascaris suum), as well as query sequences representative of theoretical FLPs. Over 1000 putative FLP-encoding ESTs were identified from multiple nematode species. A total of 969 ESTs representing sequelogs of the 23 known C. elegans flp genes were identified in 32 species, from clades I, III, IV and V. Numerical analysis of EST numbers suggests that flp-1, flp-11 and flp-14 are amongst the most highly expressed flp genes. Speculative BLAST searches were performed using theoretical FLP C-termini as queries, in an attempt to identify putative novel FLP sequences in the EST database. These searches yielded eight multi-species sequelogs encoding FLPs with novel signatures that are believed to identify distinct flp genes. These novel genes encode 25 distinct previously unidentified FLPs, and raise the current total of known nematode flp genes to 31. Additionally, software-based analyses of the presence of signal peptides were performed, with signal peptides being identified on at least one member of each group of flp ESTs, further confirming their status as secreted peptides. The data reveal that nematode FLPs encompass the most complex neuropeptide family known within the metazoa. Moreover, individual FLPs and FLP motifs are highly conserved across the nematodes with little evidence for inter-clade or inter-lifestyle variation, supporting their fundamental role in free-living and parasitic species.     

Phylogenetic analysis of family 6 glycoside hydrolases

Multiple sequence alignment separates members of glycoside hydrolase Family 6 into eight subfamilies: one of mainly actinobacterial endoglucanases (EGs), one of ascomycotal EGs, one of chytridiomycotal EGs and cellobiohydrolases (CBHs), one of actinobacterial and proteobacterial CBHs, one of chytridiomycotal CBHs, two of ascomycotal CBHs, and one of basidiomycotal CBHs. Each also has some proteins of unknown function. Multiple sequence alignment also extends to all of Family 6 the observation that lengths of loops that form the active-site tunnel in CBHs vary among subfamilies, and along with loop conformations, determine enzyme function.     

Biochemical problems of regulation by oligopeptides

Some problems are considered which arise in biochemical studies on structure and function of natural oligopeptides consisting of 2–50 amino acid residues. The problems under consideration include the generation of oligopeptides from precursors, chemical structure, the role of functionally important radicals and spatial configuration, and structure-function relationships. Different types of regulation are shown mainly for oligopeptides involved in muscle contraction.Translated from Biokhimiya, Vol. 69, No. 11, 2004, pp. 1565–1573.Original Russian Text Copyright © 2004 by Zamyatnin.     

Testing for marginal independence between two categorical variables with multiple responses

Questions that ask respondents to "choose all that apply" from a set of items occur frequently in surveys. Categorical variables that summarize this type of survey data are called both pick any/c variables and multiple-response categorical variables. It is often of interest to test for independence between two categorical variables. When both categorical variables can have multiple responses, traditional Pearson chi-square tests for independence should not be used because of the within-subject dependence among responses. An intuitively constructed version of the Pearson statistic is proposed to perform the test using bootstrap procedures to approximate its sampling distribution. First- and second-order adjustments to the proposed statistic are given in order to use a chi-square distribution approximation. A Bonferroni adjustment is proposed to perform the test when the joint set of responses for individual subjects is unavailable. Simulations show that the bootstrap procedures hold the correct size more consistently than the other procedures.     

Soret Spectral and Bioinformatic Approaches Provide Evidence for a Critical Role of the α -Subunit in Assembly of Tetrameric Hemoglobin

Soret spectral contributions of the α-subunit heme pocket have been evaluated by performing static titrations of apohemoglobin A with CNProtohemin under varied experimental conditions. Increasing the temperature from 5 to 30°C in 0.05 M potassium phosphate buffer, pH 7, resulted in a decreasingly prominent hypsochromic shifts reflecting altered the vinyl–globin interactions. Studies at 10°C in over pH range of 6.7–8.0 revealed a profile for the spectral shifts approximating the side chain pK value (7.4) a histidyl residue. These overall spectral changes correspond to ΔE of ≤7 kJ/mol indicative of electrostatic noncovalent interactions. Further our current molecular modeling studies indicate that the spatial arrangement and critical noncovalent interactions of tyrosine 42 and histidine 45 (aromatic residues unique to the α-subunit) make significant contribution to the Soret spectra. Most interestingly, phylogenetic analyses have revealed the presence of a histidyl triad in the α-chain of all vertebrates that form heterotetramers.     

A data-mining approach for multiple structural alignment of proteins

Comparing the 3D structures of proteins is an important but computationally hard problem in bioinformatics. In this paper, we propose studying the problem when much less information or assumptions are available. We model the structural alignment of proteins as a combinatorial problem. In the problem, each protein is simply a set of points in the 3D space, without sequence order information, and the objective is to discover all large enough alignments for any subset of the input. We propose a data-mining approach for this problem. We first perform geometric hashing of the structures such that points with similar locations in the 3D space are hashed into the same bin in the hash table. The novelty is that we consider each bin as a coincidence group and mine for frequent patterns, which is a well-studied technique in data mining. We observe that these frequent patterns are already potentially large alignments. Then a simple heuristic is used to extend the alignments if possible. We implemented the algorithm and tested it using real protein structures. The results were compared with existing tools. They showed that the algorithm is capable of finding conserved substructures that do not preserve sequence order, especially those existing in protein interfaces. The algorithm can also identify conserved substructures of functionally similar structures within a mixture with dissimilar ones. The running time of the program was smaller or comparable to that of the existing tools.     

Model-Guided Mutagenesis Drives Functional Studies of Human NHA2, Implicated in Hypertension

Human NHA2 is a poorly characterized Na⁺/H⁺ antiporter recently implicated in essential hypertension. We used a range of computational tools and evolutionary conservation analysis to build and validate a three-dimensional model of NHA2 based on the crystal structure of a distantly related bacterial transporter, NhaA. The model guided mutagenic evaluation of transport function, ion selectivity, and pH dependence of NHA2 by phenotype screening in yeast. We describe a cluster of essential, highly conserved titratable residues located in an assembly region made of two discontinuous helices of inverted topology, each interrupted by an extended chain. Whereas in NhaA, oppositely charged residues compensate for partial dipoles generated within this assembly, in NHA2, polar but uncharged residues suffice. Our findings led to a model for transport mechanism that was compared to the well-known electroneutral NHE1 and electrogenic NhaA subtypes. This study establishes NHA2 as a prototype for the poorly understood, yet ubiquitous, CPA2 antiporter family recently recognized in plants and metazoans and illustrates a structure-driven approach to derive functional information on a newly discovered transporter.