NBLAST: a cluster variant of BLAST for NxN comparisons

Background  

The BLAST algorithm compares biological sequences to one another in order to determine shared motifs and common ancestry. However, the comparison of all non-redundant (NR) sequences against all other NR sequences is a computationally intensive task. We developed NBLAST as a cluster computer implementation of the BLAST family of sequence comparison programs for the purpose of generating pre-computed BLAST alignments and neighbour lists of NR sequences.     

Physical descriptions of experimental selectivity measurements in ion channels

Three experiments that quantify the amount of selectivity exhibited by a biological ion channel are examined with Poisson-Nernst-Planck (PNP) theory. Conductance ratios and the conductance mole fraction experiments are examined by considering a simple model ion channel for which an approximate solution to the PNP equations with Donnan boundary conditions is derived. A more general result is derived for the Goldman-Hodgkin-Katz permeability ratio. The results show that (1) the conductance ratio measures the ratio of the diffusion coefficients of the ions inside the channel, (2) the mole fraction experiment measures the difference of the excess chemical potentials of the ions inside the channel, and (3) the permeability ratio measures both diffusion coefficients and excess chemical potentials. The results are used to divide selectivity into two components: partitioning, an equilibrium measure of how well the ions enter the channel, and diffusion, a nonequilibrium measure of how well the ions move through the channel.     

GXXXG and AXXXA: common alpha-helical interaction motifs in proteins, particularly in extremophiles

The GXXXG motif is a frequently occurring sequence of residues that is known to favor helix-helix interactions in membrane proteins. Here we show that the GXXXG motif is also prevalent in soluble proteins whose structures have been determined. Some 152 proteins from a non-redundant PDB set contain at least one alpha-helix with the GXXXG motif, 41 +/- 9% more than expected if glycine residues were uniformly distributed in those alpha-helices. More than 50% of the GXXXG-containing alpha-helices participate in helix-helix interactions. In fact, 26 of those helix-helix interactions are structurally similar to the helix-helix interaction of the glycophorin A dimer, where two transmembrane helices associate to form a dimer stabilized by the GXXXG motif. As for the glycophorin A structure, we find backbone-to-backbone atomic contacts of the C alpha-H...O type in each of these 26 helix-helix interactions that display the stereochemical hallmarks of hydrogen bond formation. These glycophorin A-like helix-helix interactions are enriched in the general set of helix-helix interactions containing the GXXXG motif, suggesting that the inferred C alpha-H...O hydrogen bonds stabilize the helix-helix interactions. In addition to the GXXXG motif, some 808 proteins from the non-redundant PDB set contain at least one alpha-helix with the AXXXA motif (30 +/- 3% greater than expected). Both the GXXXG and AXXXA motifs occur frequently in predicted alpha-helices from 24 fully sequenced genomes. Occurrence of the AXXXA motif is enhanced to a greater extent in thermophiles than in mesophiles, suggesting that helical interaction based on the AXXXA motif may be a common mechanism of thermostability in protein structures. We conclude that the GXXXG sequence motif stabilizes helix-helix interactions in proteins, and that the AXXXA sequence motif also stabilizes the folded state of proteins.     

Multicopy crystallographic refinement of a relaxed glutamine synthetase from Mycobacterium tuberculosis highlights flexible loops in the enzymatic mechanism and its regulation

The crystal structure of glutamine synthetase (GS) from Mycobacterium tuberculosis determined at 2.4 A resolution reveals citrate and AMP bound in the active site. The structure was refined with strict 24-fold noncrystallographic symmetry (NCS) constraints and has an R-factor of 22.7% and an R-free of 25.5%. Multicopy refinement using 10 atomic models and strict 24-fold NCS constraints further reduced the R-factor to 20.4% and the R-free to 23.2%. The multicopy model demonstrates the range of atomic displacements of catalytic and regulatory loops in glutamine synthesis, simulating loop motions. A comparison with loop positions in substrate complexes of GS from Salmonella typhimurium shows that the Asp50 and Glu327 loops close over the active site during catalysis. These loop closures are preceded by a conformational change of the Glu209 beta-strand upon metal ion or ATP binding that converts the enzyme from a relaxed to a taut state. We propose a model of the GS regulatory mechanism based on the loop motions in which adenylylation of the Tyr397 loop reverses the effect of metal ion binding, and regulates intermediate formation by preventing closure of the Glu327 loop.     

Computational methods of analysis of protein-protein interactions

Computational methods play an important role at all stages of the process of determining protein-protein interactions. They are used to predict potential interactions, to validate the results of high-throughput interaction screens and to analyze the protein networks inferred from interaction databases.     

MHC class II-bound self peptides from autoimmune MRL/lpr mice reveal potential T cell epitopes for autoantibody production in murine systemic lupus erythematosus

The systemic lupus erythematosus-like syndrome in MRL/lpr mice involves high-titered IgG autoantibodies, particularly antinuclear Abs that target histones, DNA, and RNA particles. Although T cell help is required for the generation of antinuclear Abs, the epitopes recognized by such helper T cells are unknown. To address this question, we isolated and sequenced self peptides bound by MHC class II molecules from MRL/lpr mice. We identified a number of peptides that are not seen in similar preparations from nonautoimmune C3H animals. The "abnormal" peptide donors include histone, a protein component of a small nuclear ribonucleoprotein, ribosomal proteins, and RNA processing enzymes. We postulate that the peptides from these donors are T cell epitopes required for the generation of the most frequent antinuclear Abs specificities seen in MRL/lpr mice.     

DIP: the database of interacting proteins

The Database of Interacting Proteins (DIP; http://dip.doe-mbi.ucla.edu) is a database that documents experimentally determined protein-protein interactions. This database is intended to provide the scientific community with a comprehensive and integrated tool for browsing and efficiently extracting information about protein interactions and interaction networks in biological processes. Beyond cataloging details of protein-protein interactions, the DIP is useful for understanding protein function and protein-protein relationships, studying the properties of networks of interacting proteins, benchmarking predictions of protein-protein interactions, and studying the evolution of protein-protein interactions.     

Role of cyclin G-associated kinase in uncoating clathrin-coated vesicles from non-neuronal cells

Auxilin is a brain-specific DnaJ homolog that is required for Hsc70 to dissociate clathrin from bovine brain clathrin-coated vesicles. However, Hsc70 is also involved in uncoating clathrin-coated vesicles formed at the plasma membrane of non-neuronal cells suggesting that an auxilin homolog may be required for uncoating in these cells. One candidate is cyclin G-associated kinase (GAK), a 150-kDa protein expressed ubiquitously in various tissues. GAK has a C-terminal domain with high sequence similarity to auxilin; like auxilin this C-terminal domain consists of three subdomains, an N-terminal tensin-like domain, a clathrin-binding domain, and a C-terminal J-domain. Western blot analysis shows that GAK is present in rat liver, bovine testes, and bovine brain clathrin-coated vesicles. More importantly, liver clathrin-coated vesicles, which contain GAK but not auxilin, are uncoated by Hsc70, suggesting that GAK acts as an auxilin homolog in non-neuronal cells. In support of this view, the clathrin-binding domain of GAK alone induces clathrin polymerization into baskets and the combined clathrin-binding domain and J-domain of GAK supports uncoating of AP180-clathrin baskets by Hsc70 at pH 7 and induces Hsc70 binding to clathrin baskets at pH 6. Immunolocalization studies suggest that GAK is a cytosolic protein that is concentrated in the perinuclear region; it appears to be highly associated with the trans-Golgi where the budding of clathrin-coated vesicles occurs. We propose that GAK is a required cofactor for the uncoating of clathrin-coated vesicles by Hsc70 in non-neuronal cells.     

Heregulin reverses the oligomerization of HER3

We analyzed the propensity of the HER3 receptor and its extracellular domain (ECD) to undergo ligand-independent self-association. The HER3-ECD, purified from Drosophila S2 cells, binds the EGF-like domain of heregulin (hrg) with a K(d) of 1.9 nM as measured by surface plasmon resonance (SPR) studies. In a gel shift assay, the HER3-ECD self-associates into a uniform, slowly migrating species in a concentration-dependent manner, starting at concentrations of <10 nM. In contrast to the HER3-ECD, the ECD from the related HER2 receptor does not oligomerize under the same conditions. The direct interaction of HER3-ECDs was also demonstrated by pull-down assays and SPR measurements under physiological salt conditions. This self-association of the HER3-ECD was reversed by the addition of hrg but not by EGF. The apparent equilibrium dissociation constant for the HER3-ECD self-association is 15 nM, based on SPR measurements. In this analysis, hrg blocks HER3-ECD self-association, and the addition of hrg during the dissociation phase resulted in an accelerated off rate. This finding suggests that hrg can bind to and disrupt preexisting HER3-ECD oligomers. Full-length HER3 likewise exhibited self-association. Under conditions where co-immunoprecipitation and cross-linking of HER2 and HER3 were stimulated by hrg, HER3 self-association and cross-linking were disrupted by hrg. The implication is that the self-association of HER3-ECD favors the formation of catalytically inactive complexes of the HER3 receptor. Binding of hrg releases HER3 which may then form signaling-competent HER3-HER2 heterodimers.