Biomphalaria snails are instrumental in transmission of the human blood fluke Schistosoma mansoni. With the World Health Organization's goal to eliminate schistosomiasis as a global health problem by 2025, there is now renewed emphasis on snail control. Here, we characterize the genome of Biomphalaria glabrata, a lophotrochozoan protostome, and provide timely and important information on snail biology. We describe aspects of phero-perception, stress responses, immune function and regulation of gene expression that support the persistence of B. glabrata in the field and may define this species as a suitable snail host for S. mansoni. We identify several potential targets for developing novel control measures aimed at reducing snail-mediated transmission of schistosomiasis.

Motivation: Identifying protein enzymatic or pharmacological activities are important areas of research in biology and chemistry. Biological and chemical databases are increasingly being populated with linkages between protein sequences and chemical structures.There is now sufficient information to apply machine-learning techniques to predict interactions between chemicals and proteins at a genome scale. Current machine-learning techniques use as input either protein sequences and structures or chemical information. We propose here a method to infer protein-chemical interactions using heterogeneous input consisting of both protein sequence and chemical information. Results: Our method relies on expressing proteins and chemicals with a common cheminformatics representation. We demonstrate our approach by predicting whether proteins can catalyze reactions not present in training sets. We also predict whether a given drug can bind a target, in the absence of prior binding information for that drug and target. Such predictions cannot be made with current machine-learning techniques requiring binding information for individual reactions or individual targets. Availability and Contact: For questions, paper reprints, please contact Jean-Loup Faulon at jfaulon@sandia.gov. Additional information on the signature molecular descriptor and codes can be downloaded at:

Quantitative structure-property relationships (QSPRs) have been developed and assessed for predicting the reorganization energy of polycyclic aromatic hydrocarbons (PAHs). Preliminary QSPR models, based on a combination of molecular signature and electronic eigenvalue difference descriptors, have been trained using more than 200 PAHs. Monte Carlo cross-validation systematically improves the performance of the models through progressive reduction of the training set and selection of best performing training subsets. The final biased QSPR model yields correlation coefficients q(2) and r(2) of 0.7 and 0.8, respectively, and an estimated error in predicting reorganization energy of ±0.014 eV.

Several pathogenic bacteria produce adenylyl cyclase toxins, such as the edema factor (EF) of Bacillus anthracis. These disturb cellular metabolism by catalyzing production of excessive amounts of the regulatory molecule cAMP. Here, a structure-based method, where a 3D-pharmacophore that fit the active site of EF was constructed from fragments, was used to identify non-nucleotide inhibitors of EF. A library of small molecule fragments was docked to the EF-active site in existing crystal structures and those with the highest HINT scores were assembled into a 3D-pharmacophore. About 10,000 compounds, from over 2.7 million compounds in the ZINC database, had a similar molecular framework. These were ranked according to their docking scores, using methodology that was shown to achieve maximum accuracy (i.e., how well the docked position matched the experimentally determined site for ATP analogues in crystal structures of the complex). Finally, 19 diverse compounds with the best AutoDock binding/docking scores were assayed in a cell based assay for their ability to reduce cAMP secretion induced by EF. Four of the test compounds, from different structural groups, inhibited in the low micromolar range. One of these has a core structure common to phosphatase inhibitors previously identified by high-throughput assays of a diversity library. Thus, the fragment based pharmacophore identified a small number of diverse compounds for assay, and greatly enhanced the selection process of advanced lead compounds for combinatorial design.

No abstract

This paper describes a novel clustering methodology for classifying over 700 conformations of a flexible analogue of GBR 12909, a dopamine reuptake inhibitor that has completed phase I clinical trials as a treatment for cocaine abuse. The major aspect of the clustering methodology includes an efficient data-conditioning scheme where a systematic feature extraction procedure based on the structural properties of the molecule was used to reduce the associated feature space. This allowed region-specific clustering that focused on individual pharmacophore elements of the molecule. For clustering of the reduced feature set, the fuzzy clustering partitional method was utilized. Due to the relational nature of the feature data, fuzzy relational clustering was employed, and it successfully detected natural groups defined by rotational minima around N(sp(3))-C(sp(3)), O(sp(3))-C(sp(3)), and C(sp(3))-C(sp(2)) bonds. The proposed clustering methodology also employed several cluster validity measures, which corroborated the partitions produced by the clustering technique and agreed with the results of hierarchical clustering using the XCluster program. Representative structures which exhibited a reasonable spread of energies and showed good spatial coverage of the conformational space were identified for use as putative bioactive conformations in a future Comparative Molecular Field Analysis of GBR 12909 analogues. The clustering methodology developed here is capable of handling other computational chemistry problems, and the feature extraction technique can be easily generalized to other molecules.

Pheromones play an important role in coordinating male and female reproductive behavior in many aquatic species. Strikingly, however, relatively few peptide or protein pheromones have been characterized in invertebrates or vertebrates, as they are often difficult to isolate in sufficient quantities to characterize either biochemically or behaviorally. Furthermore, as we have seen with the marine mollusk Aplysia, protein pheromones require interaction with other proteins to exert their effect. Egg laying in Aplysia, which can be reliably induced by injection of egg-laying hormone, induces secretion of at least four distinct proteins from Temptin, a component of the complex of water-borne protein pheromones that stimulate attraction and mating behavior in the marine mollusk Aplysia, has sequence homology to the epidermal growth factor (EGF)-like domains of higher organisms that mediate protein-cell surface contact during fertilization and blood coagulation. In this work, recombinant temptin for structural and functional studies was produced in Escherichia coli using a cold shock promoter and purified by RP-HPLC. CD spectra confirmed a predominantly b-sheet structure. Two disulfide bonds were determined via limited proteolysis and MS. One internal disulfide (Cys57-Cys77) was predicted from initial alignments with class I EGF-like domains; the second, between Cys18 and Cys103, could protect temptin against proteolysis in seawater and stabilize its interacting surface. A three-dimensional model of temptin was prepared with our MPACK suite, based on the Ca 2+ -binding, EGF-like domain of the extracellular matrix protein fibrillin. Two temptin residues, Trp52 and Trp79, which align with cysteine residues conserved in fibrillins, lie adjacent to and could stabilize the disulfide bonds and a proposed metal-binding loop. The water-borne pheromone attractin in egg cordon eluates is complexed with other proteins. Docking results with our model and the NMR structure of attractin suggest that one face of temptin interacts with the pheromone, perhaps controlling its access to the cellular receptors. Gel shifts confirmed that temptin complexes with wild-type attractin. These results indicate that temptin, analogous to the role of fibrillin in controlling transforming growth factor-b concentration, modulates pheromone signaling by direct binding to attractin.Abbreviations EGF, epidermal growth factor; HFBA, heptafluorobutyric acid; IPTG, isopropyl thio-b-D-galactoside; TEE, translation-enhancing element; TGF-b, transforming growth factor-b.

Flavitrack contains over 590 complete flavivirus genome/protein sequences and information on known mutations and literature references. Each sequence has been manually annotated according to its date and place of isolation, phenotype and lethality. Internal tools are provided to rapidly determine relationships between viruses in Flavitrack and sequences provided by the user.