Comparison Of Binding Sites Research Articles

Comparative Study of Allosteric GPCR Binding Sites and Their Ligandability Potential.

The steadily growing number of experimental G-protein-coupled receptor (GPCR) structures has revealed diverse locations of allosteric modulation, and yet few drugs target them. This gap highlights the need for a deeper understanding of allosteric modulation in GPCR drug discovery. The current work introduces a systematic annotation scheme to structurally classify GPCR binding sites based on receptor class, transmembrane helix contacts, and, for membrane-facing sites, membrane sublocation. This GPCR specific annotation scheme was applied to 107 GPCR structures bound by small molecules contributing to 24 distinct allosteric binding sites for comparative evaluation of three binding site detection methods (BioGPS, SiteMap, and FTMap). BioGPS identified the most in 22 of 24 sites. In addition, our property analysis showed that extrahelical allosteric ligands and binding sites represent a distinct chemical space characterized by shallow pockets with low volume, and the corresponding allosteric ligands showed an enrichment of halogens. Furthermore, we demonstrated that combining receptor and ligand similarity can be a viable method for ligandability assessment. One challenge regarding site prediction is the ligand shaping effect on the observed binding site, especially for extrahelical sites where the ligand-induced effect was most pronounced. To our knowledge, this is the first study presenting a binding site annotation scheme standardized for GPCRs, and it allows a comparison of allosteric binding sites across different receptors in an objective way. The insight from this study provides a framework for future GPCR binding site studies and highlights the potential of targeting allosteric sites for drug development.

Structural Comparison of Substrate Binding Sites in Dehaloperoxidase A and B.

Dehalperoxidase (DHP) has diverse catalytic activities depending on the substrate binding conformation, pH, and dynamics in the distal pocket above the heme. According to our hypothesis, the molecular structure of the substrate and binding orientation in DHP guide enzymatic function. Enzyme kinetic studies have shown that the catalytic activity of DHP B is significantly higher than that of DHP A despite 96% sequence homology. There are more than 30 substrate-bound structures with DHP B, each providing insight into the nature of enzymatic binding at the active site. By contrast, the only X-ray crystallographic structures of small molecules in a complex with DHP A are phenols. This study is focused on investigating substrate binding in DHP A to compare with DHP B structures. Fifteen substrates were selected that were known to bind to DHP B in the crystal to test whether soaking substrates into DHP A would yield similar structures. Five of these substrates yielded X-ray crystal structures of substrate-bound DHP A, namely, 2,4-dichlorophenol (1.48 Å, PDB: 8EJN), 2,4-dibromophenol (1.52 Å, PDB: 8VSK), 4-nitrophenol (2.03 Å, PDB: 8VKC), 4-nitrocatechol (1.40 Å, PDB: 8VKD), and 4-bromo-o-cresol (1.64 Å, PDB: 8VZR). For the remaining substrates that bind to DHP B, such as cresols, 5-bromoindole, benzimidazole, 4,4-biphenol, 4.4-ethylidenebisphenol, 2,4-dimethoxyphenol, and guaiacol, the electron density maps in DHP A are not sufficient to determine the presence of the substrates, much less their orientation. In our hands, only phenols, 4-Br-o-cresol, and 4-nitrocatechol can be soaked into crystalline DHP A. None of the larger substrates were observed to bind. A minimum of seven hanging drops were selected for soaking with more than 50 crystals screened for each substrate. The five high-quality examples of direct comparison of modes of binding in DHP A and B for the same substrate provide further support for the hypothesis that the substrate-binding conformation determines the enzyme function of DHP.

Deciphering the Genetic Variation: A Comparative Analysis of Parental and Attenuated Strains of the QXL87 Vaccine for Infectious Bronchitis.

The QXL87 live attenuated vaccine strain for infectious bronchitis represents the first approved QX type (GI-19 lineage) vaccine in China. This strain was derived from the parental strain CK/CH/JS/2010/12 through continuous passage in SPF chicken embryos. To elucidate the molecular mechanism behind its attenuation, whole-genome sequencing was conducted on both the parental and attenuated strains. Analysis revealed 145 nucleotide mutations in the attenuated strain, leading to 48 amino acid mutations in various proteins, including Nsp2 (26), Nsp3 (14), Nsp4 (1), S (4), 3a (1), E (1), and N (1). Additionally, a frameshift mutation caused by a single base insertion in the ORFX resulted in a six-amino-acid extension. Subsequent comparison of post-translational modification sites, protein structure, and protein-protein binding sites between the parental and attenuated strains identified three potential virulence genes: Nsp2, Nsp3, and S. The amino acid mutations in these proteins not only altered their conformation but also affected the distribution of post-translational modification sites and protein-protein interaction sites. Furthermore, three potential functional mutation sites-P106S, A352T, and L472F, all located in the Nsp2 protein-were identified through PROVEAN, PolyPhen, and I-Mutant. Overall, our findings suggest that Nsp2, Nsp3, and S proteins may play a role in modulating IBV pathogenicity, with a particular focus on the significance of the Nsp2 protein. This study contributes to our understanding of the molecular mechanisms underlying IBV attenuation and holds promise for the development of safer live attenuated IBV vaccines using reverse genetic approaches.

An Intriguing Purview on the Design of Macrocyclic Inhibitors for Unexplored Protein Kinases through Their Binding Site Comparison.

Kinases play an important role in regulating various intracellular signaling pathways that control cell proliferation, differentiation, survival, and other cellular processes, and their deregulation causes more than 400 diseases. Consequently, macrocyclization can be considered a noteworthy approach to developing new therapeutic agents for human diseases. Macrocyclization has emerged as an effective drug discovery strategy over the past decade to improve target selectivity and potency of small molecules. Small compounds with linear structures upon macrocyclization can lead to changes in their physicochemical and biological properties by firmly reducing conformational flexibility. A number of distinct protein kinases exhibit similar binding sites. Comparison of protein binding sites provides crucial insights for drug discovery and development. Binding site similarities are helpful in understanding polypharmacology, identifying potential off-targets, and repurposing known drugs. In this review, we focused on comparing the binding sites of those kinases for which macrocyclic inhibitors are available/studied so far. Furthermore, we calculated the volume of the binding site pocket for each targeted kinase and then compared it with the binding site pocket of the kinase for which only acyclic inhibitors were designed to date. Our review and analysis of several explored kinases might be useful in targeting new protein kinases for macrocyclic drug discovery.

Estimating the Similarity between Protein Pockets.

With the exponential increase in publicly available protein structures, the comparison of protein binding sites naturally emerged as a scientific topic to explain observations or generate hypotheses for ligand design, notably to predict ligand selectivity for on- and off-targets, explain polypharmacology, and design target-focused libraries. The current review summarizes the state-of-the-art computational methods applied to pocket detection and comparison as well as structural druggability estimates. The major strengths and weaknesses of current pocket descriptors, alignment methods, and similarity search algorithms are presented. Lastly, an exhaustive survey of both retrospective and prospective applications in diverse medicinal chemistry scenarios illustrates the capability of the existing methods and the hurdle that still needs to be overcome for more accurate predictions.

Methylation pattern of polymorphically imprinted nc886 is not conserved across mammalia.

BackgroundIn humans, the nc886 locus is a polymorphically imprinted metastable epiallele. Periconceptional conditions have an effect on the methylation status of nc886, and further, this methylation status is associated with health outcomes in later life, in line with the Developmental Origins of Health and Disease (DOHaD) hypothesis. Animal models would offer opportunities to study the associations between periconceptional conditions, nc886 methylation status and metabolic phenotypes further. Thus, we set out to investigate the methylation pattern of the nc886 locus in non-human mammals.DataWe obtained DNA methylation data from the data repository GEO for mammals, whose nc886 gene included all three major parts of nc886 and had sequency similarity of over 80% with the human nc886. Our final sample set consisted of DNA methylation data from humans, chimpanzees, bonobos, gorillas, orangutangs, baboons, macaques, vervets, marmosets and guinea pigs.ResultsIn human data sets the methylation pattern of nc886 locus followed the expected bimodal distribution, indicative of polymorphic imprinting. In great apes, we identified a unimodal DNA methylation pattern with 50% methylation level in all individuals and in all subspecies. In Old World monkeys, the between individual variation was greater and methylation on average was close to 60%. In guinea pigs the region around the nc886 homologue was non-methylated. Results obtained from the sequence comparison of the CTCF binding sites flanking the nc886 gene support the results on the DNA methylation data.ConclusionsOur results indicate that unlike in humans, nc886 is not a polymorphically imprinted metastable epiallele in non-human primates or in guinea pigs, thus implying that animal models are not applicable for nc886 research. The obtained data suggests that the nc886 region may be classically imprinted in great apes, and potentially also in Old World monkeys, but not in guinea pigs.

In Silico Analysis and Comparison of Pathogenesis-related Transcription Factor Binding Sites in Arabidopsis DNA Repair Gene Promoters

In Silico Analysis and Comparison of Pathogenesis-related Transcription Factor Binding Sites in Arabidopsis DNA Repair Gene Promoters

Analysis of role of aromatic residues in extracellular loop 2 of Prokineticin receptor 2 in ligand binding probed with genetically encoded photo-crosslinkers

Prokineticin 2 (PK2) and Prokineticin 2 beta (PK2β), products of alternative splicing of pk2 gene, are chemokine-like proteins. While PK2 mediates its biological activities by signaling with the same efficiency through two homologous G protein coupled receptors, prokineticin receptor 1 (PKR1) and prokineticin receptor 2 (PKR2), PK2β is able to bind specifically PKR1.Extracellular loop 2 (ECL2) of chemokine receptors is a part of a transmembrane (TM) ligand binding site. In the ECL2 of PKR2 is present, as well as in all chemokine receptors, an aromatic residue cluster, involving tryptophan 212 localized four residues after an ECL2 conserved cysteine, and Phenylalanine 198 located in the top of TM 4.In this work, the photoactivatable unnatural amino acid p-benzoyl-L-phenylalanine is incorporated by amber codon suppression technology into PKR2 in position 212. Experiments of photoactivatable cross-linking demonstrated the role of tryptophan in position 212 for binding the ligand contacting Tryptophan in position 24. We also analyzed the role of Phenylalanine 198 in the specificity of PKRs binding. The comparison of TM-bundle binding sites between PKR1 and PKR2 revealed that they are completely conserved except for one residue: valine 207 in human PKR1, which is phenylalanine 198 in human PKR2. The F198V mutation in PKR2 permits to obtain a receptor able to bind more efficiently PK2β, a ligand highly specific for PKR1.

The genome analysis of halotolerant Sphingobium yanoikuyae YC-XJ2 with aryl organophosphorus flame retardants degrading capacity and characteristics of related phosphotriesterase

The genome of Sphingobium yanoikuyae YC-XJ2 with aryl organophosphorus flame retardants (aryl-OPFRs) degrading activity was sequenced and analyzed systematically. The gene sy-pte was cloned and expressed, and the molecular mechanism of aryl-OPFRs degradation was elucidated. Sy-PTE exhibited enzymatic properties of alkali, high temperature and salt resistance. The hydrolyzation of EHDPP by phosphotriesterase was first reported. The specific activity of Sy-PTE was 67.0 ± 1.4 (U mg−1) for TPhP with kcat/Km and Vmax values of 4.8 ± 0.6 × 103 (mM−1·s−1) and 0.29 ± 0.02 (μM·s−1), Km of 15.6 ± 0.34 (μM), kcat of 75.02 ± 1.5 s−1, and need no metal ion activator to assist. The structural comparison of the metal binding sites with Sb-PTE (PDB ID 5HRM) showed that the transition from uncharged Met-259 to positively charged Asn-259 in Sy-PET might affect the binding of Mn2+ to His-258 because of the same electrical properties, and the transition from Gly-473 without side chains in Sb-PTE to Ala-473 might cause steric hindrance in binding of Mn2+ to His-475. Based on enzymatic advantages, Sy-PTE should be a promising engineered enzyme in the future for bioremediation of aryl-OPFRs polluted sites.

DeeplyTough: Learning Structural Comparison of Protein Binding Sites.

Protein pocket matching, or binding site comparison, is of importance in drug discovery. Identification of similar binding pockets can help guide efforts for hit-finding, understanding polypharmacology, and characterization of protein function. The design of pocket matching methods has traditionally involved much intuition and has employed a broad variety of algorithms and representations of the input protein structures. We regard the high heterogeneity of past work and the recent availability of large-scale benchmarks as an indicator that a data-driven approach may provide a new perspective. We propose DeeplyTough, a convolutional neural network that encodes a three-dimensional representation of protein pockets into descriptor vectors that may be compared efficiently in an alignment-free manner by computing pairwise Euclidean distances. The network is trained with supervision (i) to provide similar pockets with similar descriptors, (ii) to separate the descriptors of dissimilar pockets by a minimum margin, and (iii) to achieve robustness to nuisance variations. We evaluate our method using three large-scale benchmark datasets, on which it demonstrates excellent performance for held-out data coming from the training distribution and competitive performance when the trained network is required to generalize to datasets constructed independently. DeeplyTough is available at https://github.com/BenevolentAI/DeeplyTough.

Ehrlichia chaffeensis TRP47 enters the nucleus via a MYND-binding domain-dependent mechanism and predominantly binds enhancers of host genes associated with signal transduction, cytoskeletal organization, and immune response.

Ehrlichia chaffeensis is an obligately intracellular bacterium that establishes infection in mononuclear phagocytes through largely undefined reprogramming strategies including modulation of host gene transcription. In this study, we demonstrate that the E. chaffeensis effector TRP47 enters the host cell nucleus and binds regulatory regions of host genes relevant to infection. TRP47 was observed in the nucleus of E. chaffeensis-infected host cells, and nuclear localization was dependent on a variant MYND-binding domain. An electrophoretic mobility shift assay (EMSA) demonstrated that TRP47 directly binds host DNA via its tandem repeat domain. Utilizing chromatin immunoprecipitation followed by high-throughput DNA sequencing (ChIP-seq) with E. chaffeensis-infected cells, TRP47 was found to bind at multiple sites in the human genome (n = 2,051 at p < 10−30). Ontology analysis identified genes involved in functions such as immune response, cytoskeletal organization, and signal transduction. TRP47-bound genes included RNA-coding genes, many of these linked to cell proliferation or apoptosis. Comparison of TRP47 binding sites with those of previously-identified E. chaffeensis nucleomodulins identified multiple genes and gene functional categories in common including intracellular transport, cell signaling, and transcriptional regulation. Further, motif analysis followed by EMSA with synthetic oligonucleotides containing discovered motifs revealed a conserved TRP47 DNA-binding motif. This study reveals that TRP47 is a nucleomodulin that enters the nucleus via a MYND-binding domain and appears to play a role in host cell reprogramming by regulation of transcription.

Interspecies comparison of putative ligand binding sites of human, rat and mouse P-glycoprotein

Prior to the clinical phases of testing, safety, efficacy and pharmacokinetic profiles of lead compounds are evaluated in animal studies. These tests are primarily performed in rodents, such as mouse and rats. In order to reduce the number of animal experiments, computational models that predict the outcome of these studies and thus aid in prioritization of preclinical candidates are heavily needed. However, although computational models for human off-target interactions with decent quality are available, they cannot easily be transferred to rodents due to lack of respective data. In this study, we assess the transferability of human P-glycoprotein activity data for development of in silico models to predict in vivo effects in rats and mouse using a structure-based approach. P-glycoprotein (P-gp) is an ATP-dependent efflux transporter that transports xenobiotic compounds such as toxins and drugs out of cells and has a broad substrate and inhibitor specificity. Being mostly expressed at barriers, it influences the bioavailability of drugs and thus contributes also to toxicity. Comparison of the binding site interaction profiles of human, rat and mouse P-gp derived from docking studies with a set of common inhibitors suggests that the inhibitors share potentially similar binding modes. These findings encourage the use of in vitro human P-gp data for predicting in vivo effects in rodents and thus contributes to the 3Rs (Replace, Reduce and Refine) of animal experiments.

Superposition-free comparison and clustering of antibody binding sites: implications for the prediction of the nature of their antigen

We describe here a superposition free method for comparing the surfaces of antibody binding sites based on the Zernike moments and show that they can be used to quickly compare and cluster sets of antibodies. The clusters provide information about the nature of the bound antigen that, when combined with a method for predicting the number of direct antibody antigen contacts, allows the discrimination between protein and non-protein binding antibodies with an accuracy of 76%. This is of relevance in several aspects of antibody science, for example to select the framework to be used for a combinatorial antibody library.

The Production and Utilization of GDP-glucose in the Biosynthesis of Trehalose 6-Phosphate by Streptomyces venezuelae

Trehalose-6-phosphate synthase OtsA from streptomycetes is unusual in that it uses GDP-glucose as the donor substrate rather than the more commonly used UDP-glucose. We now confirm that OtsA from Streptomyces venezuelae has such a preference for GDP-glucose and can utilize ADP-glucose to some extent too. A crystal structure of the enzyme shows that it shares twin Rossmann-like domains with the UDP-glucose-specific OtsA from Escherichia coli. However, it is structurally more similar to Streptomyces hygroscopicus VldE, a GDP-valienol-dependent pseudoglycosyltransferase enzyme. Comparison of the donor binding sites reveals that the amino acids associated with the binding of diphosphoribose are almost all identical in these three enzymes. By contrast, the amino acids associated with binding guanine in VldE (Asn, Thr, and Val) are similar in S. venezuelae OtsA (Asp, Ser, and Phe, respectively) but not conserved in E. coli OtsA (His, Leu, and Asp, respectively), providing a rationale for the purine base specificity of S. venezuelae OtsA. To establish which donor is used in vivo, we generated an otsA null mutant in S. venezuelae. The mutant had a cell density-dependent growth phenotype and accumulated galactose 1-phosphate, glucose 1-phosphate, and GDP-glucose when grown on galactose. To determine how the GDP-glucose is generated, we characterized three candidate GDP-glucose pyrophosphorylases. SVEN_3027 is a UDP-glucose pyrophosphorylase, SVEN_3972 is an unusual ITP-mannose pyrophosphorylase, and SVEN_2781 is a pyrophosphorylase that is capable of generating GDP-glucose as well as GDP-mannose. We have therefore established how S. venezuelae can make and utilize GDP-glucose in the biosynthesis of trehalose 6-phosphate.

Core binding factor (CBF) is required for Epstein-Barr virus EBNA3 proteins to regulate target gene expression.

ChIP-seq performed on lymphoblastoid cell lines (LCLs), expressing epitope-tagged EBNA3A, EBNA3B or EBNA3C from EBV-recombinants, revealed important principles of EBNA3 binding to chromatin. When combined with global chromatin looping data, EBNA3-bound loci were found to have a singular character, each directly associating with either EBNA3-repressed or EBNA3-activated genes, but not with both. EBNA3A and EBNA3C showed significant association with repressed and activated genes. Significant direct association for EBNA3B loci could only be shown with EBNA3B-repressed genes. A comparison of EBNA3 binding sites with known transcription factor binding sites in LCL GM12878 revealed substantial co-localization of EBNA3s with RUNX3—a protein induced by EBV during B cell transformation. The beta-subunit of core binding factor (CBFβ), that heterodimerizes with RUNX3, could co-immunoprecipitate robustly EBNA3B and EBNA3C, but only weakly EBNA3A. Depletion of either RUNX3 or CBFβ with lentivirus-delivered shRNA impaired epitope-tagged EBNA3B and EBNA3C binding at multiple regulated gene loci, indicating a requirement for CBF heterodimers in EBNA3 recruitment during target-gene regulation. ShRNA-mediated depletion of CBFβ in an EBNA3C-conditional LCL confirmed the role of CBF in the regulation of EBNA3C-induced and -repressed genes. These results reveal an important role for RUNX3/CBF during B cell transformation and EBV latency that was hitherto unexplored.

Comparative electrostatic analysis of adenylyl cyclase for isoform dependent regulation properties

The enzyme adenylyl cyclase (AC) plays a pivotal role in a variety of signal transduction pathways inside the cell, where it catalyzes the cyclization of adenosine triphosphate (ATP) into the second-messenger cyclic adenosine monophosphate (cAMP). Among other roles, AC regulates processes involved in neural plasticity, innervation of smooth muscles of the heart and the endocrine system of the pancreas. The functional diversity of AC is manifested in its different isoforms, each having a specific regulation pattern. There is an increasing amount of data available concerning the regulatory properties of AC isoforms, however little is known about the interactions on a structural level. Here, we conducted a comparative electrostatic analysis of the catalytic domains of all nine transmembrane AC isoforms with the aim of detecting, verifying and predicting the binding sites of molecular regulators on AC. The results provide support for the positioning of the binding site of the inhibitory protein Gi α at a pseudo-symmetric position to the stimulatory Gs α binding site. They also provide a structural interpretation of the Gβγ interaction with ACs 2, 4, and 7 and suggest a new binding site for RGS2. Comparison of the small molecule binding sites on AC shows that overall they have high electrostatic similarity, but regions of electrostatic differences are identified. These could provide a basis for the development of novel compounds with isoform-specific modulatory effects on AC. Proteins 2016; 84:1844-1858. © 2016 Wiley Periodicals, Inc.

CavSimBase: A Database for Large Scale Comparison of Protein Binding Sites

CavBase is a database containing information about the three-dimensional geometry and the physicochemical properties of putative protein binding sites. Analyzing CavBase data typically involves computing the similarity of pairs of binding sites. In contrast to sequence alignment, however, a structural comparison of protein binding sites is a computationally challenging problem, making large scale studies difficult or even infeasible. One possibility to overcome this obstacle is to precompute pairwise similarities in an all-against-all comparison, and to make these similarities subsequently accessible to data analysis methods. Pairwise similarities, once being computed, can also be used to equip CavBase with a neighborhood structure. Taking advantage of this structure, methods for problems such as similarity retrieval can be implemented efficiently. In this paper, we tackle the problem of performing an all-against-all comparison using CavBase, consisting of more than 200,000 protein cavities, by means of parallel computation and cloud computing techniques. We present the conceptual design and technical realization of a large-scale study to create a similarity database called CavSimBase. We illustrate how CavSimBase is constructed, is accessed, and is used to answer biological questions by data analysis and similarity retrieval.

Scientists Explore a Master Genetic Switch That Plays a Key Role in Energy Metabolism and Human Brain Evolution and Function

Scientists have long used comparative animal studies to better understand the nuances of human evolution, from making diverse body plans to the emergence of entirely powerful and unique features structures, including the human brain. Now, in a new study appearing in the advanced online edition of Molecular Biology and Evolution, corresponding authors Perdomo-Sabogal et al. (2016) have explored the global gene regulator GABPa to better understand its influence as a master switch. Their work honed in on human-specific DNA regions potentially contributing to human evolution in key functions including cell energy, division and death, and diseases ranging from brain disorders to diabetes. Many of the GABPa gene regulatory sites were found to be significantly enriched at genes important for brain, central nervous system and spinal cord functions that are unique to humans, an area ripe for future research, and GABPa deficits linked to Alzheimer’s and Parkinson diseases, breast cancer, and autism. GABPa exerts its profound genomic influence by a tandem repeat of four highly conserved DNA letters. The GGAA consensus motif is found to control genes in all primates, mice, dogs, and cows, underscoring its importance in these animals. By performing ChIP-Seq experiments in a human cell line, the group found that GABPa exercises its influence across the genome at 11,619 putative GABPa binding sites representing nearly 4,000 genes. Through a tour de force sequence comparison of the human GABPa binding sites with similar sequences from 34 mammals, they identified 224 GABPa binding sites unique to humans. Comparing chimpanzee, macaque, and human cell lines and mutating the GGAA sequence to more ancestral states to determine how these changes affected GABPa function, they further narrowed down their function. For instance, when the human GABPa binding sites were introduced into primate cells, gene expression was increased up to five times in 17 out of 18 cases, indicating that a recent human mutation in GABPa resulted in higher gene expression and forces behind the evolution of novel, uniquely human functions. “Mutations that cause changes in the regulation of gene activities are one of the major factors in shaping species during evolution,” said corresponding author Robert Querfurth. “Our study demonstrates how, out of the millions of DNA regions, in which we differ from other apes, we can sift out those that, in response to a specific regulatory protein, cause gene activity changes in human cells. In this set, we find genes involved in brain and breast development and also in diseases like Alzheimer’s and Parkinson’s.” Reference

A web server for analysis, comparison and prediction of protein ligand binding sites.

BackgroundOne of the major challenges in the field of system biology is to understand the interaction between a wide range of proteins and ligands. In the past, methods have been developed for predicting binding sites in a protein for a limited number of ligands.ResultsIn order to address this problem, we developed a web server named ‘LPIcom’ to facilitate users in understanding protein-ligand interaction. Analysis, comparison and prediction modules are available in the “LPIcom’ server to predict protein-ligand interacting residues for 824 ligands. Each ligand must have at least 30 protein binding sites in PDB. Analysis module of the server can identify residues preferred in interaction and binding motif for a given ligand; for example residues glycine, lysine and arginine are preferred in ATP binding sites. Comparison module of the server allows comparing protein-binding sites of multiple ligands to understand the similarity between ligands based on their binding site. This module indicates that ATP, ADP and GTP ligands are in the same cluster and thus their binding sites or interacting residues exhibit a high level of similarity. Propensity-based prediction module has been developed for predicting ligand-interacting residues in a protein for more than 800 ligands. In addition, a number of web-based tools have been integrated to facilitate users in creating web logo and two-sample between ligand interacting and non-interacting residues.ConclusionsIn summary, this manuscript presents a web-server for analysis of ligand interacting residue. This server is available for public use from URL http://crdd.osdd.net/raghava/lpicom.ReviewersThis article was reviewed by Prof Michael Gromiha, Prof Vladimir Poroikov and Prof Zlatko Trajanoski.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-016-0118-5) contains supplementary material, which is available to authorized users.

First Genome-Wide Mutational Map of Macaques a Boon for Evolutionary Studies.

Scientists have long used comparative animal studies to better understand the nuances of human evolution, from making diverse body plans to the emergence of entirely powerful and unique features structures, including the human brain. Now, in a new study appearing in the advanced online edition of Molecular Biology and Evolution, corresponding authors Perdomo-Sabogal et al. (2016) have explored the global gene regulator GABPa to better understand its influence as a master switch. Their work honed in on human-specific DNA regions potentially contributing to human evolution in key functions including cell energy, division and death, and diseases ranging from brain disorders to diabetes. Many of the GABPa gene regulatory sites were found to be significantly enriched at genes important for brain, central nervous system and spinal cord functions that are unique to humans, an area ripe for future research, and GABPa deficits linked to Alzheimer’s and Parkinson diseases, breast cancer, and autism. GABPa exerts its profound genomic influence by a tandem repeat of four highly conserved DNA letters. The GGAA consensus motif is found to control genes in all primates, mice, dogs, and cows, underscoring its importance in these animals. By performing ChIP-Seq experiments in a human cell line, the group found that GABPa exercises its influence across the genome at 11,619 putative GABPa binding sites representing nearly 4,000 genes. Through a tour de force sequence comparison of the human GABPa binding sites with similar sequences from 34 mammals, they identified 224 GABPa binding sites unique to humans. Comparing chimpanzee, macaque, and human cell lines and mutating the GGAA sequence to more ancestral states to determine how these changes affected GABPa function, they further narrowed down their function. For instance, when the human GABPa binding sites were introduced into primate cells, gene expression was increased up to five times in 17 out of 18 cases, indicating that a recent human mutation in GABPa resulted in higher gene expression and forces behind the evolution of novel, uniquely human functions. “Mutations that cause changes in the regulation of gene activities are one of the major factors in shaping species during evolution,” said corresponding author Robert Querfurth. “Our study demonstrates how, out of the millions of DNA regions, in which we differ from other apes, we can sift out those that, in response to a specific regulatory protein, cause gene activity changes in human cells. In this set, we find genes involved in brain and breast development and also in diseases like Alzheimer’s and Parkinson’s.”