Using RNase sequence specificity to refine the identification of RNA-protein binding regions

Massively parallel pyrosequencing is a high-throughput technology that can sequence hundreds of thousands of DNA/RNA fragments in a single experiment. Combining it with immunoprecipitation-based biochemical assays, such as cross-linking immunoprecipitation (CLIP), provides a genome-wide method to detect the sites at which proteins bind DNA or RNA. In a CLIP-pyrosequencing experiment, the resolutions of the detected protein binding regions are partially determined by the length of the detected RNA fragments (CLIP amplicons) after trimming by RNase digestion. The lengths of these fragments usually range from 50-70 nucleotides. Many genomic regions are marked by multiple RNA fragments. In this paper, we report an empirical approach to refine the localization of protein binding regions by using the distribution pattern of the detected RNA fragments and the sequence specificity of RNase digestion. We present two regions to which multiple amplicons map as examples to demonstrate this approach.     

ChIP-chip: considerations for the design, analysis, and application of genome-wide chromatin immunoprecipitation experiments

Chromatin immunoprecipitation (ChIP) is a well-established procedure to investigate interactions between proteins and DNA. Coupled with whole-genome DNA microarrays, ChIPS allow one to determine the entire spectrum of in vivo DNA binding sites for any given protein. The design and analysis of ChIP-microarray (also called ChIP-chip) experiments differ significantly from the conventions used for locus ChIP approaches and ChIP-chip experiments, and these differences require new methods of analysis. In this light, we review the design of DNA microarrays, the selection of controls, the level of repetition required, and other critical parameters for success in the design and analysis of ChIP-chip experiments, especially those conducted in the context of mammalian or other relatively large genomes.     

ChIPpeakAnno: a Bioconductor package to annotate ChIP-seq and ChIP-chip data

Background  

Chromatin immunoprecipitation (ChIP) followed by high-throughput sequencing (ChIP-seq) or ChIP followed by genome tiling array analysis (ChIP-chip) have become standard technologies for genome-wide identification of DNA-binding protein target sites. A number of algorithms have been developed in parallel that allow identification of binding sites from ChIP-seq or ChIP-chip datasets and subsequent visualization in the University of California Santa Cruz (UCSC) Genome Browser as custom annotation tracks. However, summarizing these tracks can be a daunting task, particularly if there are a large number of binding sites or the binding sites are distributed widely across the genome.     

Mapping convulsants' binding to the GABA-A receptor chloride ionophore: a proposed model for channel binding sites

Gamma-aminobutyric acid (GABA) type A receptors play a key role in brain inhibitory neurotransmission, and are ligand-activated chloride channels blocked by numerous convulsant ligands. Here we summarize data on binding of picrotoxin, tetrazoles, beta-lactams, bicyclophosphates, butyrolactones and neurotoxic pesticides to GABA-A ionophore, and discuss functional and structural overlapping of their binding sites. The paper reviews data on convulsants' binding sensitivity to different point mutations in ionophore-lining second trans-membrane domains of GABA-A subunits, and maps possible location of convulsants' sites within the chloride ionophore. We also discuss data on inhibition of glycine, glutamate, serotonin (5-HT3) and N-acetylcholine receptors by GABA-A channel blockers, and examine the applicability of this model to other homologous ionotropic receptors. Positioning various convulsant-binding sites within ionophore of GABA-A receptors, this model enables a better understanding of complex architectonics of ionotropic receptors, and may be used for developing new channel-modulating drugs.     

Chemoreception by a cell with age-dependent receptor binding sites

The authors determine the time-dependent ligand current into a spherical cell that is covered with a large number of age-dependent receptors. These receptors can be in either of two states: active (i.e., available for ligand binding) or inactive. An active receptor turns inactive upon binding a ligand, and it can reappear as active at some later time. The transition inactive----active is treated as a probabilistic process. The ligand distribution around the cell is determined analytically in terms of this distribution at the cell surface. A set of nonlinear integral equations is derived for the distribution at the cell surface, which is solved numerically. In this way the time-dependent ligand current into the cell as well as the average active receptor population at the cell surface are determined.     

Direct coupling of a G-protein to dihydropyridine binding sites

Electrophysiological data support the existence of GTP-binding proteins interacting with voltage dependent calcium channels. Along this line the present study investigates the effect of GMP-PNP, a stable GTP analogue, on the displacement of [3H]-PN 200-110 binding by agonist and antagonist dihydropyridines in synaptic membranes prepared from rat cortex. The results show that GMP-PNP increases the ability of the agonist dihydropyridine BAY K 8644 to displace [3H]-PN 200-110 binding. The in vivo treatment with Pertussis Toxin abolishes the effect produced by the non-hydrolysable GTP analogue.     

Heparin binding to antithrombin III: Variation in binding sites and affinity

Heparin fractions of different molecular weights and anticoagulant activities were prepared by chromatography on protamine-Sepharose, and the association constants and stoichiometry for binding to antithrombin III were determined by measurement of enhancement of tryptophan fluorescence. A 7,900 molecular weight heparin preparation bound to antithrombin III with a stoichiometry of close to 2:1, whereas 14,300 and 21,600 molecular weight fractions bound at approximately 1:1 with the protein. Apparent association constants were 0.66 × 10⁶ M^?1 for the low molecular weight preparation and 2.89 × 10⁶ M^?1 for the high molecular weight material. Maximal fluorescence enhancement was greater with the higher molecular weight heparin. These results suggest a model of heparin-antithrombin III binding in which two sites on antithrombin III can accommodate one large heparin molecule with high affinity or two smaller molecules with low affinity.     

ChIPOTle: a user-friendly tool for the analysis of ChIP-chip data

ChIPOTle (Chromatin ImmunoPrecipitation On Tiled arrays) takes advantage of two unique properties of ChIP-chip data: the single-tailed nature of the data, caused by specific enrichment but not specific depletion of genomic fragments; and the predictable enrichment of DNA fragments adjacent to sites of direct protein-DNA interaction. Implemented as a Microsoft Excel macro written in Visual Basic, ChIPOTle uses a sliding window approach that yields improvements in the identification of bona fide sites of protein-DNA interaction.     

Homologous RNA polymerase binding to Escherichia coli DNA: a study of the distribution of stable binding sites.

The number and the distribution of the sites of Escherichia coli DNA that form stable complexes with the homologous RNA polymerase (class A sites according to Hinkle and Chamberlin [3]) have been investigated. Almost all the DNA can bind RNA polymerase, even when fragmented at short (undergenic) size; this general non-promoter-specific binding is highly labile and is not temperature-dependent. The range of RNA polymerase/DNA ratios that give rise to the stable temperature-dependent complexes was examined. The amount and the distribution of class A complexes were studied analysing the dissociation of complexes formed by RNA polymerase on DNA fragments of various length. The E. coli genome appears to form 3.8 X 10(3) stable complexes; the majority of these complexes shows a short-range distribution (800-1200 base pairs). The rest is more widely spaced (1200-6000 base pairs).     

ProMate: a structure based prediction program to identify the location of protein-protein binding sites

Is the whole protein surface available for interaction with other proteins, or are specific sites pre-assigned according to their biophysical and structural character? And if so, is it possible to predict the location of the binding site from the surface properties? These questions are answered quantitatively by probing the surfaces of proteins using spheres of radius of 10 A on a database (DB) of 57 unique, non-homologous proteins involved in heteromeric, transient protein-protein interactions for which the structures of both the unbound and bound states were determined. In structural terms, we found the binding site to have a preference for beta-sheets and for relatively long non-structured chains, but not for alpha-helices. Chemically, aromatic side-chains show a clear preference for binding sites. While the hydrophobic and polar content of the interface is similar to the rest of the surface, hydrophobic and polar residues tend to cluster in interfaces. In the crystal, the binding site has more bound water molecules surrounding it, and a lower B-factor already in the unbound protein. The same biophysical properties were found to hold for the unbound and bound DBs. All the significant interface properties were combined into ProMate, an interface prediction program. This was followed by an optimization step to choose the best combination of properties, as many of them are correlated. During optimization and prediction, the tested proteins were not used for data collection, to avoid over-fitting. The prediction algorithm is fully automated, and is used to predict the location of potential binding sites on unbound proteins with known structures. The algorithm is able to successfully predict the location of the interface for about 70% of the proteins. The success rate of the predictor was equal whether applied on the unbound DB or on the disjoint bound DB. A prediction is assumed correct if over half of the predicted continuous interface patch is indeed interface. The ability to predict the location of protein-protein interfaces has far reaching implications both towards our understanding of specificity and kinetics of binding, as well as in assisting in the analysis of the proteome.     

Drug binding sites on P-glycoprotein are altered by ATP binding prior to nucleotide hydrolysis

P-glycoprotein (P-gp) confers multiple drug resistance on cancer cells by acting as a plasma membrane localized ATP-dependent drug efflux pump. Currently, there is little information on the nature of the communication between the energy-providing nucleotide binding domains (NBDs) and the drug binding sites of P-gp to generate transport of substrate. Many substrates and modulators cause alterations in ATP hydrolysis, but what effect do the various stages of the catalytic cycle have on drug interaction with P-gp? Vanadate trapping of Mg.ADP caused a reversible decrease in the binding capacity of the transported substrate [(3)H]-vinblastine and the nontransported modulator [(3)H]XR9576 to P-gp in CH(r)B30 cell membranes. The non-hydrolyzable nucleotide analogue ATP-gamma-S also caused a reduction in the binding capacity of [(3)H]-vinblastine but not for the modulator [(3)H]XR9576. This indicates that signaling to the NBDs following binding of a nontransported modulator is different to that transmitted upon interaction of a transported substrate. Second, it appears that the binding of nucleotide, rather than its hydrolysis, causes the initial conformational shift in the drug-binding site during a transport cycle.     

Evidence for two distinct effector-binding sites in threonine deaminase by site-directed mutagenesis, kinetic, and binding experiments

A three-dimensional structure comparison between the dimeric regulatory serine-binding domain of Escherichia coli D-3-phosphoglycerate dehydrogenase [Schuller, D. J., Grant, G. A., and Banaszak, L. J. (1995) Nat. Struct. Biol. 2, 69-76] and the regulatory domain of E. coli threonine deaminase [Gallagher, D. T., Gilliland, G. L., Xiao, G., Zondlo, J., Fisher, K. E., Chinchilla, D. , and Eisenstein, E. (1998) Structure 6, 465-475] led us to make the hypothesis that threonine deaminase could have two binding sites per monomer. To test this hypothesis about the corresponding plant enzyme, site-directed mutagenesis was carried out on the recombinant Arabidopsis thaliana threonine deaminase. Kinetic and binding experiments demonstrated for the first time that each regulatory domain of the monomers of A. thaliana threonine deaminase possesses two different effector-binding sites constituted in part by Y449 and Y543. Our results demonstrate that Y449 belongs to a high-affinity binding site whose interaction with a first isoleucine induces conformational modifications yielding a conformer displaying a higher activity and with enhanced ability to bind a second isoleucine on a lower-affinity binding site containing Y543. Isoleucine interaction with this latter binding site is responsible for conformational modifications leading to final inhibition of the enzyme. Y449 interacts with both regulators, isoleucine and valine. However, interaction of valine with the high-affinity binding site induces different conformational modifications leading to reversal of isoleucine binding and reversal of inhibition.     

High and low affinity opioid binding sites: relationship to mu and delta sites

Binding and pharmacological studies suggest a common opiate and enkephalin binding site in addition to their previously reported selective sites. This common high affinity site has tentatively been named mu1, distinguishing it from the morphine-selective site (mu2) and enkephalin-selective site (delta). The existence of this additional common high affinity site and its association with opiate and opioid peptide analgesia may help explain some pharmacological observations, such as the cross tolerance between morphine and enkephalin analgesia and the lack of cross tolerance between them in the guinea pig ileum and mouse vas deferens bioassays.     

Guided docking: first step to locate potential binding sites

The long-range electrostatic forces of the targets in round 2 of the Critical Assessment of PRediction of Interactions (CAPRI) experiment were examined and a simple guided docking method, based on these forces, was applied. The method described consists of calculating an initial rigid body trajectory and an optional final, fully flexible refinement stage. Although only limited success was found in predicting the final complexes, some interesting information was discovered. In particular, the long-range forces seem to give some insight into the unusual binding mode of target 4 while raising some questions about target 7, which warrant further investigation.     

Repression of dpp targets by binding of brinker to mad sites

Signaling by decapentaplegic (Dpp), a Drosophila member of the transforming growth factor (TGF) beta superfamily of growth factors, has recently been shown to activate targets such as vestigial (vg) indirectly through negative regulation of brinker (brk). Here we show that the Brk protein functions as a repressor by binding to Dpp response elements. The Brk DNA binding activity was localized to an amino-terminal region containing a putative homeodomain. Brk bound to a Dpp response element of the Ultrabithorax (Ubx) midgut enhancer at a sequence that overlaps a binding site for the Smad protein, Mothers Against Dpp (Mad). Furthermore, Brk was able to compete with Mad for occupancy of this binding site. This recognition of overlapping binding sites provides a potential explanation for why the G/C-rich Mad binding site consensus differs the Smad3/Smad4 binding site consensus. We also found that the Dpp response element from Ubx was more sensitive than the vg quadrant enhancer to repression by Brk. This difference correlates with short-range activation of Ubx by Dpp in the visceral mesoderm, whereas vg exhibits a long-range response to Dpp in the wing imaginal disc, indicating that Brk binding sites may play a critical role in limiting thresholds for activation by Dpp. Finally, we provide evidence that Brk is capable of functioning as an active repressor. Thus, whereas Brk and Mad compete for regulation of Ubx and vg, Brk may regulate other Dpp targets without direct involvement of Mad.     

Flavonoid-serum albumin complexation: determination of binding constants and binding sites by fluorescence spectroscopy

After a meal rich in plant products, dietary flavonols can be detected in plasma as serum albumin-bound conjugates. Flavonol-albumin binding is expected to modulate the bioavailability of flavonols. In this work, the binding of structurally different flavonoids to human and bovine serum albumins is investigated by fluorescence spectroscopy using three methods: the quenching of the albumin fluorescence, the enhancement of the flavonoid fluorescence, the quenching of the fluorescence of the quercetin-albumin complex by a second flavonoid. The latter method is extended to probes whose high-affinity binding sites are known to be located in one of the two major subdomains (warfarin and dansyl-L-asparagine for subdomain IIA, ibuprofen and diazepam for subdomain IIIA). Overall, flavonoids display moderate affinities for albumins (binding constants in the range 1-15 x 10(4) M(-1)), flavones and flavonols being most tightly bound. Glycosidation and sulfation could lower the affinity to albumin by one order of magnitude depending on the conjugation site. Despite multiple binding of both quercetin and site probes, it can be proposed that the binding of flavonols primarily takes place in subdomain IIA. Significant differences in affinity and binding location are observed for the highly homologous HSA and BSA.