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.     

Distribution of stable DnaA-binding sites on the Bacillus subtilis genome detected using a modified ChIP-chip method.

We developed a modified ChIP-chip method, designated ChAP-chip (Chromatin Affinity Precipitation coupled with tiling chip). The binding sites of Bacillus subtilis Spo0J determined using this technique were consistent with previous findings. A DNA replication initiator protein, DnaA, formed stable complexes at eight intergenic regions on the B. subtilis genome. Characterization of the binding sequences suggested that two factors -- the local density of DnaA boxes and their affinities for DnaA -- are critical for stable binding. We further showed that in addition to autoregulation, DnaA directly modulate the expression of sda in a positive, and ywlC and yydA in a negative manner. Examination of possible stable DnaA-binding sequences in other Bacillus species suggested that DnaA-dependent regulation of those genes is maintained in most bacteria examined, supporting their biological significance. In addition, a possible stable DnaA-binding site downstream of gcp is also suggested to be conserved. Furthermore, potential DnaA-binding sequences specific for each bacterium have been identified, generally in close proximity to oriC. These findings suggest that DnaA plays several additional roles, such as control of the level of effective initiator, ATP-DnaA, and/or stabilization of the domain structure of the genome around oriC for the proper initiation of chromosome replication.     

Alpha 2-macroglobulin binding to cultured fibroblasts: identification by affinity chromatography of high-affinity binding sites

Binding sites having the properties of high-affinity receptors for activated alpha 2-macroglobulin (alpha 2M) have been purified over 100-fold from membranes of spontaneously transformed NIH-3T3 cells (J. A. Hanover, S.-y. Cheng, M. C. Willingham, and I. H. Pastan [1983] J. Biol. Chem. 258, 370-377). To identify the molecular species involved in high-affinity binding, the solubilized receptor has been purified 500-fold by conventional procedures and further purified by affinity chromatography. After radioiodination of the 500-fold-purified preparation, the detergent-solubilized extract was applied to alpha 2M-Sepharose and an 85,000 +/- 5000 Mr species was selectively retained by the column. Binding of the 85,000 +/- 5000 Mr species to the affinity resin was inhibited by EDTA and by excess alpha 2M. Elution from the affinity column could be accomplished with bacitracin, a competitive inhibitor of alpha 2M binding, or with EDTA. Consistent with the previously reported characteristics of the high-affinity alpha 2M receptor, the 85,000 Mr species bound much more efficiently to methylamine-activated alpha 2M-Affigel than to alpha 2M-Affigel which had not been amine-activated. The present data suggest that a protein with a subunit Mr of 85,000 +/- 5000 may represent a component of the high-affinity alpha 2M receptor present on cultured fibroblasts.     

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.     

IgG binding to cytoskeletal intermediate filaments activates the complement cascade

The cellular plasma membrane becomes permeable to macromolecules during the cell injury process. This results in exposure of the interior of the cell to plasma proteins and to high-affinity binding of the Fc part of IgG to intermediate filaments (Hansson, G K, Starkebaum, G A, Benditt, E P & Schwartz, S M, Proc natl acad sci USA 81 (1984) 3103). Such IgG binding could be an early step in a process that serves to eliminate the injured cell. We have now identified its effect on the complement system. Intermediate filaments were reconstituted in vitro from purified vimentin, and incubated with plasma proteins. Cross-linker experiments showed binding of the heavy chain of IgG to vimentin, indicating that the vimentin protein carries an Fc-binding site. In contrast, no direct binding of complement factor Clq to vimentin could be detected. Binding of both IgG and Clq could, however, be detected by immunofluorescence when cytoskeletons of cultured endothelial cells were incubated with fresh serum. Therefore, IgG binding to filaments in the presence of serum is accompanied by Clq binding to IgG. This was in turn followed by fixation of C4 and C3 to intermediate filaments in a process that was dependent on both Ca2+, Mg2+ and Clq, indicating that it was part of a complement activation via the classical pathway. Exposure of fresh serum to intermediate filaments also resulted in production of the anaphylatoxic complement cleavage fragment. C3a, with a dose-response relationship between the amount of filaments present and the amount of C3a generated. Chemotactic activity towards granulocytes and monocytes was also generated by exposure of serum to intermediate filaments, and this activity was dependent on the presence of complement factor C5 and on the classical complement activation cascade, implying that it was due to the C5a peptide. Exposure of the interior of the cell to plasma proteins thus results in binding of IgG to intermediate filaments and activation of the complement cascade via the classical pathway. This, in turn generates bioactive mediators which may recruit leukocytes to the injured cell (C5a) and have profound effects on vascular permeability (C3a, C5a). We propose that this is part of a scavenger mechanism for the elimination of damaged cells.     

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.