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.

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.

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.

Exon structure conservation despite low sequence similarity: a relic of dramatic events in evolution?

The evolutionary significance of introns remains a mystery. The current availability of several complete eukaryotic genomes permits new studies to probe the possible function of these peculiar genomic features. Here we investigate the degree to which gene structure (intron position, phase and length) is conserved between homologous protein domains. We find that for certain extracellular-signalling and nuclear domains, gene structures are similar even when protein sequence similarity is low or not significant and sequences can only be aligned with a knowledge of protein tertiary structure. In contrast, other domains, including most intracellular signalling modules, show little gene structure conservation. Intriguingly, many domains with conserved gene structures, such as cytokines, are involved in similar biological processes, such as the immune response. This suggests that gene structure conservation may be a record of key events in evolution, such as the origin of the vertebrate immune system or the duplication of nuclear receptors in nematodes. The results suggest ways to detect new and potentially very remote homologues, and to construct phylogenies for proteins with limited sequence similarity.

On the evolution of protein folds: are similar motifs in different protein folds the result of convergence, insertion, or relics of an ancient peptide world?

This paper presents and discusses evidence suggesting how the diversity of domain folds in existence today might have evolved from peptide ancestors. We apply a structure similarity detection method to detect instances where localized regions of different protein folds contain highly similar sequences and structures. Results of performing an all-on-all comparison of known structures are described and compared with other recently published findings. The numerous instances of local sequence and structure similarities within different protein folds, together with evidence from proteins containing sequence and structure repeats, argues in favor of the evolution of modern single polypeptide domains from ancient short peptide ancestors (antecedent domain segments (ADSs)). In this model, ancient protein structures were formed by self-assembling aggregates of short polypeptides. Subsequently, and perhaps concomitantly with the evolution of higher fidelity DNA replication and repair systems, single polypeptide domains arose from the fusion of ADSs genes. Thus modern protein domains may have a polyphyletic origin.

Adeno-associated virus is associated with a lower risk of high-grade cervical neoplasia.

Adeno-associated virus (AAV) is a ubiquitous human helper-dependent parvovirus which may interact with human papillomaviruses (HPV) to modify a woman's risk of cervical neoplasia. This analysis was nested in a cohort study of low-income women receiving Pap smears as part of their family planning services. We selected cases (55 with high-grade cervical squamous intraepithelial lesions (HSIL) and 162 with low-grade LSIL) and controls (96 women with normal cervical cytology) and analyzed cervical DNA for AAV, using PCR amplification/dot blot hybridization, and HPV, using hybrid capture I. AAV positivity was associated with a significantly reduced risk of HSIL (age and HPV-adjusted odds ratio (aOR) = 0.32) yet not with LSIL (aOR = 0.78); 53.8% of HSIL, 66.9% of LSIL, and 70.7% of controls were AAV+. AAV appears to interact with HPV to reduce SIL risk; relative to the HPV-/AAV+ exposure, the respective aORs for HSIL and HPV+/AAV-, HPV+/AAV+, and HPV-/AAV+ were 17.0, 6.9, and 3.5. AAV+ was not associated with age, race, HPV status, or sexual or reproductive risk factors. These results strongly suggest that AAV may play a protective or inhibitory role in late stage cervical carcinogenesis. This conclusion needs to be verified in additional epidemiologic studies.

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.

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.

Predicting function from structure: examples of the serine protease inhibitor canonical loop conformation found in extracellular proteins.

The prediction of protein function from structure is becoming of growing importance in the age of structural genomics. We have focused 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. Also a test is performed to see if the loop can dock into the binding site of the serine protease trypsin without unacceptable steric clashes. By comparing different canonical loop structures to the protein structure database we show that the method was 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.

Fold recognition from sequence comparisons.

We applied a new protocol based on PSI-Blast to predict the structures of fold recognition targets during CASP4. The protocol used a back-validation step to infer biologically significant connections between sequences with PSI-Blast E-values up to 10. If connections were found to proteins of known structure, alignments were generated by using HMMer. The protocol was implemented in a fully automated version (SBauto) and in a version that allowed manual intervention (SBfold). We found that the automated version made 17 predictions for target domains, of which 8 identified the correct fold with an average alignment accuracy of 24% for alignable residues and 43% for equivalent secondary structure elements. The manual version improved predictions somewhat, with 10 of 15 predictions identifying the correct fold with alignment accuracies of 33% for alignable residues and 64% for equivalent secondary structure elements. We describe successes and failures of our approach and discuss future developments of fold recognition.

Changes in mitochondrial genetic codes as phylogenetic characters: two examples from the flatworms.

Shared molecular genetic characteristics other than DNA and protein sequences can provide excellent sources of phylogenetic information, particularly if they are complex and rare and are consequently unlikely to have arisen by chance convergence. We have used two such characters, arising from changes in mitochondrial genetic code, to define a clade within the Platyhelminthes (flatworms), the Rhabditophora. We have sampled 10 distinct classes within the Rhabditophora and find that all have the codon AAA coding for the amino acid Asn rather than the usual Lys and AUA for Ile rather than the usual Met. We find no evidence to support claims that the codon UAA codes for Tyr in the Platyhelminthes rather than the standard stop codon. The Rhabditophora are a very diverse group comprising the majority of the free-living turbellarian taxa and the parasitic Neodermata. In contrast, three other classes of turbellarian flatworm, the Acoela, Nemertodermatida, and Catenulida, have the standard invertebrate assignments for these codons and so are convincingly excluded from the rhabditophoran clade. We have developed a rapid computerized method for analyzing genetic codes and demonstrate the wide phylogenetic distribution of the standard invertebrate code as well as confirming already known metazoan deviations from it (ascidian, vertebrate, echinoderm/hemichordate).

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.

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.

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.

Using focus groups to develop a culturally sensitive videotape intervention for HIV-positive women.

Research-based interventions for women with the human immunodeficiency virus (HIV) are usually developed without input from the women who receive the intervention. An exploratory study was performed using focus group methodology to develop a culturally sensitive videotape intervention for educating HIV-positive women about pregnancy and antiretroviral use. Women who met the study criteria were HIV-positive and of childbearing age. These women volunteered to participate in the focus groups to provide information on decisions concerning pregnancy and antiretroviral use during pregnancy to decrease perinatal transmission. A total of five focus groups were conducted in 1998. Responses to three questions that were relevant to the video are presented in this article. Information gained from the focus groups was used successfully to develop a videotape currently being used in a multisite intervention study. Focus group methodology is a useful strategy to develop culturally and content relevant educational interventions for research and practice.

Fold recognition using sequence and secondary structure information.

We applied a succession of sequence search and structure prediction methods to the targets in the fold recognition part of the CASP3 experiment. For each target, we expanded an initial sequence space, obtained through PSI-BLAST, by searching for statistically significant relationships to low-scoring sequences and then by searching for conserved sequence patterns. We then divided the proteins in the sequence space into families and built an alignment hierarchically, using the multiple alignment program MACAW. If no significant similarity to a protein of known structure was apparent at this point, we submitted the alignment to the Jpred server for consensus secondary structure prediction and searched the structure space using the secondary structure mapping program MAP. Failing this, we compared the structural properties that we believed we recognized in the aligned proteins to the folds in the SCOP database, using visual inspection. If all these methods failed to uncover a plausible match, we predicted that the target would adopt a novel fold. This procedure yielded correct answers for seven of twenty-one targets and a partly correct answer for one. A retrospective analysis shows that automating the sequence search procedures would have represented a significant improvement, with at least three additional correct predictions.

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.