Hb M Dothan [beta 25/26 (B7/B8)/(GGT/GAG-->GAG//Gly/Glu-->Glu]; a new mechanism of unstable methemoglobin variant and molecular characteristics.

A new unstable beta globin chain variant associated with methemoglobin (Met-Hb) phenotype was found in a Caucasian infant. Molecular analysis of the beta globin gene using polymerase chain reaction (PCR) amplification and sequencing led to the detection of a new in frame deletion in exon-1. Direct sequencing of the PCR product revealed a 3 bp deletion (-GTG) between codons 25/26, which resulted in the loss of a single amino acid (-Gly). We propose that this newly discovered unstable M-hemoglobin (M-Hb) variant, named Hb Dothan [GGT/GAG-->GAG//Gly/Glu-->Glu], is caused by a shift in the amino acid sequence and altered packing of the B and E helices during beta globin synthesis, and also changes the orientation of the critical proximal and distal histidine in the F and E helices respectively. Phenotype/Genotype features and molecular characteristics of this new beta chain are presented in this communication.

Structure prediction. How good are we?

Recent successes show that, in certain circumstances, protein secondary structures can be predicted with high accuracy. How far are we from being able to predict the complete structure of a protein from its sequence?

Functional analysis of H-Ryk, an atypical member of the receptor tyrosine kinase family.

H-Ryk is an atypical receptor tyrosine kinase which differs from other members of this family at a number of conserved residues in the activation and nucleotide binding domains. Using a chimeric receptor approach, we demonstrate that H-Ryk has impaired catalytic activity. Despite the receptor's inability to undergo autophosphorylation or phosphorylate substrates, we demonstrate that ligand stimulation of the chimeric receptor results in activation of the mitogen-activated protein kinase pathway. The ability to transduce signals is abolished by mutation of the invariant lysine (K334A) in subdomain II of H-Ryk. Further, by in vitro mutagenesis, we show that the amino acid substitutions in the activation domain of H-Ryk account for the loss of catalytic activity. In addition to the essential aspartate residue, either phenylalanine or glycine is required in the activation domain to maintain proper conformation of the catalytic domain and thus ensure receptor autophosphorylation. Homology modelling of the catalytic domain of H-Ryk provides a rationale for these findings. Thus, the signalling properties of H-Ryk are divergent from those of other classical receptor tyrosine kinases.

Protein fold irregularities that hinder sequence analysis.

The detection of homologous protein sequences frequently provides useful predictions of function and structure. Methods for homology searching have continued to improve, such that very distant evolutionary relationships can now be detected. Little attention has been paid, however, to the problems of detecting homology when domains are inserted or permuted. Here we review recent occurrences of these phenomena and discuss methods that permit their detection.

GABA(B2) is essential for g-protein coupling of the GABA(B) receptor heterodimer.

GABA(B) receptors are unique among G-protein-coupled receptors (GPCRs) in their requirement for heterodimerization between two homologous subunits, GABA(B1) and GABA(B2), for functional expression. Whereas GABA(B1) is capable of binding receptor agonists and antagonists, the role of each GABA(B) subunit in receptor signaling is unknown. Here we identified amino acid residues within the second intracellular domain of GABA(B2) that are critical for the coupling of GABA(B) receptor heterodimers to their downstream effector systems. Our results provide strong evidence for a functional role of the GABA(B2) subunit in G-protein coupling of the GABA(B) receptor heterodimer. In addition, they provide evidence for a novel "sequential" GPCR signaling mechanism in which ligand binding to one heterodimer subunit can induce signal transduction through the second partner of a heteromeric complex.

The C-terminal domains of the GABA(b) receptor subunits mediate intracellular trafficking but are not required for receptor signaling.

GABA(B) receptors are G-protein-coupled receptors that mediate slow synaptic inhibition in the brain and spinal cord. These receptors are heterodimers assembled from GABA(B1) and GABA(B2) subunits, neither of which is capable of producing functional GABA(B) receptors on homomeric expression. GABA(B1,) although able to bind GABA, is retained within the endoplasmic reticulum (ER) when expressed alone. In contrast, GABA(B2) is able to access the cell surface when expressed alone but does not couple efficiently to the appropriate effector systems or produce any detectable GABA-binding sites. In the present study, we have constructed chimeric and truncated GABA(B1) and GABA(B2) subunits to explore further GABA(B) receptor signaling and assembly. Removal of the entire C-terminal intracellular domain of GABA(B1) results in plasma membrane expression without the production of a functional GABA(B) receptor. However, coexpression of this truncated GABA(B1) subunit with either GABA(B2) or a truncated GABA(B2) subunit in which the C terminal has also been removed is capable of functional signaling via G-proteins. In contrast, transferring the entire C-terminal tail of GABA(B1) to GABA(B2) leads to the ER retention of the GABA(B2) subunit when expressed alone. These results indicate that the C terminal of GABA(B1) mediates the ER retention of this protein and that neither of the C-terminal tails of GABA(B1) or GABA(B2) is an absolute requirement for functional coupling of heteromeric receptors. Furthermore although GABA(B1) is capable of producing GABA-binding sites, GABA(B2) is of central importance in the functional coupling of heteromeric GABA(B) receptors to G-proteins and the subsequent activation of effector systems.

Detection of protein three-dimensional side-chain patterns: new examples of convergent evolution.

Detection of recurring three-dimensional side-chain patterns is a potential means of inferring protein function. This paper presents a new method for detecting such patterns and discusses various implications. The method allows detection of side-chain patterns without any prior knowledge of function, requiring only protein structure data and associated multiple sequence alignments. A recursive, depth-first search algorithm finds all possible groups of identical amino acids common to two protein structures independent of sequence order. The search is highly constrained by distance constraints, and by ignoring amino acids unlikely to be involved in protein function. A weighted root-mean-square deviation (RMSD) between equivalenced groups of amino acids is used as a measure of similarity. The statistical significance of any RMSD is assigned by reference to a distribution fitted to simulated data. Searches with the Ser/His/Asp catalytic triad, a His/His porphyrin binding pattern, and the zinc-finger Cys/Cys/His/His pattern are performed to test the method on known examples. An all-against-all comparison of representatives from the structural classification of proteins (SCOP) is performed, revealing several new examples of evolutionary convergence to common patterns of side-chains within different tertiary folds and in different orders along the sequence. These include a di-zinc binding Asp/Asp/His/His/Ser pattern common to alkaline phosphatase/bacterial aminopeptidase, and an Asp/Glu/His/His/Asn/Asn pattern common to the active sites of DNase I and endocellulase E1. Implications for protein evolution, function prediction and the rational design of functional regulators are discussed.

The limits of protein secondary structure prediction accuracy from multiple sequence alignment.

The expected best residue-by-residue accuracies for secondary structure prediction from multiple protein sequence alignment have been determined by an analysis of known protein structural families. The results show substantial variation is possible among homologous protein structures, and that 100% agreement is unlikely between a consensus prediction and one member of a protein structural family. The study provides the range of agreement to be expected between a perfect secondary structure prediction from a multiple alignment and each protein within the alignment. The results of this study overcome the difficulties inherent in the use of residue-by-residue accuracy for assessing the quality of consensus secondary structure predictions. The accuracies of recent consensus predictions for the annexins, SH2 domains and SH3 domains fall within the expected range for a perfect prediction.

Structural features can be unconserved in proteins with similar folds. An analysis of side-chain to side-chain contacts secondary structure and accessibility.

Side-chain to side-chain contacts, accessibility, secondary structure and RMS deviation were compared within 607 pairs of proteins having similar three-dimensional (3D) structures. Three types of protein 3D structural similarities were defined: type A having sequence and usually functional similarity; type B having functional, but no sequence similarity; and type C having only 3D structural similarity. Within proteins having little or no sequence similarity (types B and C), structural features frequently had a degree of conservation comparable to dissimilar 3D structures. Despite similar protein folds, as few as 30% of residues within similar protein 3D structures can form a common core. RMS deviations on core C alpha atoms can be as high as 3.2 A. Similar protein structures can have secondary structure identities as low as 41%, which is equivalent to that expected by chance. By defining three categories of amino acid accessibility (buried, half buried and exposed), some similar protein 3D structures have as few as 30% of positions in the same category, making them indistinguishable from pairs of dissimilar protein structures. Similar structures can also have as few as 12% of common side-chain to side-chain contacts, and virtually no similar energetically favourable side-chain to side-chain interactions. Complementary changes are defined as structurally equivalent pairs of interacting residues in two structures with energetically favourable but different side-chain interactions. For many proteins with similar three-dimensional structures, the proportion of complementary changes is near to that expected by chance, suggesting that many similar structures have fundamentally different stabilising interactions. All of the results suggest that proteins having similar 3D structures can have little in common apart from a scaffold of core secondary structures. This has profound implications for methods of protein fold detection, since many of the properties assumed to be conserved across similar protein 3D structures (e.g. accessibility, side-chain to side-chain contacts, etc.) are often unconserved within weakly similar (i.e. type B and C) protein 3D structures. Little difference was found between type B and C similarities suggesting that the structure of similar proteins can evolve beyond recognition even when function is conserved. Our findings suggest that it is more general features of protein structure, such as the requirements for burial of hydrophobic residues and exposure of polar residues, rather than specific residue-residue interactions that determine how well a particular sequence adopts a particular fold.(ABSTRACT TRUNCATED AT 400 WORDS)

Analysis and prediction of functional sub-types from protein sequence alignments.

The increasing number and diversity of protein sequence families requires new methods to define and predict details regarding function. Here, we present a method for analysis and prediction of functional sub-types from multiple protein sequence alignments. Given an alignment and set of proteins grouped into sub-types according to some definition of function, such as enzymatic specificity, the method identifies positions that are indicative of functional differences by comparison of sub-type specific sequence profiles, and analysis of positional entropy in the alignment. Alignment positions with significantly high positional relative entropy correlate with those known to be involved in defining sub-types for nucleotidyl cyclases, protein kinases, lactate/malate dehydrogenases and trypsin-like serine proteases. We highlight new positions for these proteins that suggest additional experiments to elucidate the basis of specificity. The method is also able to predict sub-type for unclassified sequences. We assess several variations on a prediction method, and compare them to simple sequence comparisons. For assessment, we remove close homologues to the sequence for which a prediction is to be made (by a sequence identity above a threshold). This simulates situations where a protein is known to belong to a protein family, but is not a close relative of another protein of known sub-type. Considering the four families above, and a sequence identity threshold of 30 %, our best method gives an accuracy of 96 % compared to 80 % obtained for sequence similarity and 74 % for BLAST. We describe the derivation of a set of sub-type groupings derived from an automated parsing of alignments from PFAM and the SWISSPROT database, and use this to perform a large-scale assessment. The best method gives an average accuracy of 94 % compared to 68 % for sequence similarity and 79 % for BLAST. We discuss implications for experimental design, genome annotation and the prediction of protein function and protein intra-residue distances.

Identification of distant homologues of fibroblast growth factors suggests a common ancestor for all beta-trefoil proteins.

Determination of the structures of fibroblast growth factors and interleukin-1s has previously revealed that they both adopt a beta-trefoil fold, similar to those found in Kunitz soybean trypsin inhibitors, ricin-like toxins, plant agglutinins and hisactophilin. These families possess distinct functions and occur in different subcellular localisations, and they appear to lack significant similarities in their sequences, ligands and modes of ligand binding. We have analysed the significance of sequence identities observed after structure alignment and provide statistical evidence that these beta-trefoil proteins are all homologues, having arisen from a common ancestor. In addition, we have explored the sequence space of all beta-trefoil proteins and have determined that the actin-binding proteins fascins, and other proteins of unknown function, are beta-trefoil family homologues. Unlike other beta-trefoil proteins, the triplicated repeats in each of the four beta-trefoil domains of fascins are significantly similar in sequence. This hints at how the beta-trefoil fold arose from the duplication of an ancestral gene encoding a homotrimeric single-repeat protein. The combined analysis of structure and sequence databases for detecting significant similarities is suggested as a highly sensitive approach to determining the common ancestry of extremely divergent homologues.

The serine protease inhibitor canonical loop conformation: examples found in extracellular hydrolases, toxins, cytokines and viral proteins.

Methods for the prediction of protein function from structure are of growing importance in the age of structural genomics. Here, we focus on the problem of identifying sites of potential serine protease inhibitor interactions on the surface of proteins of known structure. Given that there is no sequence conservation within canonical loops from different inhibitor families, we first compare representative loops to all fragments of equal length among proteins of known structure by calculating main-chain RMS deviation. Fragments with RMS deviation below a certain threshold (hits) are removed if residues have solvent accessibilities appreciably lower than those observed in the search structure. These remaining hits are further filtered to remove those occurring largely within secondary structure elements. Likely functional significance is restricted further by considering only extracellular protein domains. By comparing different canonical loop structures to the protein structure database, we show that the method is able to detect previously known inhibitors. In addition, we discuss potentially new canonical loop structures found in secreted hydrolases, toxins, viral proteins, cytokines and other proteins. We discuss the possible functional significance of several of the examples found, and comment on implications for the prediction of function from protein 3D structure.

Protein fold recognition by mapping predicted secondary structures.

A strategy is presented for protein fold recognition from secondary structure assignments (alpha-helix and beta-strand). The method can detect similarities between protein folds in the absence of sequence similarity. Secondary structure mapping first identifies all possible matches (maps) between a query string of secondary structures and the secondary structures of protein domains of known three-dimensional structure. The maps are then passed through a series of structural filters to remove those that do not obey simple rules of protein structure. The surviving maps are ranked by scores from the alignment of predicted and experimental accessibilities. Searches made with secondary structure assignments for a test set of 11 fold-families put the correct sequence-dissimilar fold in the first rank 8/11 times. With cross-validated predictions of secondary structure this drops to 4/11 which compares favourably with the widely used THREADER program (1/11). The structural class is correctly predicted 10/11 times by the method in contrast to 5/11 for THREADER. The new technique obtains comparable accuracy in the alignment of amino acid residues and secondary structure elements. Searches are also performed with published secondary structure predictions for the von-Willebrand factor type A domain, the proteasome 20 S alpha subunit and the phosphotyrosine interaction domain. These searches demonstrate how the method can find the correct fold for a protein from a carefully constructed secondary structure prediction, multiple sequence alignment and distant restraints. Scans with experimentally determined secondary structures and accessibility, recognise the correct fold with high alignment accuracy (86% on secondary structures). This suggests that the accuracy of mapping will improve alongside any improvements in the prediction of secondary structure or accessibility. Application to NMR structure determination is also discussed.

Supersites within superfolds. Binding site similarity in the absence of homology.

A method is presented to assess the significance of binding site similarities within superimposed protein three-dimensional (3D) structures and applied to all similar structures in the Protein Data Bank. For similarities between 3D structures lacking significant sequence similarity, the important distinction was made between remote homology (an ancient common ancestor) and analogy (likely convergence to a folding motif) according to the structural classification of proteins (SCOP) database. Supersites were defined as structural locations on groups of analogous proteins (i.e. superfolds) showing a statistically significant tendency to bind substrates despite little evidence of a common ancestor for the proteins considered. We identify three potentially new superfolds containing supersites: ferredoxin-like folds, four-helical bundles and double-stranded beta helices. In addition, the method quantifies binding site similarities within homologous proteins and previously identified supersites such as that found in the beta/alpha (TIM) barrels. For the nine superfolds, the accuracy of predictions of binding site locations is assessed. Implications for protein evolution, and the prediction of protein function either through fold recognition or tertiary structure comparison, are discussed.

Recognition of analogous and homologous protein folds: analysis of sequence and structure conservation.

An analysis was performed on 335 pairs of structurally aligned proteins derived from the structural classification of proteins (SCOP http://scop.mrc-lmb.cam.ac.uk/scop/) database. These similarities were divided into analogues, defined as proteins with similar three-dimensional structures (same SCOP fold classification) but generally with different functions and little evidence of a common ancestor (different SCOP superfamily classification). Homologues were defined as pairs of similar structures likely to be the result of evolutionary divergence (same superfamily) and were divided into remote, medium and close sub-divisions based on the percentage sequence identity. Particular attention was paid to the differences between analogues and remote homologues, since both types of similarities are generally undetectable by sequence comparison and their detection is the aim of fold recognition methods. Distributions of sequence identities and substitution matrices suggest a higher degree of sequence similarity in remote homologues than in analogues. Matrices for remote homologues show similarity to existing mutation matrices, providing some validity for their use in previously described fold recognition methods. In contrast, matrices derived from analogous proteins show little conservation of amino acid properties beyond broad conservation of hydrophobic or polar character. Secondary structure and accessibility were more conserved on average in remote homologues than in analogues, though there was no apparent difference in the root-mean-square deviation between these two types of similarities. Alignments of remote homologues and analogues show a similar number of gaps, openings (one or more sequential gaps) and inserted/deleted secondary structure elements, and both generally contain more gaps/openings/deleted secondary structure elements than medium and close homologues. These results suggest that gap parameters for fold recognition should be more lenient than those used in sequence comparison. Parameters were derived from the analogue and remote homologue datasets for potential used in fold recognition methods. Implications for protein fold recognition and evolution are discussed.

The mutational pattern of primary lymphoma of the central nervous system determined by whole-exome sequencing.

To decipher the mutational pattern of primary CNS lymphoma (PCNSL), we performed whole-exome sequencing to a median coverage of 103 × followed by mutation verification in 9 PCNSL and validation using Sanger sequencing in 22 PCNSL. We identified a median of 202 (range: 139-251) potentially somatic single nucleotide variants (SNV) and 14 small indels (range: 7-22) with potentially protein-changing features per PCNSL. Mutations affected the B-cell receptor, toll-like receptor, and NF-κB and genes involved in chromatin structure and modifications, cell-cycle regulation, and immune recognition. A median of 22.2% (range: 20.0-24.7%) of somatic SNVs in 9 PCNSL overlaps with the RGYW motif targeted by somatic hypermutation (SHM); a median of 7.9% (range: 6.2-12.6%) affects its hotspot position suggesting a major impact of SHM on PCNSL pathogenesis. In addition to the well-known targets of aberrant SHM (aSHM) (PIM1), our data suggest new targets of aSHM (KLHL14, OSBPL10, and SUSD2). Among the four most frequently mutated genes was ODZ4 showing protein-changing mutations in 4/9 PCNSL. Together with mutations affecting CSMD2, CSMD3, and PTPRD, these findings may suggest that alterations in genes having a role in CNS development may facilitate diffuse large B-cell lymphoma manifestation in the CNS. This may point to intriguing mechanisms of CNS tropism in PCNSL.

Sialidase-like Asp-boxes: sequence-similar structures within different protein folds.

Sequence similarity is the most common measure currently used to infer homology between proteins. Typically, homologous protein domains show sequence similarity over their entire lengths. Here we identify Asp box motifs, initially found as repeats in sialidases and neuraminidases, in new structural and sequence contexts. These motifs represent significantly similar sequences, localized to beta hairpins within proteins that are otherwise different in sequence and three-dimensional structure. By performing a combined sequence- and structure-based analysis we detect Asp boxes in more than nine protein families, including bacterial ribonucleases, sulfite oxidases, reelin, netrins, some lipoprotein receptors, and a variety of glycosyl hydrolases. Although the function common to each of these proteins, if any, remains unclear, we discuss possible functions of Asp boxes on the basis of previously determined experimental results and discuss different evolutionary scenarios for the origin of Asp-box containing proteins.

Conservation analysis and structure prediction of the SH2 family of phosphotyrosine binding domains.

Src homology 2 (SH2) regions are short (approximately 100 amino acids), non-catalytic domains conserved among a wide variety of proteins involved in cytoplasmic signaling induced by growth factors. It is thought that SH2 domains play an important role in the intracellular response to growth factor stimulation by binding to phosphotyrosine containing proteins. In this paper we apply the techniques of multiple sequence alignment, secondary structure prediction and conservation analysis to 67 SH2 domain amino acid sequences. This combined approach predicts seven core secondary structure regions with the pattern beta-alpha-beta-beta-beta-beta-alpha, identifies those residues most likely to be buried in the hydrophobic core of the native SH2 domain, and highlights patterns of conservation indicative of secondary structural elements. Residues likely to be involved in phosphotyrosine binding are shown and orientations of the predicted secondary structures suggested which could enable such residues to cooperate in phosphate binding. We propose a consensus pattern that encapsulates the principal conserved features of the SH2 domains. Comparison of the proposed SH2 domain of akt to this pattern shows only 12/40 matches, suggesting that this domain may not exhibit SH2-like properties.

Antimalarial activity of some novel derivatives of 2,4-diamino-5(p-chlorophenyl)-6-ethylpyrimidine (pyrimethamine).

Thirteen new analogs of 2,4-diamino-5(p-chlorophenyl)-6-ethylpyrimidine (Daraprim, pyrimethamine) in which the alph position of the 6-ethyl substituent was modified were prepared. The respective oxygens analogs (ketals, ketone, alcohol), the dimethyl hydrazone, and the nitrone displayed activities in the range of 1/4 to 1/16 that of pyrimethamine toward Plasmodium berghei in mice. The therapeutic ratios of some of these compounds may be slightly better than that of pyrimethamine.