Fold assignments for amino acid sequences of the CASP2 experiment

Interaction of the red pigment-concentrating hormone of the crustacean Daphnia pulex, with its cognate receptor, Dappu-RPCHR: A nuclear magnetic resonance and modeling study

The primary sequence of the red pigment-concentrating hormone (RPCH) receptor of the water flea, Daphnia pulex, was used in homology modeling to construct the first 3D model of a crustacean G-protein coupled receptor, Dappu-RPCHR. This receptor was found to belong to the class A subfamily of GPCRs with a disulfide bridge between Cys⁷² and Cys¹⁵⁰ and an ionic lock between Arg⁹⁷ and Thr²²⁴ and Thr²²⁰. NMR restrained molecular dynamics was used to determine the structure of an agonist, Dappu-RPCH, in a membrane-mimicking environment. The agonist was found to be flexible but has two main conformations in solution, both having β-turns. Docking of the predominant structure was used to find a binding pocket on the receptor. The pocket's spatial location was similar to that of the AKH receptor of Anopheles gambiae. The binding affinity was -69kcalmol^-1 with the N-terminus of Dappu-RPCH inserted between helices 4 and 6, and the C-terminus interacting with extra-cellular loop, ECL2. Upon binding, H-bonding to the peptide may activate the receptor. This development of the first Dappu-RPCH/Dappu-RPCHR model could be useful for understanding ligand-receptor interactions in crustaceans.

Structural genomics: computational methods for structure analysis

The success of structural genomics initiatives requires the development and application of tools for structure analysis, prediction, and annotation. In this paper we review recent developments in these areas; specifically structure alignment, the detection of remote homologs and analogs, homology modeling and the use of structures to predict function. We also discuss various rationales for structural genomics initiatives. These include the structure-based clustering of sequence space and genome-wide function assignment. It is also argued that structural genomics can be integrated into more traditional biological research if specific biological questions are included in target selection strategies.

Structure of bovine alpha-1,3-galactosyltransferase and its complexes with UDP and DPGal inferred from molecular modeling

A homology model of alpha-1,3-galactosyltransferase (alpha-1,3-GalT), the retaining enzyme responsible for the formation of alpha-galactosyl epitopes, has been developed by means of molecular modeling using the SpsA glycosyltransferase structure. A protein-ligand docking approach was used to model alpha-1,3-GalT complexed with UDP and UDP-Gal. The comparison of structural features found in the alpha-1,3-GalT homology model with available structural data on this class of enzymes revealed similarities in the UDP-binding pocket. In the predicted structure of the complexes, the pyrophosphate interacts with the DVD motif (Asp-225, Val-226, and Asp-227) of alpha-1,3-GalT through the Mn(2+) cation. The uridine part of the UDP binds into the well-defined cavity that consists of Phe-134, Tyr-139, Ile-140, Val-136, Arg-194, Arg-202, Lys-209, Asp-173, His-218, and Thr-137 in a conformation that is generally observed in the crystal structures of other glycosyltransferase complexes.

PRED-CLASS: cascading neural networks for generalized protein classification and genome-wide applications

A cascading system of hierarchical, artificial neural networks (named PRED-CLASS) is presented for the generalized classification of proteins into four distinct classes-transmembrane, fibrous, globular, and mixed-from information solely encoded in their amino acid sequences. The architecture of the individual component networks is kept very simple, reducing the number of free parameters (network synaptic weights) for faster training, improved generalization, and the avoidance of data overfitting. Capturing information from as few as 50 protein sequences spread among the four target classes (6 transmembrane, 10 fibrous, 13 globular, and 17 mixed), PRED-CLASS was able to obtain 371 correct predictions out of a set of 387 proteins (success rate approximately 96%) unambiguously assigned into one of the target classes. The application of PRED-CLASS to several test sets and complete proteomes of several organisms demonstrates that such a method could serve as a valuable tool in the annotation of genomic open reading frames with no functional assignment or as a preliminary step in fold recognition and ab initio structure prediction methods. Detailed results obtained for various data sets and completed genomes, along with a web sever running the PRED-CLASS algorithm, can be accessed over the World Wide Web at http://o2.biol.uoa.gr/PRED-CLASS.

Activity of human IgG and IgA subclasses in immune defense against Neisseria meningitidis serogroup B

Both IgG and IgA Abs have been implicated in host defense against bacterial infections, although their relative contributions remain unclear. We generated a unique panel of human chimeric Abs of all human IgG and IgA subclasses with identical V genes against porin A, a major subcapsular protein Ag of Neisseria meningitidis and a vaccine candidate. Chimeric Abs were produced in baby hamster kidney cells, and IgA-producing clones were cotransfected with human J chain and/or human secretory component. Although IgG (isotypes IgG1-3) mediated efficient complement-dependent lysis, IgA was unable to. However, IgA proved equally active to IgG in stimulating polymorphonuclear leukocyte respiratory burst. Remarkably, although porin-specific monomeric, dimeric, and polymeric IgA triggered efficient phagocytosis, secretory IgA did not. These studies reveal unique and nonoverlapping roles for IgG and IgA Abs in defense against meningococcal infections.

Use of residue pairs in protein sequence-sequence and sequence-structure alignments

Two new sets of scoring matrices are introduced: H2 for the protein sequence comparison and T2 for the protein sequence-structure correlation. Each element of H2 or T2 measures the frequency with which a pair of amino acid types in one protein, k-residues apart in the sequence, is aligned with another pair of residues, of given amino acid types (for H2) or in given structural states (for T2), in other structurally homologous proteins. There are four types, corresponding to the k-values of 1 to 4, for both H2 and T2. These matrices were set up using a large number of structurally homologous protein pairs, with little sequence homology between the pair, that were recently generated using the structure comparison program SHEBA. The two scoring matrices were incorporated into the main body of the sequence alignment program SSEARCH in the FASTA package and tested in a fold recognition setting in which a set of 107 test sequences were aligned to each of a panel of 3,539 domains that represent all known protein structures. Six procedures were tested; the straight Smith-Waterman (SW) and FASTA procedures, which used the Blosum62 single residue type substitution matrix; BLAST and PSI-BLAST procedures, which also used the Blosum62 matrix; PASH, which used Blosum62 and H2 matrices; and PASSC, which used Blosum62, H2, and T2 matrices. All procedures gave similar results when the probe and target sequences had greater than 30% sequence identity. However, when the sequence identity was below 30%, a similar structure could be found for more sequences using PASSC than using any other procedure. PASH and PSI-BLAST gave the next best results.

A structural model for the rolA protein and its interaction with DNA

The study of the plant oncogene rolA has been hampered by a lack of structural information. Here we show that, despite a lack of significant sequence similarity to proteins of known structure, the rolA sequence adopts a known fold; that of the papillomavirus E2 DNA-binding domain. This fold is reliably identified by modern threading programs, which consider predicted secondary structure, but not by others. Although the rolA sequence is only around 16% identical to those of the available template structures, a structural model could be built that performed well against protein structure verification programs. The adopted strategy involved alignment corrections, justified by multiple model building and evaluation, with particular attention paid to the hydrophobic core residues. We find that rolA protein is predicted to resemble the template proteins in two key aspects; existence as a dimer and ability to bind DNA. rolA protein has recently been shown experimentally to possess DNA binding ability. This model predicts Lys 24 and Arg 27 to be involved in sequence-specific interactions and eight other residues to hydrogen-bond phosphate groups of the DNA.

Factors limiting the performance of prediction-based fold recognition methods

In the past few years, a new generation of fold recognition methods has been developed, in which the classical sequence information is combined with information obtained from secondary structure and, sometimes, accessibility predictions. The results are promising, indicating that this approach may compete with potential-based methods (Rost B et al., 1997, J Mol Biol 270:471-480). Here we present a systematic study of the different factors contributing to the performance of these methods, in particular when applied to the problem of fold recognition of remote homologues. Our results indicate that secondary structure and accessibility prediction methods have reached an accuracy level where they are not the major factor limiting the accuracy of fold recognition. The pattern degeneracy problem is confirmed as the major source of error of these methods. On the basis of these results, we study three different options to overcome these limitations: normalization schemes, mapping of the coil state into the different zones of the Ramachandran plot, and post-threading graphical analysis.

Structure and function prediction of the Brucella abortus P39 protein by comparative modeling with marginal sequence similarities

A methodology is proposed to solve a difficult modeling problem related to the recently sequenced P39 protein. This sequence shares no similarity with any known 3D structure, but a fold is proposed by several threading tools. The difficulty in aligning the target sequence on one of the proposed template structures is overcome by combining the results of several available prediction methods and by refining a rational consensus between them. In silico validation of the obtained model and a preliminary cross-check with experimental features allow us to state that this borderline prediction is at least reasonable. This model raises relevant hypotheses on the main structural features of the protein and allows the design of site-directed mutations. Knowing the genetic context of the P39 reading frame, we are now able to suggest a function for the P39 protein: it would act as a periplasmic substrate-binding protein.