The GRID: The General Repository for Interaction Datasets

We have developed a relational database, called the General Repository for Interaction Datasets (The GRID) to archive and display physical, genetic and functional interactions. The GRID displays data-rich interaction tables for any protein of interest, combines literature-derived and high-throughput interaction datasets, and is readily accessible via the web. Interactions parsed in The GRID can be viewed in graphical form with a versatile visualization tool called Osprey.     

GRID: a novel Grb-2-related adapter protein that interacts with the activated T cell costimulatory receptor CD28

Adapter proteins such as Grb2 play a central role in the formation of signaling complexes through their association with multiple protein binding partners. These interactions are mediated by specialized domains such as the well-characterized Src homology SH2 and SH3 motifs. Using yeast three-hybrid technology, we have identified a novel adapter protein, expressed predominantly in T lymphocytes, that associates with the activated form of the costimulatory receptor, CD28. The protein is a member of the Grb2 family of adapter proteins and contains an SH3-SH2-SH3 domain structure. A unique glutamine/proline-rich domain (insert domain) of unknown function is situated between the SH2 and N-terminal SH3 domains. We term this protein GRID for Grb2-related protein with insert domain. GRID coimmunoprecipitates with CD28 from Jurkat cell lysates following activation of CD28. Using mutants of CD28 and GRID, we demonstrate that interaction between the proteins is dependent on phosphorylation of CD28 at tyrosine 173 and integrity of the GRID SH2 domain, although there are also subsidiary stabilizing contacts between the PXXP motifs of CD28 and the GRID C-terminal SH3 domain. In addition to CD28, GRID interacts with a number of other T cell signaling proteins, including SLP-76 (SH2 domain-containing leukocyte protein of 76 kDa), p62dok, and RACK-1 (receptor for activated protein kinase C-1). These findings suggest that GRID functions as an adapter protein in the CD28-mediated costimulatory pathway in T cells.     

ClusFCM: an algorithm for predicting protein functions using homologies and protein interactions

We introduce a new algorithm, called ClusFCM, which combines techniques of clustering and fuzzy cognitive maps (FCM) for prediction of protein functions. ClusFCM takes advantage of protein homologies and protein interaction network topology to improve low recall predictions associated with existing prediction methods. ClusFCM exploits the fact that proteins of known function tend to cluster together and deduce functions not only through their direct interaction with other proteins, but also from other proteins in the network. We use ClusFCM to annotate protein functions for Saccharomyces cerevisiae (yeast), Caenorhabditis elegans (worm), and Drosophila melanogaster (fly) using protein-protein interaction data from the General Repository for Interaction Datasets (GRID) database and functional labels from Gene Ontology (GO) terms. The algorithm's performance is compared with four state-of-the-art methods for function prediction--Majority, chi(2) statistics, Markov random field (MRF), and FunctionalFlow--using measures of Matthews correlation coefficient, harmonic mean, and area under the receiver operating characteristic (ROC) curves. The results indicate that ClusFCM predicts protein functions with high recall while not lowering precision. Supplementary information is available at http://www.egr.vcu.edu/cs/dmb/ClusFCM/.     

Spatially fractionated radiation induces cytotoxicity and changes in gene expression in bystander and radiation adjacent murine carcinoma cells

Radiation-induced bystander effects have been extensively studied at low doses, since evidence of bystander induced cell killing and other effects on unirradiated cells were found to be predominant at doses up to 0.5 Gy. Therefore, few studies have examined bystander effects induced by exposure to higher doses of radiation, such as spatially fractionated radiation (GRID) treatment. In the present study, we evaluate the ability of GRID treatment to induce changes in GRID adjacent (bystander) regions, in two different murine carcinoma cell lines following exposure to a single irradiation dose of 10 Gy. Murine SCK mammary carcinoma cells and SCCVII squamous carcinoma cells were irradiated using a brass collimator to create a GRID pattern of nine circular fields 12 mm in diameter with a center-to-center distance of 18 mm. Similar to the typical clinical implementation of GRID, this is approximately a 50:50 ratio of direct and bystander exposure. We also performed experiments by irradiating separate cultures and transferring the medium to unirradiated bystander cultures. Clonogenic survival was evaluated in both cell lines to determine the occurrence of radiation-induced bystander effects. For the purpose of our study, we have defined bystander cells as GRID adjacent cells that received approximately 1 Gy scatter dose or unirradiated cells receiving conditioned medium from irradiated cells. We observed significant bystander killing of cells adjacent to the GRID irradiated regions compared to sham treated controls. We also observed bystander killing of SCK and SCCVII cells cultured in conditioned medium obtained from cells irradiated with 10 Gy. Therefore, our results confirm the occurrence of bystander effects following exposure to a high-dose of radiation and suggest that cell-to-cell contact is not required for these effects. In addition, the gene expression profile for DNA damage and cellular stress response signaling in SCCVII cells after GRID exposure was studied. The occurrence of GRID-induced bystander gene expression changes in significant numbers of DNA damage and cellular stress response signaling genes, providing molecular evidence for possible mechanisms of bystander cell killing.     

The eNMR platform for structural biology

The e-NMR project is a European cooperation initiative that aims at providing the bio-NMR user community with a software platform integrating and streamlining the computational approaches necessary for the analysis of bio-NMR data. The e-NMR platform is based on a Grid computational infrastructure. A main focus of the current implementation of the e-NMR platform is on streamlining structure determination protocols. Indeed, to facilitate the use of NMR spectroscopy in the life sciences, the eNMR consortium has set out to provide protocolized services through easy-to-use web interfaces, while still retaining sufficient flexibility to handle specific requests by expert users. Various programs relevant for structural biology applications are already available through the e-NMR portal, including HADDOCK, XPLOR-NIH, CYANA and csRosetta. The implementation of these services, and in particular the distribution of calculations to the GRID infrastructure, has required the development of specific tools. However, the GRID infrastructure is maintained completely transparent to the users. With more than 150 registered users, eNMR is currently the second largest European Virtual Organization in the life sciences.     

An innovative application of the "flexible" GRID/PCA computational method: study of differences in selectivity between PGAs from Escherichia coli and a Providentia rettgeri mutant

The original GRID/PCA technique was adapted for the development of a tool potentially useful for the plan of a research strategy in rational enzyme design. The use of the MOVE directive of GRID made it possible to partially take into account protein flexibility, and the multivariate analysis was used as an instrument for focusing only on relevant information related to the differences in enzyme substrate selectivities. The comparison of two different penicillin G acylases, from Escherichia coli and from Providentia rettgeri, was used as a case study; these enzymes are very similar and their reported selectivities differ only for a couple of mutations around the active site. The "flexible" GRID/PCA method was able to correctly predict the observed selectivity differences caused not only by mutations of residues of the active site but also by long range effects on substrate selectivity due to sequence mutations on residues not directly involved in substrate recognition.     

Lurcher GRID2-induced death and depolarization can be dissociated in cerebellar Purkinje cells

The Lurcher mutation transforms the GRID2 receptor into a constitutively opened channel. In Lurcher heterozygous mice, cerebellar Purkinje cells are permanently depolarized, a characteristic that has been thought to be the primary cause of their death, which occurs from the second postnatal week onward. The more dramatic phenotype of Lurcher homozygotes is thought to be due to a simple gene dosage effect of the mutant allele. We have analyzed the phenotype of Lurcher/hotfoot heteroallelic mutants bearing only one copy of the Lurcher allele and no wild-type Grid2. Our results show that the absence of wild-type GRID2 receptors in these heteroallelic mutants induces an early and massive Purkinje cell death that is correlated with early signs of autophagy. This neuronal death is independent of depolarization and can be explained by the direct activation of autophagy by Lurcher GRID2 receptors through the recently discovered signaling pathway formed by GRID2, n-PIST, and Beclin1.     

Substrate-enzyme interactions and catalytic mechanism in phospholipase C: a molecular modeling study using the GRID program.

Based on the high-resolution X-ray crystallographic structure of phospholipase C from Bacillus cereus, the orientation of the phosphatidylcholine substrate in the active site of the enzyme is proposed. The proposal is based on extensive calculations using the GRID program and molecular mechanics geometry relaxations. The substrate model has been constructed by successively placing phosphate, choline and diacylglycerol moieties in the positions indicated from GRID calculations. On the basis of the resulting orientation of a complete phosphatidylcholine molecule, we propose a mechanism for the hydrolysis of the substrate.     

In silico screening for antibiotic escort molecules to overcome efflux

Resistance to antibiotics is a growing problem worldwide and occurs in part due to the overexpression of efflux pumps responsible for the removal of antibiotics from bacterial cells. The current study examines complex formation between efflux pump substrates and escort molecules as a criterion for an in silico screening method for molecules that are able to potentiate antibiotic activities. Initially, the SUPERDRUG database was queried to select molecules that were similar to known multidrug resistance (MDR) modulators. Molecular interaction fields generated by GRID and the docking module GLUE were used to calculate the interaction energies between the selected molecules and the antibiotic norfloxacin. Ten compounds forming the most stable complexes with favourable changes to the norfloxacin molecular properties were tested for their potentiation ability by efflux pump modulation assays. Encouragingly, two molecules were proven to act as efflux pump modulators, and hence provide evidence that complex formation between a substrate and a drug can be used for in silico screening for novel escort molecules.     

Binding affinity prediction of novel estrogen receptor ligands using receptor-based 3-D QSAR methods

We have recently reported the development of a 3-D QSAR model for estrogen receptor ligands showing a significant correlation between calculated molecular interaction fields and experimentally measured binding affinity. The ligand alignment obtained from docking simulations was taken as basis for a comparative field analysis applying the GRID/GOLPE program. Using the interaction field derived with a water probe and applying the smart region definition (SRD) variable selection procedure, a significant and robust model was obtained (q²_LOO=0.921, SDEP=0.345). To further analyze the robustness and the predictivity of the established model several recently developed estrogen receptor ligands were selected as external test set. An excellent agreement between predicted and experimental binding data was obtained indicated by an external SDEP of 0.531. Two other traditionally used prediction techniques were applied in order to check the performance of the receptor-based 3-D QSAR procedure. The interaction energies calculated on the basis of receptor–ligand complexes were correlated with experimentally observed affinities. Also ligand-based 3-D QSAR models were generated using program FlexS. The interaction energy-based model, as well as the ligand-based 3-D QSAR models yielded models with lower predictivity. The comparison with the interaction energy-based model and with the ligand-based 3-D QSAR models, respectively, indicates that the combination of receptor-based and 3-D QSAR methods is able to improve the quality of prediction.     

Path PPI: an integrated dataset of human pathways and protein-protein interactions

Integration of pathway and protein-protein interaction(PPI) data can provide more information that could lead to new biological insights. PPIs are usually represented by a simple binary model, whereas pathways are represented by more complicated models. We developed a series of rules for transforming protein interactions from pathway to binary model, and the protein interactions from seven pathway databases, including PID, Bio Carta, Reactome, Net Path, INOH, SPIKE and KEGG, were transformed based on these rules. These pathway-derived binary protein interactions were integrated with PPIs from other five PPI databases including HPRD, Int Act, Bio GRID, MINT and DIP, to develop integrated dataset(named Path PPI). More detailed interaction type and modification information on protein interactions can be preserved in Path PPI than other existing datasets. Comparison analysis results indicate that most of the interaction overlaps values(OAB) among these pathway databases were less than 5%, and these databases must be used conjunctively. The Path PPI data was provided at http://proteomeview. hupo.org.cn/Path PPI/Path PPI.html.     

Study on improving the selectivity of compounds that inhibit two PI3Ks (gamma and delta)

A desirable characteristic of PI3K inhibitors is their selectivity. Up to now, there has been no report that describes the 3 D-structure differences between two PI3Ks (δ and γ) and applies them to designing selective compounds. In the present study, we used an approach combining protein-structure modeling, GRID/PCA (Principal Component Analysis) and docking methods to investigate the detail interactions of the two PI3Ks with various chemical groups. At first, we constructed a 3 D-model of the PI3Kδ catalytic subunit with the program Modeller7.0 based on the high resolution X-ray structure of the PI3Kγ catalytic subunit, and then employed GRID and PCA to reveal the most relevant structural and physicochemical differences between the two PI3Ks related to their selectivity. As a result, the analysis unveiled the most important regions on the two PI3Ks that should be taken into account for the design of selective inhibitors. Finally, based on activity data of 10 PI3Kδ-selective compounds, a docking study validated the results of the GRID/PCA method, which suggested that the approach could provide clear guidelines for selective drug design.     

Experimental and computational mapping of the binding surface of a crystalline protein

Multiple Solvent Crystal Structures (MSCS) is a crystallographic technique to identify energetically favorable positions and orientations of small organic molecules on the surface of proteins. We determined the high-resolution crystal structures of thermolysin (TLN), generated from crystals soaked in 50--70% acetone, 50--80% acetonitrile and 50 mM phenol. The structures of the protein in the aqueous-organic mixtures are essentially the same as the native enzyme and a number of solvent interaction sites were identified. The distribution of probe molecules shows clusters in the main specificity pocket of the active site and a buried subsite. Within the active site, we compared the experimentally determined solvent positions with predictions from two computational functional group mapping techniques, GRID and Multiple Copy Simultaneous Search (MCSS). The experimentally determined small molecule positions are consistent with the structures of known protein--ligand complexes of TLN.     

The hemibiotrophic lifestyle of Colletotrichum species

Colletotrichum species infect several economically important crop plants. To establish a compatible parasitic interaction, a specialized infection cell, the melanized appressorium, is differentiated on the cuticle of the host. After penetration, an infection vesicle and primary hyphae are formed. These structures do not kill the host cell and show some similarities with haustoria formed by powdery mildews and rust fungi. Therefore, this stage of infection is called biotrophic. Later in the infection process, necrotrophic secondary hyphae spread within and kill the host tissue. The lifestyle of Colletotrichum species is called hemibiotrophic, as biotrophic and necrotrophic developmental stages are sequentially established. As most Colletotrichum species are accessible to molecular techniques, genes can be identified and functionally characterized. Here we demonstrate that Agrobacterium tumefaciens-mediated transformation is a well-suited method for tagging of genes mediating compatibility in the Colletotrichum graminicola-maize interaction.     

Bipolar I disorder and schizophrenia: a 440-single-nucleotide polymorphism screen of 64 candidate genes among Ashkenazi Jewish case-parent trios

Bipolar, schizophrenia, and schizoaffective disorders are common, highly heritable psychiatric disorders, for which familial coaggregation, as well as epidemiological and genetic evidence, suggests overlapping etiologies. No definitive susceptibility genes have yet been identified for any of these disorders. Genetic heterogeneity, combined with phenotypic imprecision and poor marker coverage, has contributed to the difficulty in defining risk variants. We focused on families of Ashkenazi Jewish descent, to reduce genetic heterogeneity, and, as a precursor to genomewide association studies, we undertook a single-nucleotide polymorphism (SNP) genotyping screen of 64 candidate genes (440 SNPs) chosen on the basis of previous linkage or of association and/or biological relevance. We genotyped an average of 6.9 SNPs per gene, with an average density of 1 SNP per 11.9 kb in 323 bipolar I disorder and 274 schizophrenia or schizoaffective Ashkenazi case-parent trios. Using single-SNP and haplotype-based transmission/disequilibrium tests, we ranked genes on the basis of strength of association (P<.01). Six genes (DAO, GRM3, GRM4, GRIN2B, IL2RB, and TUBA8) met this criterion for bipolar I disorder; only DAO has been previously associated with bipolar disorder. Six genes (RGS4, SCA1, GRM4, DPYSL2, NOS1, and GRID1) met this criterion for schizophrenia or schizoaffective disorder; five replicate previous associations, and one, GRID1, shows a novel association with schizophrenia. In addition, six genes (DPYSL2, DTNBP1, G30/G72, GRID1, GRM4, and NOS1) showed overlapping suggestive evidence of association in both disorders. These results may help to prioritize candidate genes for future study from among the many suspected/proposed for schizophrenia and bipolar disorders. They provide further support for shared genetic susceptibility between these two disorders that involve glutamate-signaling pathways.     

Knowledge-based modeling of a legume lectin and docking of the carbohydrate ligand: The Ulex europaeus lectin I and its interaction with fucose

Ulex europaeus isolectin I is specific for fucose-containing oligosaccharide such as H type 2 trisaccharide α-l-Fuc (1→2) β-d-Gal (1→4) β-d-GlcNAc. Several legume lectins have been crystallized and modeled, but no structural data are available concerning such fucose-binding lectin. The three-dimensional structure of Ulex europaeus isolectin I has been constructed using seven legume lectins for which high-resolution crystal structures were available. Some conserved water molecules, as well as the structural cations, were taken into account for building the model. In the predicted binding site, the most probable locations of the secondary hydroxyl groups were determined using the GRID method. Several possible orientations could be determined for a fucose residue. All of the four possible conformations compatible with energy calculations display several hydrogen bonds with Asp-87 and Ser-132 and a stacking interaction with Tyr-220 and Phe-136. In two orientations, the O-3 and O-4 hydroxyl groups of fucose are the most buried ones, whereas two other, the O-2 and O-3 hydroxyl groups are at the bottom of the site. Possible docking modes are also studied by analysis of the hydrophobic and hydrophilic surfaces for both the ligand and the protein. The SCORE method allows for a quantitative evaluation of the complementarity of these surfaces, on the basis of molecular lipophilicity calculations. The predictions presented here are compared with known biochemical data.