Search PubMedSearch

Biomedical subjects

S H Bryant

Publications and source records attributed to S H Bryant.

At least 19 recordsLinked to original sources

MMDB: Entrez's 3D structure database.

The three dimensional structures for representatives of nearly half of all protein families are now available in public databases. Thus, no matter which protein one investigates, it is increasingly likely that the 3D structure of a homolog will be known and may reveal unsuspected structure-function relationships. The goal of Entrez's 3D-structure database is to make this information accessible and usable by molecular biologists (http://www.ncbi.nlm.nih.gov/Entrez). To this end Entrez provides two major analysis tools, a search engine based on sequence and structure 'neighboring' and an integrated visualization system for sequence and structure alignments. From a protein's sequence 'neighbors' one may rapidly identify other members of a protein family, including those where 3D structure is known. By comparing aligned sequences and/or structures in detail, using the visualization system, one may identify conserved features and perhaps infer functional properties. Here we describe how these analysis tools may be used to investigate the structure and function of newly discovered proteins, using the PTEN gene product as an example.

Amino Acid Sequence

Structure databases.

Explore the source record for details and available documents.

Computational Biology

Measures of threading specificity and accuracy.

Threading predictions for CASP2 target proteins were compared to their true structures using a series of precisely defined measures of agreement, calculated in a fully automatic way. Fold recognition specificity was calculated as the proportion of a predictor's "bet" that was placed on previously-known structures similar to the prediction target, as identified by a "jury" of well-tested structure-structure comparison methods. Values approaching 100% indicate that a prediction correctly identified the structural and/or evolutionary family to which a target belongs. Alignment specificity was calculated as the proportion of aligned residue paris in the predicted target-to-known-structure alignment that also occur in the structure-structure alignments produced by the "jury" methods. Contact specificity was calculated as the proportion of nonlocal residue contacts in the molecular model implied by threading alignment, that also occur in the experimental structure of the target. Alignment specificity and contact specificity measure the accuracy of a predicted 3-dimensional model. Values approaching 100% indicate that target residues have been assigned to the correct spatial locations and that the model is as accurate as possible for a threading prediction.

Amino Acid Sequence

A retrospective analysis of CASP2 threading predictions.

Analysis of CASP2 protein threading results shows that the success rate of structure predictions varies widely among prediction targets. We set "critical" thresholds in fold recognition specificity and threading model accuracy at the points where "incorrect" CASP2 predictions just outnumber "correct" predictions. Using these thresholds we find that correct predictions were made for all of those targets and for only those targets where more than 50% of target residues may be superimposed on previously known structures. Three-fourths of these correct predictions were furthermore made for targets with greater than 12% residue identity in structural alignment, where characteristic sequence motifs are also present. Based on these observations we suggest that the sustained performance of threading methods is best characterized by counting the numbers of correct predictions for targets of increasing "difficulty." We suggest that target difficulty may be assigned, once the true structure of the target is known, according to the fraction of residues superimposable onto previously known structures and the fraction of identical residues in those structural alignments.

Models, Molecular

Extent and nature of contacts between protein molecules in crystal lattices and between subunits of protein oligomers.

A survey was compiled of several characteristics of the intersubunit contacts in 58 oligomeric proteins, and of the intermolecular contracts in the lattice for 223 protein crystal structures. The total number of atoms in contact and the secondary structure elements involved are similar in the two types of interfaces. Crystal contact patches are frequently smaller than patches involved in oligomer interfaces. Crystal contacts result from more numerous interactions by polar residues, compared with a tendency toward nonpolar amino acids at oligomer interfaces. Arginine is the only amino acid prominent in both types of interfaces. Potentials of mean force for residue-residue contacts at both crystal and oligomer interfaces were derived from comparison of the number of observed residue-residue interactions with the number expected by mass action. They show that hydrophobic interactions at oligomer interfaces favor aromatic amino acids and methionine over aliphatic amino acids; and that crystal contacts form in such a way as to avoid inclusion of hydrophobic interactions. They also suggest that complex salt bridges with certain amino acid compositions might be important in oligomer formation. For a protein that is recalcitrant to crystallization, substitution of lysine residues with arginine or glutamine is a recommended strategy.

Crystallization

A measure of success in fold recognition.

Prediction of protein structure by fold recognition, or threading, was recently put to the test in a 'blind' structure prediction experiment, CASP2. Thirty-two teams from around the world participated, preparing predictions for 22 different 'target' proteins whose structures were soon to be determined. As experimental structures became available, we, as organizers of the threading competition, computed objective measures of fold-recognition specificity and model accuracy, to identify and characterize successful predictions. Here, we present a brief summary of these prediction evaluations, a tally of 'correct' predictions and a discussion of factors associated with correct predictions. We find that threading produced specific recognition and accurate models whenever the structural database contained a template spanning a large fraction of target sequence. Presence of conserved sequence motifs was helpful, but not required, and it would appear that threading can succeed whenever similarity to a known structure is sufficiently extensive.

Computer Simulation

Evaluation of threading specificity and accuracy.

Threading experiments with proteins from the globin family provide an indication of the nature of the structural similarity required for successful fold recognition and accurate sequence-structure alignment. Threading scores are found to rise above the noise of false positives whenever roughly 60% of residues from a sequence can be aligned with analogous sites in the structure of a remote homolog. Fold recognition specificity thus appears to be limited by the extent of structural similarity, regardless of the degree of sequence similarity. Threading alignment accuracy is found to depend more critically on the degree of structural similarity. Alignments are accurate, placing the majority of residues exactly as in structural alignment, only when superposition residuals are less than 2.5 A. These criteria for successful recognition and sequence-structure alignment appear to be consistent with the successes and failures of threading methods in blind structure prediction. They also suggest a direct assay for improved threading methods: Potentials and alignment models should be tested for their ability to detect less extensive structural similarities, and to produce accurate alignments when superposition residuals for this conserved "core" fall in the range characteristic of remote homologs.

Evaluation Studies as Topic

Surprising similarities in structure comparison.

Examination of a protein's structural 'neighbors' can reveal distant evolutionary relationships that are otherwise undetectable, and perhaps suggest unsuspected functional properties. In the past, such analyses have often required specialized software and computer skills, but new structural comparison methods, developed in the past two years, increasingly offer this opportunity to structural and molecular biologists in general. These methods are based on similarity-search algorithms that are fast enough to have effectively removed the computer-time limitation for structure-structure search and alignment, and have made it possible for several groups to conduct systematic comparisons of all publicly available structures, and offer this information via the World Wide Web. Furthermore, and perhaps surprisingly given the difficulty of the structure-comparison problem, these groups seem to have converged on quite similar approaches with respect to both fast search algorithms and the identification of statistically significant similarities.

Computer Communication Networks

Threading analysis suggests that the obese gene product may be a helical cytokine.

The ob gene encodes a protein that, in mutant form, is associated with obesity and type II diabetes in mice. Sequence analysis has revealed no similarities to other proteins, however, and no clues as to possible functions. The possibility nonetheless remains that ob is functionally or ancestrally related to other proteins, whose sequences are divergent to the point that only a comparison of three-dimensional structures might detect relationship. To explore this possibility, we conduct a 'threading' search of a 3-dimensional structure database, to determine whether the ob protein might adopt a fold similar to any known structure. This search reveals that the ob sequence is compatible, at a significance level of P < 0.05, with structures from the family of helical cytokines that includes interleukin-2 and growth hormone. A structural model of ob based upon these results is physically and biologically plausible and leads to testable predictions, including the prediction that ob may activate the JAK-STAT pathway, via binding to a receptor resembling those of the cytokine family.

Amino Acid Sequence

Homology model building of the HMG-1 box structural domain.

Nucleoproteins belonging to the HMG-1/2 family possess homologous domains approximately 75 amino acids in length. These domains, termed HMG-1 boxes, are highly structured, compact, and mediate the interaction between HMG-1 box-containing proteins and DNA in a variety of biological contexts. Homology model building experiments on HMG-1 box sequences 'threaded' through the 1H-NMR structure of an HMG-1 box from rat indicate that the domain does not have rigid sequence requirements for its formation. Energy calculations indicate that the structure of all HMG-1 box domains is stabilized primarily through hydrophobic interactions. We have found structural relationships in the absence of statistically significant sequence similarity, identifying several candidate proteins which could possibly assume the same three-dimensional conformation as the rat HMG-1 box motif. The threading technique provides a method by which significant structural similarities in a diverse protein family can be efficiently detected, and the 'structural alignment' derived by this method provides a rational basis through which phylogenetic relationships and the precise sites of interaction between HMG-1 box proteins and DNA can be deduced.

Amino Acid Sequence

A proposed structural model of domain 1 of fasciclin III neural cell adhesion protein based on an inverse folding algorithm.

Fasciclin III is an integral membrane protein expressed on a subset of axons in the developing Drosophila nervous system. It consists of an intracellular domain, a transmembrane region, and an extracellular region composed of three domains, each predicted to form an immunoglobulin-like fold. The most N-terminal of these domains is expected to be important in mediating cell-cell recognition events during nervous system development. To learn more about the structure/function relationships in this cellular recognition molecule, a model structure of this domain was built. A sequence-to-structure alignment algorithm was used to align the protein sequence of the fasciclin III first domain to the immunoglobulin McPC603 structure. Based on this alignment, a model of the domain was built using standard homology modeling techniques. Side-chain conformations were automatically modeled using a rotamer search algorithm and the model was minimized to relax atomic overlaps. The resulting model is compact and has chemical characteristics consistent with related globular protein structures. This model is a de novo test of the sequence-to-structure alignment algorithm and is currently being used as the basis for mutagenesis experiments to discern the parts of the fasciclin III protein that are necessary for homophilic molecular recognition in the developing Drosophila nervous system.

Algorithms

Threading a database of protein cores.

We present an analysis of 10 blind predictions prepared for a recent conference, "Critical Assessment of Techniques for Protein Structure Prediction." The sequences of these proteins are not detectably similar to those of any protein in the structure database then available, but we attempted, by a threading method, to recognize similarity to known domain folds. Four of the 10 proteins, as we subsequently learned, do indeed show significant similarity to then-known structures. For 2 of these proteins the predictions were accurate, in the sense that a similar structure was at or near the top of the list of threading scores, and the threading alignment agreed well with the corresponding structural alignment. For the best predicted model mean alignment error relative to the optimal structural alignment was 2.7 residues, arising entirely from small "register shifts" of strands or helices. In the analysis we attempt to identify factors responsible for these successes and failures. Since our threading method does not use gap penalties, we may readily distinguish between errors arising from our prior definition of the "cores" of known structures and errors arising from inherent limitations in the threading potential. It would appear from the results that successful substructure recognition depends most critically on accurate definition of the "fold" of a database protein. This definition must correctly delineate substructures that are, and are not, likely to be conserved during protein evolution.

Amino Acid Sequence

Statistics of sequence-structure threading.

The past two years have seen the rapid development of new recognition methods for protein structure prediction. These algorithms 'thread' the sequence of one protein through the known structure of another, looking for an alignment that corresponds to an energetically favorable model structure. Because they are based on energy calculation, rather than evolutionary distance, these methods extend the possibility of structure prediction by comparative modeling to a larger class of new sequences, where similarity to known structures is recognizable by no other means. The strength of the evidence they offer should be judged by objective statistical tests, however, so as to rule out the possibility that favorable scores arise from chance factors such as similarity of length, composition, or the consideration of a large number of alternative alignments. Calculation of objective p-values by analytical means is not yet possible, but it would appear that approximate values may be obtained by simulation, as they are in gapped, global sequence alignment. We propose that the results of threading experiments should include Z-scores relative to the composition-corrected score distribution obtained for shuffled and optimally aligned sequences.

Algorithms

Eukaryotic translation elongation factor 1 gamma contains a glutathione transferase domain--study of a diverse, ancient protein superfamily using motif search and structural modeling.

Using computer methods for multiple alignment, sequence motif search, and tertiary structure modeling, we show that eukaryotic translation elongation factor 1 gamma (EF1 gamma) contains an N-terminal domain related to class theta glutathione S-transferases (GST). GST-like proteins related to class theta comprise a large group including, in addition to typical GSTs and EF1 gamma, stress-induced proteins from bacteria and plants, bacterial reductive dehalogenases and beta-etherases, and several uncharacterized proteins. These proteins share 2 conserved sequence motifs with GSTs of other classes (alpha, mu, and pi). Tertiary structure modeling showed that in spite of the relatively low sequence similarity, the GST-related domain of EF1 gamma is likely to form a fold very similar to that in the known structures of class alpha, mu, and pi GSTs. One of the conserved motifs is implicated in glutathione binding, whereas the other motif probably is involved in maintaining the proper conformation of the GST domain. We predict that the GST-like domain in EF1 gamma is enzymatically active and that to exhibit GST activity, EF1 gamma has to form homodimers. The GST activity may be involved in the regulation of the assembly of multisubunit complexes containing EF1 and aminoacyl-tRNA synthetases by shifting the balance between glutathione, disulfide glutathione, thiol groups of cysteines, and protein disulfide bonds. The GST domain is a widespread, conserved enzymatic module that may be covalently or noncovalently complexed with other proteins. Regulation of protein assembly and folding may be 1 of the functions of GST.

Amino Acid Sequence

Effects of veratridine on Na and Ca currents in frog skeletal muscle.

1. Voltage-clamp experiments were performed to determine the effects of veratridine on Na and Ca currents in frog skeletal muscle fibres. 2. Veratridine (1 microM) did not affect the kinetics of the fast Na current but it did induce a slowly inactivating tetrodotoxin-sensitive inward current that was apparent after Na current inactivation. This slow current had a peak amplitude of 6.7 +/- 0.7 microA/cm2 at -20 mV and decayed monoexponentially with a time constant of 606 +/- 77 ms. 3. The slow current had a voltage-dependence for activation that was similar to that of the fast Na current. Single depolarizing prepulses that induced complete inactivation of the fast Na channels, prevented development of the slow current. Trains of brief depolarizations at increasing frequencies increased the amplitude of the slow current. These results suggest that the slow current may be mediated by veratridine modified Na channels that must be in the open position. 4. The low concentration of veratridine (1 microM) did not affect the Ca current, while 100 microM veratridine reversibly suppressed the Ca current and shifted its peak current-voltage relation towards more negative potentials. Thus, veratridine appears not to be a selective fast Na channel modifier as it may also alter Ca channel gating properties in skeletal muscle fibres.

Animals

Effects of veratrine on ion currents in single rabbit cardiomyocytes.

1. Voltage-clamp experiments were performed to determine the effects of veratrine (1 microgram/ml) on Na and K currents in isolated rabbit ventricular cardiomyocytes. 2. Veratrine did not affect the inward rectifier K current, increased the inactivation time constant of the transient outward current (I(to)) and induced a slowly decaying inward current component (Iv), which was sensitive to tetrodotoxin. 3. Inactivation of fast Na channels by application of short depolarizing prepulses to potentials between -90 and -50 mV prevented the development of Iv.Iv decayed biexponentially with time constants equal to 139 +/- 9.0 ms and 776 +/- 47 ms. The net amplitude of Iv and the time constants for its rapidly and slowly inactivating components were little affected by trains of conditioning prepulses to 0 mV. The contributions, however, of the fast and slow components to the net current were significantly altered by repetitive depolarizations. 4. These components of Iv are likely due to modification of open cardiac Na channels by veratrum alkaloids.

Animals