PocketMatch: A new algorithm to compare binding sites in protein structures

Recognizing similarities and deriving relationships among protein molecules is a fundamental requirement in present-day biology. Similarities can be present at various levels which can be detected through comparison of protein sequences or their structural folds. In some cases similarities obscure at these levels could be present merely in the substructures at their binding sites. Inferring functional similarities between protein molecules by comparing their binding sites is still largely exploratory and not as yet a routine protocol. One of the main reasons for this is the limitation in the choice of appropriate analytical tools that can compare binding sites with high sensitivity. To benefit from the enormous amount of structural data that is being rapidly accumulated, it is essential to have high throughput tools that enable large scale binding site comparison. Results: Here we present a new algorithm PocketMatch for comparison of binding sites in a frame invariant manner. Each binding site is represented by 90 lists of sorted distances capturing shape and chemical nature of the site. The sorted arrays are then aligned using an incremental alignment method and scored to obtain PMScores for pairs of sites. A comprehensive sensitivity analysis and an extensive validation of the algorithm have been carried out. A comparison with other site matching algorithms is also presented. Perturbation studies where the geometry of a given site was retained but the residue types were changed randomly, indicated that chance similarities were virtually non-existent. Our analysis also demonstrates that shape information alone is insufficient to discriminate between diverse binding sites, unless combined with chemical nature of amino acids. Conclusion: A new algorithm has been developed to compare binding sites in accurate, efficient and high-throughput manner. Though the representation used is conceptually simplistic, we demonstrate that along with the new alignment strategy used, it is sufficient to enable binding comparison with high sensitivity. Novel methodology has also been presented for validating the algorithm for accuracy and sensitivity with respect to geometry and chemical nature of the site. The method is also fast and takes about 1/250(th) second for one comparison on a single processor. A parallel version on BlueGene has also been implemented.

Perception of ultraviolet light by crab spiders and its role in selection of hunting sites

The perception of ultraviolet (UV) light by spiders has so far been only demonstrated in salticids. Crab spiders (Thomisidae) hunt mostly on flowers and need to find appropriate hunting sites. Previous studies have shown that some crab spiders that reflect UV light use UV contrast to enhance prey capture. The high UV contrast can be obtained either by modulation of body colouration or active selection of appropriate backgrounds for foraging. We show that crab spiders (Thomisus sp.)hunting on Spathiphyllum plants use chromatic contrast, especially UV contrast, to make themselves attractive to hymenopteran prey. Apart from that, they are able to achieve high UV contrast by active selection of non-UV reflecting surfaces when given a choice of UV-reflecting and non-UV reflecting surfaces in the absence of odour cues. Honeybees (Apis cerana) approached Spathiphyllum plants bearing crab spiders on which the spiders were high UV-contrast targets with greater frequency than those plants on which the UV-contrast of the spiders was low. Thus, crab spiders can perceive UV and may use it to choose appropriate backgrounds to enhance prey capture, by exploiting the attraction of prey such as honeybees to UV.

Cholinergic Muscarinic Receptors in Human Fetal Brain: Ontogeny of [3H]Quinuclidinyl Benzilate Binding Sites in Corpus Striatum, Brainstem, and Cerebellum

The ontogeny of muscarinic receptors was studied in human fetal striatum, brainstem, and cerebellum to investigate general principles of synaptogenesis as well as the physiological balance between various chemical synapses during development in a given region of the brain. [3H]Quinuclidinyl benzilate ([-'H]QNB) binding was assayed in total particulate fraction (TPF) from various parts of brain. In the corpus striatum, QNB binding sites are present at 16 weeks of gestation (average concentration 180 fmol/mg protein of TPF), slowly increase up to 24 weeks (average concentration 217 fmol/mg protein), and rapidly increase during the third trimester to 480 fmol/mg protein of TPF. In contrast, dopaminergic receptors exist as two subpopulations. one with low affinity and the other with high affinity up to the 24th week of gestation; all of them acquire the highaffinity characteristic during the third trimester. In brainstem, the muscarinic receptors show maximum concentration by 16 weeks of age (360 fmolimg protein of TPF). Subsequently the muscarinic receptor concentration shows a gradual decline in the brainstem. In cerebellum, except for a slight increase at 24 weeks (average concentration 90 fmol/mg protein of TPF), the receptor concentration remained nearly constant at about 60-70 fmolimg protein of TPF throughout fetal life. This study demonstrates that the ontogeny of muscarinic receptors varies among the different regions, and the patterns observed suggest that receptor formation occurs principally in the third trimester. Also noteworthy is the finding that the QNB binding sites decreased in all regions of the human brain during adult life. Key Words: Cholinergic muscarinic receptors-Quinuclidinyl benzilate-Corpus striaturn-Brainstem-Cerebellum. Ravikumar B. V. and Sastry P. S. Cholinergic muscarinic receptors in human fetal brain: Ontogeny of [3H]quinuclidinyl benzilate binding sites in corpus striatum, brainstem, and cerebellum. J. Neurochem. 45, 1948- 1950 (1985).

Fluorescent labelling of strychnine: A novel approach for recognition of strychnine binding sites on neuronal membrane

trychnine was coupled to fluorescein isothiocyanate to mark strychnine binding sites in spinal cord of rat. Specific binding of strychnine could be demonstrated in synaptosomal fraction. Addition of glycine to the strychninised membrane led to a decrease in fluorescence indicating same receptor loci.

A measure of logical complexity of programs

The research in software science has so far been concentrated on three measures of program complexity: (a) software effort; (b) cyclomatic complexity; and (c) program knots. In this paper we propose a measure of the logical complexity of programs in terms of the variable dependency of sequence of computations, inductive effort in writing loops and complexity of data structures. The proposed complexity mensure is described with the aid of a graph which exhibits diagrammatically the dependence of a computation at a node upon the computation of other (earlier) nodes. Complexity measures of several example programs have been computed and the related issues have been discussed. The paper also describes the role played by data structures in deciding the program complexity.

Synergistic Execution of Stream Programs on Multicores with Accelerators

The StreamIt programming model has been proposed to exploit parallelism in streaming applications oil general purpose multicore architectures. The StreamIt graphs describe task, data and pipeline parallelism which can be exploited on accelerators such as Graphics Processing Units (GPUs) or CellBE which support abundant parallelism in hardware. In this paper, we describe a novel method to orchestrate the execution of if StreamIt program oil a multicore platform equipped with an accelerator. The proposed approach identifies, using profiling, the relative benefits of executing a task oil the superscalar CPU cores and the accelerator. We formulate the problem of partitioning the work between the CPU cores and the GPU, taking into account the latencies for data transfers and the required buffer layout transformations associated with the partitioning, as all integrated Integer Linear Program (ILP) which can then be solved by an ILP solver. We also propose an efficient heuristic algorithm for the work-partitioning between the CPU and the GPU, which provides solutions which are within 9.05% of the optimal solution on an average across the benchmark Suite. The partitioned tasks are then software pipelined to execute oil the multiple CPU cores and the Streaming Multiprocessors (SMs) of the GPU. The software pipelining algorithm orchestrates the execution between CPU cores and the GPU by emitting the code for the CPU and the GPU, and the code for the required data transfers. Our experiments on a platform with 8 CPU cores and a GeForce 8800 GTS 512 GPU show a geometric mean speedup of 6.94X with it maximum of 51.96X over it single threaded CPU execution across the StreamIt benchmarks. This is a 18.9% improvement over it partitioning strategy that maps only the filters that cannot be executed oil the GPU - the filters with state that is persistent across firings - onto the CPU.

Probing into the role of conserved N-glycosylation sites in the Tyrosinase glycoprotein family

N-linked glycosylation has a profound effect on the proper folding, oligomerization and stability of glycoproteins. These glycans impart many properties to proteins that may be important for their proper functioning, besides having a tendency to exert a chaperone-like effect on them. Certain glycosylation sites in a protein however, are more important than other sites for their function and stability. It has been observed that some N-glycosylation sites are conserved over families of glycoproteins over evolution, one such being the tyrosinase related protein family. The role of these conserved N-glycosylation sites in their trafficking, sorting, stability and activity has been examined here. By scrutinizing the different glycosylation sites on this family of glycoproteins it was inferred that different sites in the same family of polypeptides can perform distinct functions and conserved sites across the paralogues may perform diverse functions.

Prediction of bubble growth rates and departure volumes in nucleate boiling at isolated sites

A model has been developed to predict heat transfer rates and sizes of bubbles generated during nucleate pool boiling. This model assumes conduction and a natural convective heat transfer mechanism through the liquid layer under the bubble and transient conduction from the bulk liquid. The temperature of the bulk liquid in the vicinity of the bubble is obtained by assuming a turbulent natural convection process from the hot plate to the liquid bulk. The shape of the bubble is obtained by equilibrium analysis. The bubble departure condition is predicted by a force balance equation. Good agreement has been found between the bubble radii predicted by the present theory and the ones obtained experimentally.

Studies on bicuculline binding sites on neuronal membrane using fluorescent antibody technique: Comparative binding of gaba and bicuculline

Highly purified fluorescent labelled anti-bicuculline antibodies were used to mark bicuculline binding sites in cerebral cortex of monkey brain. Specific binding of bicuculline could be demonstrated in the synaptosomal fraction, when bicuculline was added both Image and Image . Addition of γ-aminobutyric acid (GABA) to the bicucullinised membrane led to a decrease in fluorescence indicating same receptor loci and establishing GABA-bicuculline antagonism at a molecular level.

Inhibition of Human Two-Pore Domain K+ Channel TREK1 by Local Anesthetic Lidocaine: Negative Cooperativity and Half-of-Sites Saturation Kinetics

TWIK-related K+ channel TREK1, a background leak K+ channel, has been strongly implicated as the target of several general and local anesthetics. Here, using the whole-cell and single-channel patch-clamp technique, we investigated the effect of lidocaine, a local anesthetic, on the human (h) TREK1 channel heterologously expressed in human embryonic kidney 293 cells by an adenoviral-mediated expression system. Lidocaine, at clinical concentrations, produced reversible, concentration-dependent inhibition of hTREK1 current, with IC50 value of 180 mu M, by reducing the single-channel open probability and stabilizing the closed state. We have identified a strategically placed unique aromatic couplet (Tyr352 and Phe355) in the vicinity of the protein kinase A phosphorylation site, Ser348, in the C-terminal domain (CTD) of hTREK1, that is critical for the action of lidocaine. Furthermore, the phosphorylation state of Ser348 was found to have a regulatory role in lidocaine-mediated inhibition of hTREK1. It is interesting that we observed strong intersubunit negative cooperativity (Hill coefficient = 0.49) and half-of-sites saturation binding stoichiometry (half-reaction order) for the binding of lidocaine to hTREK1. Studies with the heterodimer of wild-type (wt)-hTREK1 and Delta 119 C-terminal deletion mutant (hTREK1(wt)-Delta 119) revealed that single CTD of hTREK1 was capable of mediating partial inhibition by lidocaine, but complete inhibition necessitates the cooperative interaction between both the CTDs upon binding of lidocaine. Based on our observations, we propose a model that explains the unique kinetics and provides a plausible paradigm for the inhibitory action of lidocaine on hTREK1.

Additional binding sites in lysozyme. X-ray analysis of lysozyme complexes with bromophenol red and bromophenol blue

The binding sites in hen egg-white lysozyme for neutral bromophenol red (BPR) and ionized bromophenol blue (BPB) have been characterized at 2 Å resolution. In either case, the dye-bound enzyme is active against the polysaccharide, but not against the cell wall. Both binding sites are outside, but close to, the hexasaccharide binding cleft in the enzyme. The binding site of BPR made up of Arg5, Lys33, Phe34, Asn37, Phe38, Ala122, Trp123 and possibly Arg125, is dose to subsite F while that of BPB made up of Tyr20, Arg21, Asn93, Lys96, Lys97 and Ser100, is close to subsites A and B. The binding sites of the neutral dye and the ionized dye are thus spatially far apart. The peptide component of the bacterial cell wall probably interacts with these cells during enzyme action. Such interactions are perhaps necessary for appropriately positioning the enzyme molecule on the bacterial cell wall.

Structural Modes of Stabilization of Permissive Phosphorylation Sites in Protein Kinases: Distinct Strategies in Ser/Thr and Tyr Kinases

Protein kinases phosphorylate several cellular proteins providing control mechanisms for various signalling processes. Their activity is impeded in a number of ways and restored by alteration in their structural properties leading to a catalytically active state. Most protein kinases are subjected to positive and negative regulation by phosphorylation of Ser/Thr/Tyr residues at specific sites within and outside the catalytic core. The current review describes the analysis on 3D structures of protein kinases that revealed features distinct to active states of Ser/Thr and Tyr kinases. The nature and extent of interactions among well-conserved residues surrounding the permissive phosphorylation sites differ among the two classes of enzymes. The network of interactions of highly conserved Arg preceding the catalytic base that mediates stabilization of the activation segment exemplifies such diverse interactions in the two groups of kinases. The N-terminal and the C-terminal lobes of various groups of protein kinases further show variations in their extent of coupling as suggested from the extent of interactions between key functional residues in activation segment and the N-terminal αC-helix. We observe higher similarity in the conformations of ATP bound to active forms of protein kinases compared to ATP conformations in the inactive forms of kinases. The extent of structural variations accompanying phosphorylation of protein kinases is widely varied. The comparison of their crystal structures and the distinct features observed are hoped to aid in the understanding of mechanisms underlying the control of the catalytic activity of distinct subgroups of protein kinases.

Solution NMR studies of acetohydroxy acid synthase I: Identification of the sites of inter-subunit interactions using multidimensional NMR methods

The novel multidomain organization in the multimeric Escherichia coli AHAS I (ilvBN) enzyme has been dissected to generate polypeptide fragments. These fragments when cloned, expressed and purified reassemble in the presence of cofactors to yield a catalytically competent enzyme. Structural characterization of AHAS has been impeded due to the fact that the holoenzyme is prone to dissociation leading to heterogeneity in samples. Our approach has enabled the structural characterization using high-resolution nuclear magnetic resonance methods. Near complete sequence specific NMR assignments for backbone H-N, N-15, C-13 alpha and C-13(beta) atoms of the FAD binding domain of ilvB have been obtained on samples isotopically enriched in H-2, C-13 and N-15. The secondary structure determined on the basis of observed C-13(alpha) secondary chemical shifts and sequential NOEs indicates that the secondary structure of the FAD binding domain of E. coli AHAS large Subunit (ilvB) is similar to the structure of this domain in the catalytic subunit of yeast AHAS. Protein-protein interactions involving the regulatory subunit (ilvN) and the domains of the catalytic subunit (ilvB) were studied using circular dichroic and isotope edited solution nuclear magnetic resonance spectroscopic methods. Observed changes in circular dichroic spectra indicate that the regulatory subunit (ilvN) interacts with ilvB alpha and ilvB beta domains of the catalytic subunit and not with the ilvB gamma domain. NMR chemical shift mapping methods show that ilvN binds close to the FAD binding site in ilvB beta and proximal to the intrasubunit ilvB alpha/ilvB beta domain interface. The implication of this interaction on the role of the regulatory subunit oil the activity of the holoenzyme is discussed. NMR studies of the regulatory domains show that these domains are structured in solution. Preliminary evidence for the interaction of ilvN with the metabolic end product of the pathway, viz., valine is also presented.

The RecA intein of Mycobacterium tuberculosis promotes cleavage of ectopic DNA sites. Implications for the dispersal of inteins in natural populations

The RecA intein of Mycobacterium tuberculosis, a novel double-stranded DNA endonuclease, requires both Mn(2+) and ATP for efficient cleavage of the inteinless recA allele. In this study, we show that Mg(2+) alone was sufficient to stimulate PI-MtuI to cleave double-stranded DNA at ectopic sites. In the absence of Mg(2+), PI-MtuI formed complexes with topologically different forms of DNA containing ectopic recognition sequences with equal affinity but failed to cleave DNA. We observed that PI-MtuI was able to inflict double-strand breaks robustly within the ectopic recognition sequence to generate either a blunt end or 1-2-nucleotide 3'-hydroxyl overhangs. Mutational analyses of the presumptive metal ion-binding ligands (Asp(122), Asp(222), and Glu(220)) together with immunoprecipitation assays provided compelling evidence to link both the Mg(2+)- and Mn(2+) and ATP-dependent endonuclease activities to PI-MtuI. The kinetic mechanism of PI-MtuI promoted cleavage of ectopic DNA sites proceeded through a sequential mechanism with transient accumulation of nicked circular duplex DNA as an intermediate. Together, these data suggest that PI-MtuI, like group II introns, might mediate ectopic DNA transposition and hence its lateral transfer in natural populations.

Software Pipelined Execution of Stream Programs on GPUs

The StreamIt programming model has been proposed to exploit parallelism in streaming applications on general purpose multi-core architectures. This model allows programmers to specify the structure of a program as a set of filters that act upon data, and a set of communication channels between them. The StreamIt graphs describe task, data and pipeline parallelism which can be exploited on modern Graphics Processing Units (GPUs), as they support abundant parallelism in hardware. In this paper, we describe the challenges in mapping StreamIt to GPUs and propose an efficient technique to software pipeline the execution of stream programs on GPUs. We formulate this problem - both scheduling and assignment of filters to processors - as an efficient Integer Linear Program (ILP), which is then solved using ILP solvers. We also describe a novel buffer layout technique for GPUs which facilitates exploiting the high memory bandwidth available in GPUs. The proposed scheduling utilizes both the scalar units in GPU, to exploit data parallelism, and multiprocessors, to exploit task and pipelin parallelism. Further it takes into consideration the synchronization and bandwidth limitations of GPUs, and yields speedups between 1.87X and 36.83X over a single threaded CPU.