Phosphate binding sites identification in protein structures

Nearly half of known protein structures interact with phosphate-containing ligands, such as nucleotides and other cofactors. Many methods have been developed for the identification of metal ions-binding sites and some for bigger ligands such as carbohydrates, but none is yet available for the prediction of phosphate-binding sites. Here we describe Pfinder, a method that predicts binding sites for phosphate groups, both in the form of ions or as parts of other non-peptide ligands, in proteins of known structure. Pfinder uses the Query3D local structural comparison algorithm to scan a protein structure for the presence of a number of structural motifs identified for their ability to bind the phosphate chemical group. Pfinder has been tested on a data set of 52 proteins for which both the apo and holo forms were available. We obtained at least one correct prediction in 63% of the holo structures and in 62% of the apo. The ability of Pfinder to recognize a phosphate-binding site in unbound protein structures makes it an ideal tool for functional annotation and for complementing docking and drug design methods. The Pfinder program is available at http://pdbfun.uniroma2.it/pfinder.     

In vivo identification and characterization of binding sites for selective serotonin reuptake inhibitors in mouse brain

The present study was undertaken to identify and characterize in vivo binding sites of selective serotonin reuptake inhibitors (SSRIs) in the mouse brain by using [3H]paroxetine as radioligand. Relatively higher concentration of [3H]paroxetine was detected in the whole brain (minus cerebellum) than in the plasma of mice after the i.v. injection of the radioligand, and the half-life (t1/2) of elimination was much slower. The in vivo specific [3H]paroxetine binding in the mouse brain after the i.v. injection was defined as the difference of particulate-bound radioactivity between the whole brain and cerebellum, and it was dose-dependently attenuated by oral or intraperitoneal administration of fluoxetine (8.68-116 micromol/kg). Furthermore, oral administration of fluvoxamine, fluoxetine, paroxetine and sertraline at the pharmacologically relevant doses reduced significantly (25-94%) in vivo specific [3H]paroxetine binding in the cerebral cortex, striatum, hippocampus, thalamus and midbrain of mice, and their significant decreases were observed up to at least 8 h (fluvoxamine), 24 h (fluoxetine), and 12 h (paroxetine and sertraline) later. The value of area under the curve (AUC) for decrease in [3H]paroxetine binding vs. time in each brain region was largest for fluoxetine among these SSRIs, due to the relatively longer-lasting occupation of brain serotonin transporter. The AUC value in mouse brain after oral administration of each SSRI was 1.2-3.2 times greater in the thalamus and midbrain than in the cerebral cortex, striatum and hippocampus. Thus, the present study has revealed that [3H]paroxetine may be a suitable radioligand for in vivo characterization of brain binding sites and pharmacological effects of SSRIs.     

Functional identification of catalytic metal ion binding sites within RNA

The viability of living systems depends inextricably on enzymes that catalyze phosphoryl transfer reactions. For many enzymes in this class, including several ribozymes, divalent metal ions serve as obligate cofactors. Understanding how metal ions mediate catalysis requires elucidation of metal ion interactions with both the enzyme and the substrate(s). In the Tetrahymena group I intron, previous work using atomic mutagenesis and quantitative analysis of metal ion rescue behavior identified three metal ions (MA, MB, and MC) that make five interactions with the ribozyme substrates in the reaction's transition state. Here, we combine substrate atomic mutagenesis with site-specific phosphorothioate substitutions in the ribozyme backbone to develop a powerful, general strategy for defining the ligands of catalytic metal ions within RNA. In applying this strategy to the Tetrahymena group I intron, we have identified the pro-SP phosphoryl oxygen at nucleotide C262 as a ribozyme ligand for MC. Our findings establish a direct connection between the ribozyme core and the functionally defined model of the chemical transition state, thereby extending the known set of transition-state interactions and providing information critical for the application of the recent group I intron crystallographic structures to the understanding of catalysis.     

In vivo functional analysis of in vitro protein binding sites in the immunoglobulin heavy chain enhancer.

We have systematically investigated the functional role of protein binding sites within the mouse immunoglobulin heavy chain enhancer which we previously identified by in vitro binding studies (1,2). Each binding site was deleted, mutant enhancers were cloned 3' of the chloramphenicol acetyl transferase gene in the vector pA10CAT2 and transfected into plasmacytoma cells. We demonstrate that the newly identified site E, located at 324-338 bp, is important for enhancer function; previously identified sites B(uE1), Cl(uE2), C2(uE3) and C3 were also shown to be important for enhancer activity. Sites A and D are not required for IgH enhancer function, as assayed by our methods. Thus, including the octamer site, six protein binding sites which bind at least six different proteins are important for enhancer function in vivo.     

In vivo selection of functional variations in essential sites of ribosomal RNA

The technique of "in vivo selection of functional ribosomes" is a genetic approach to dissecting the link between the structure and function of critical sites of rRNA. This method proceeds through selection of functional variants among cells that express ribosomes from a pool of rRNA-containing randomized sites. The selection of bacterial clones with functional ribosomes is based on the use of a plasmid carrying a rRNA operon in which a site of interest has been randomized and a point mutation conferring an antibiotic resistance has been introduced. Cells expressing functional ribosomes are then selected on medium containing the antibiotic. With this approach one can isolate at once all the possible variations at a given rRNA site that are able to sustain normal ribosome function. The identification of covariations in between several nucleotides that maintain wild-type ribosome activity can thus help demonstrate the function of specific interactions in rRNA.     

In vivo detection of regulatory factor binding sites of Arabidopsis thaliana Adh

In vivo footprinting experiments have been used to analyze the binding of trans-acting regulatory factors in the 5 flanking region upstream of the alcohol dehydrogenase (Adh) gene from Arabidopsis thaliana. Protein-DNA interactions were detected by dimethyl sulfate footprinting and genomic sequencing, using an A. thaliana cell suspension culture that constitutively expressed the Adh gene. Several distinct footprinting domains have been characterized, and the potential effects of the corresponding trans-acting factors have been inferred from a comparison with data from the maize alcohol dehydrogenase-1 (Adh1) gene. One binding site is similar in sequence to one of the anaerobic response elements (ARE) of the maize gene, which has also been shown to bind to a trans-acting factor. Several of the remaining binding sites apparently represent a class of elements sharing the sequence 5-GTGG-3 within their footprint.Comparisons with maize Adh1 in vivo protein interactions reveal that the elements of Adh promoter structure are highly conserved, but the relative and absolute positions of the elements are variable.     

Two distinct binding modes of a protein cofactor with its target RNA

Like most cellular RNA enzymes, the bI5 group I intron requires binding by a protein cofactor to fold correctly. Here, we use single-molecule approaches to monitor the structural dynamics of the bI5 RNA in real time as it assembles with its CBP2 protein cofactor. These experiments show that CBP2 binds to the target RNA in two distinct modes with apparently opposite effects: a "non-specific" mode that forms rapidly and induces large conformational fluctuations in the RNA, and a "specific" mode that forms slowly and stabilizes the native RNA structure. The bI5 RNA folds though multiple pathways toward the native state, typically traversing dynamic intermediate states induced by non-specific binding of CBP2. These results suggest that the protein cofactor-assisted RNA folding involves sequential non-specific and specific protein-RNA interactions. The non-specific interaction potentially increases the local concentration of CBP2 and the number of conformational states accessible to the RNA, which may promote the formation of specific RNA-protein interactions.     

Independent binding sites of small protein B onto transfer-messenger RNA during trans-translation

Stalled bacterial ribosomes are freed by transfer-messenger RNA (tmRNA). With the help of small protein B (SmpB), protein synthesis restarts and tmRNA adds a tag to the stalled protein for destruction. The conformation of a 347 nt long tmRNA from a thermophile and its interactions with SmpB were monitored using structural probes. The RNA is highly folded, including the reading frame, with <30% of unpaired residues. Footprints between SmpB and tmRNA are in the elbow of the tRNA domain, in some pseudoknots including one essential for function and in the lower part of the stem exiting the tRNA domain. The footprints outside the tRNA domain are scattered onto the tmRNA sequence, but form a cluster onto its tertiary structure derived from cryo-EM data. Some footprints flank the first triplet to be translated in tmRNA, suggesting that SmpB participates in the insertion of the tmRNA-encoded reading frame into the decoding center. To discriminate between a conformational rearrangement of tmRNA and independent binding sites, surface plasmon resonance was used and has identified three independent binding sites of SmpB on the RNA, including the site on the tRNA domain. Accordingly, SmpB is proposed to move on the tmRNA scaffold during trans-translation.     

Using binding profiles to predict binding sites of target RNAs

Prediction of RNA-RNA interaction is a key to elucidating possible functions of small non-coding RNAs, and a number of computational methods have been proposed to analyze interacting RNA secondary structures. In this article, we focus on predicting binding sites of target RNAs that are expected to interact with regulatory antisense RNAs in a general form of interaction. For this purpose, we propose bistaRNA, a novel method for predicting multiple binding sites of target RNAs. bistaRNA employs binding profiles that represent scores for hybridized structures, leading to reducing the computational cost for interaction prediction. bistaRNA considers an ensemble of equilibrium interacting structures and seeks to maximize expected accuracy using dynamic programming. Experimental results on real interaction data validate good accuracy and fast computation time of bistaRNA as compared with several competitive methods. Moreover, we aim to find new targets given specific antisense RNAs, which provides interesting insights into antisense RNA regulation. bistaRNA is implemented in C++. The program and Supplementary Material are available at http://rna.naist.jp/program/bistarna/.     

High-throughput identification of interacting protein-protein binding sites

Background  

With the advent of increasing sequence and structural data, a number of methods have been proposed to locate putative protein binding sites from protein surfaces. Therefore, methods that are able to identify whether these binding sites interact are needed.     

In vivo DNA binding of bacteriophage GA-1 protein p6

Bacteriophage GA-1 infects Bacillus sp. strain G1R and has a linear double-stranded DNA genome with a terminal protein covalently linked to its 5′ ends. GA-1 protein p6 is very abundant in infected cells and binds DNA with no sequence specificity. We show here that it binds in vivo to the whole viral genome, as detected by cross-linking, chromatin immunoprecipitation, and real-time PCR analyses, and has the characteristics of a histone-like protein. Binding to DNA of GA-1 protein p6 shows little supercoiling dependency, in contrast to the ortholog protein of the evolutionary related Bacillus subtilis phage 29. This feature is a property of the protein rather than the DNA or the cellular background, since 29 protein p6 shows supercoiling-dependent binding to GA-1 DNA in Bacillus sp. strain G1R. GA-1 DNA replication is impaired in the presence of the gyrase inhibitors novobiocin and nalidixic acid, which indicates that, although noncovalently closed, the viral genome is topologically constrained in vivo. GA-1 protein p6 is also able to bind 29 DNA in B. subtilis cells; however, as expected, the binding is less supercoiling dependent than the one observed with the 29 protein p6. In addition, the nucleoprotein complex formed is not functional, since it is not able to transcomplement the DNA replication deficiency of a 29 sus6 mutant. Furthermore, we took advantage of 29 protein p6 binding to GA-1 DNA to find that the viral DNA ejection mechanism seems to take place, as in the case of 29, with a right to left polarity in a two-step, push-pull process.     

In vivo threonine phosphorylation of immunoglobulin binding protein (BiP) maps to its protein binding domain

lmmunoglobin binding protein (BiP) molecules exist as both monomers and oligomers and phosphorylated BiP is restricted to the oligomeric pool. Modified BiP is not bound to proteins such as immunoglobulin heavy chain and consequently, may constitute an inactive form. Unlike earlier analysis of mammalian BiP isolated by two-dimensional gel electrophoresis, results here demonstrated that immunoprecipitated BiP displayed predominantly threonine phosphorylation with only a trace of detectable phosphoserine. Like other Hsp70 family members, BiP is comprised of three domains: an amino terminal domain which binds nucleotide, an 18 kilodalton domain which binds peptide, and a carboxyl terminal variable domain of unknown function. Cyanogen bromide cleavage and enzymatic digestion experiments mapped threonine phosphorylation to a site within a 47 amino acid sequence of the peptide binding domain which contains seven threonine residues. Partial proteinase K digestion in the presence of ATP independently verified that the in vivo phosphorylation site of mammalian (BiP) is located within the peptide binding domain. Furthermore, phosphorylation did not impede BiPs ATP-induced conformational change. Thus, the peptide binding domain of BiP is phosphorylated on threonine residue(s) mapping to not more than two tryptic fragments within the peptide binding domain. This location on the molecule could explain why phosphorylated BiP is not detected bound to proteins in vivo.