Assessing the structural conservation of protein pockets to study functional and allosteric sites: implications for drug discovery

Background: With the classical, active-site oriented drug-development approach reaching its limits, protein ligand-binding sites in general and allosteric sites in particular are increasingly attracting the interest of medicinal chemists in the search for new types of targets and strategies to drug development. Given that allostery represents one of the most common and powerful means to regulate protein function, the traditional drug discovery approach of targeting active sites can be extended by targeting allosteric or regulatory protein pockets that may allow the discovery of not only novel drug-like inhibitors, but activators as well. The wealth of available protein structural data can be exploited to further increase our understanding of allosterism, which in turn may have therapeutic applications. A first step in this direction is to identify and characterize putative effector sites that may be present in already available structural data. Results: We performed a large-scale study of protein cavities as potential allosteric and functional sites, by integrating publicly available information on protein sequences, structures and active sites for more than a thousand protein families. By identifying common pockets across different structures of the same protein family we developed a method to measure the pocket's structural conservation. The method was first parameterized using known active sites. We characterized the predicted pockets in terms of sequence and structural conservation, backbone flexibility and electrostatic potential. Although these different measures do not tend to correlate, their combination is useful in selecting functional and regulatory sites, as a detailed analysis of a handful of protein families shows. We finally estimated the numbers of potential allosteric or regulatory pockets that may be present in the data set, finding that pockets with putative functional and effector characteristics are widespread across protein families. Conclusions: Our results show that structurally conserved pockets are a common feature of protein families. The structural conservation of protein pockets, combined with other characteristics, can be exploited in drug discovery procedures, in particular for the selection of the most appropriate target protein and pocket for the design of drugs against entire protein families or subfamilies (e.g. for the development of broad-spectrum antimicrobials) or against a specific protein (e.g. in attempting to reduce side effects)

dnaMATE: a consensus melting temperature prediction server for short DNA sequences

An accurate and robust large-scale melting temperature prediction server for short DNA sequences is dispatched. The server calculates a consensus melting temperature value using the nearest-neighbor model based on three independent thermodynamic data tables. The consensus method gives an accurate prediction of melting temperature, as it has been recently demonstrated in a benchmark performed using all available experimental data for DNA sequences within the length range of 16–30 nt. This constitutes the first web server that has been implemented to perform a large-scale calculation of melting temperatures in real time (up to 5000 DNA sequences can be submitted in a single run). The expected accuracy of calculations carried out by this server in the range of 50–600 mM monovalent salt concentration is that 89% of the melting temperature predictions will have an error or deviation of <5°C from experimental data. The server can be freely accessed at . The standalone executable versions of this software for LINUX, Macintosh and Windows platforms are also freely available at the same web site. Detailed further information supporting this server is available at the same web site referenced above

Evolutionary potentials: structure specific knowledge-based potentials exploiting the evolutionary record of sequence homologs

So-called ‘Evolutionary potentials’ for protein structure prediction are derived using a single experimental protein structure and all three-dimensional models of its homologous sequences

antibacTR : dynamic antibacterial-drug-target ranking integrating comparative genomics, structural analysis and experimental annotation

Background : development of novel antibacterial drugs is both an urgent healthcare necessity and a partially neglected field. The last decades have seen a substantial decrease in the discovery of novel antibiotics, which combined with the recent thrive of multi-drug-resistant pathogens have generated a scenario of general concern. The procedures involved in the discovery and development of novel antibiotics are economically challenging, time consuming and lack any warranty of success. Furthermore, the return-on-investment for an antibacterial drug is usually marginal when compared to other therapeutics, which in part explains the decrease of private investment. - Results : in this work we present antibacTR, a computational pipeline designed to aid researchers in the selection of potential drug targets, one of the initial steps in antibacterial-drug discovery. The approach was designed and implemented as part of two publicly funded initiatives aimed at discovering novel antibacterial targets, mechanisms and drugs for a priority list of Gram-negative pathogens: Acinetobacter baumannii, Escherichia coli, Helicobacter pylori, Pseudomonas aeruginosa and Stenotrophomonas maltophilia. However, at present this list has been extended to cover a total of 74 fully sequenced Gram-negative pathogens. antibacTR is based on sequence comparisons and queries to multiple databases (e.g. gene essentiality, virulence factors) to rank proteins according to their potential as antibacterial targets. The dynamic ranking of potential drug targets can easily be executed, customized and accessed by the user through a web interface which also integrates computational analyses performed in-house and visualizable on-site. These include three-dimensional modeling of protein structures and prediction of active sites among other functionally relevant ligand-binding sites. - Conclusions : given its versatility and ease-of-use at integrating both experimental annotation and computational analyses, antibacTR may effectively assist microbiologists, medicinal-chemists and other researchers working in the field of antibacterial drug-discovery. The public web-interface for antibacTR is available at 'http://bioinf.uab.cat/antibactr'

Large-scale Oscillation of Structure-Related DNA Sequence Features in Human Chromosome 21

Human chromosome 21 is the only chromosome in human genome that exhibits oscillation of (G+C)-content of cycle length of hundreds kilobases (500 kb near the right telomere). We aim at establishing the existence of similar periodicity in structure-related sequence features in order to relate this (G+C)% oscillation to other biological phenomena. The following quantities are shown to oscillate with the same 500kb periodicity in human chromosome 21: binding energy calculated by two sets of dinucleotide-based thermodynamic parameters, AA/TT and AAA/TTT bi-/tri-nucleotide density, 5'-TA-3' dinucleotide density, and signal for 10/11-base periodicity of AA/TT or AAA/TTT. These intrinsic quantities are related to structural features of the double helix of DNA molecules, such as base-pair binding, untwisting/unwinding, stiffness, and a putative tendency for nucleosome formation.Comment: submitted to Physical Review

A Random Forest model for predicting allosteric and functional sites on proteins

We thank the Scottish Universities Life Sciences Alliance (SULSA) for funding to JBOM and for PB’s PhD studentship under NJW’s supervision.We created a computational method to identify allosteric sites using a machine learning method trained and tested on protein structures containing bound ligand molecules. The Random Forest machine learning approach was adopted to build our three-way predictive model. Based on descriptors collated for each ligand and binding site, the classification model allows us to assign protein cavities as allosteric, regular or orthosteric, and hence to identify allosteric sites. 43 structural descriptors per complex were derived and were used to characterize individual protein-ligand binding sites belonging to the three classes, allosteric, regular and orthosteric. We carried out a separate validation on a further unseen set of protein structures containing the ligand 2-(N-cyclohexylamino) ethane sulfonic acid (CHES).PostprintPeer reviewe

3did Update: domain–domain and peptide-mediated interactions of known 3D structure

The database of 3D interacting domains (3did) is a collection of protein interactions for which high-resolution 3D structures are known. 3did exploits structural information to provide the crucial molecular details necessary for understanding how protein interactions occur. Besides interactions between globular domains, the new release of 3did also contains a hand-curated set of transient peptide-mediated interactions. The interactions are grouped in Interaction Types, based on the mode of binding, and the different binding interfaces used in each type are also identified and catalogued. A web-based tool to query 3did is available at http://3did.irbbarcelona.org

Development and application of the human intestinal tract chip, a phylogenetic microarray : analysis of universally conserved phylotypes in the abundant microbiota of young and elderly adults

"P>In this paper we present the in silico assessment of the diversity of variable regions of the small subunit ribosomal RNA (SSU rRNA) gene based on an ecosystem-specific curated database, describe a probe design procedure based on two hypervariable regions with minimal redundancy and test the potential of such probe design strategy for the design of a flexible microarray platform. This resulted in the development and application of a phylogenetic microarray for studying the human gastrointestinal microbiota - referred as the human intestinal tract chip (HITChip). Over 4800 dedicated tiling oligonucleotide probes were designed based on two hypervariable regions of the SSU rRNA gene of 1140 unique microbial phylotypes (In this paper we present the in silico assessment of the diversity of variable regions of the small subunit ribosomal RNA (SSU rRNA) gene based on an ecosystem-specific curated database, describe a probe design procedure based on two hypervariable regions with minimal redundancy and test the potential of such probe design strategy for the design of a flexible microarray platform. This resulted in the development and application of a phylogenetic microarray for studying the human gastrointestinal microbiota - referred as the human intestinal tract chip (HITChip). Over 4800 dedicated tiling oligonucleotide probes were designed based on two hypervariable regions of the SSU rRNA gene of 1140 unique microbial phylotypes (In this paper we present the in silico assessment of the diversity of variable regions of the small subunit ribosomal RNA (SSU rRNA) gene based on an ecosystem-specific curated database, describe a probe design procedure based on two hypervariable regions with minimal redundancy and test the potential of such probe design strategy for the design of a flexible microarray platform. This resulted in the development and application of a phylogenetic microarray for studying the human gastrointestinal microbiota - referred as the human intestinal tract chip (HITChip). Over 4800 dedicated tiling oligonucleotide probes were designed based on two hypervariable regions of the SSU rRNA gene of 1140 unique microbial phylotypes (< 98% identity) following analysis of over 16 000 human intestinal SSU rRNA sequences. These HITChip probes were hybridized to a diverse set of human intestinal samples and SSU rRNA clones to validate its fingerprinting and quantification potential. Excellent reproducibility (median Pearson's correlation of 0.99) was obtained following hybridization with T7 polymerase transcripts generated in vitro from SSU rRNA gene amplicons. A linear dose-response was observed with artificial mixtures of 40 different representative amplicons with relative abundances as low as 0.1% of total microbiota. Analysis of three consecutively collected faecal samples from ten individuals (five young and five elderly adults) revealed temporal dynamics and confirmed that the adult intestinal microbiota is an individual-specific and relatively stable ecosystem. Further analysis of the stable part allowed for the identification of a universal microbiota core at the approximate genus level (90% sequence similarity). This core consists of members of Actinobacteria, Bacteroidetes and Firmicutes. Used as a phylogenetic fingerprinting tool with the possibility for relative quantification, the HITChip has the potential to bridge the gaps in our knowledge in the quantitative and qualitative description of the human gastrointestinal microbiota composition."Peer reviewe

Thermodynamic stability and Watson–Crick base pairing in the seed duplex are major determinants of the efficiency of the siRNA-based off-target effect

Short interfering RNA (siRNA) may down-regulate many unintended genes whose transcripts possess complementarity to the siRNA seed region, which contains 7 nt. The capability of siRNA to induce this off-target effect was highly correlated with the calculated melting temperature or standard free-energy change for formation of protein-free seed duplex, indicating that thermodynamic stability of seed duplex formed between the seed and target is one of the major factor in determining the degree of off-target effects. Furthermore, unlike intended gene silencing (RNA interference), off-target effect was completely abolished by introduction of a G:U pair into the seed duplex, and this loss in activity was completely recovered by a second mutation regenerating Watson–Crick pairing, indicating that seed duplex Watson–Crick pairing is also essential for off-target gene silencing. The off-target effect was more sensitive to siRNA concentration compared to intended gene silencing, which requires a near perfect sequence match between the siRNA guide strand and target mRNA

PrimerHunter: a primer design tool for PCR-based virus subtype identification

Rapid and reliable virus subtype identification is critical for accurate diagnosis of human infections, effective response to epidemic outbreaks and global-scale surveillance of highly pathogenic viral subtypes such as avian influenza H5N1. The polymerase chain reaction (PCR) has become the method of choice for virus subtype identification. However, designing subtype-specific PCR primer pairs is a very challenging task: on one hand, selected primer pairs must result in robust amplification in the presence of a significant degree of sequence heterogeneity within subtypes, on the other, they must discriminate between the subtype of interest and closely related subtypes. In this article, we present a new tool, called PrimerHunter, that can be used to select highly sensitive and specific primers for virus subtyping. Our tool takes as input sets of both target and nontarget sequences. Primers are selected such that they efficiently amplify any one of the target sequences, and none of the nontarget sequences. PrimerHunter ensures the desired amplification properties by using accurate estimates of melting temperature with mismatches, computed based on the nearest neighbor model via an efficient fractional programming algorithm. Validation experiments with three avian influenza HA subtypes confirm that primers selected by PrimerHunter have high sensitivity and specificity for target sequences