Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Identification of eIF2Bγ and eIF2γ as cofactors of hepatitis C virus internal ribosome entry site-mediated translation using a functional genomics approach

The 5′-untranslated region of hepatitis C virus (HCV) is highly conserved, folds into a complex secondary structure, and functions as an internal ribosome entry site (IRES) to initiate translation of HCV proteins. We have developed a selection system based on a randomized hairpin ribozyme gene library to identify cellular factors involved in HCV IRES function. A retroviral vector ribozyme library with randomized target recognition sequences was introduced into HeLa cells, stably expressing a bicistronic construct encoding the hygromycin B phosphotransferase gene and the herpes simplex virus thymidine kinase gene (HSV-tk). Translation of the HSV-tk gene was mediated by the HCV IRES. Cells expressing ribozymes that inhibit HCV IRES-mediated translation of HSV-tk were selected via their resistance to both ganciclovir and hygromycin B. Two ribozymes reproducibly conferred the ganciclovir-resistant phenotype and were shown to inhibit IRES-mediated translation of HCV core protein but did not inhibit cap-dependent protein translation or cell growth. The functional targets of these ribozymes were identified as the gamma subunits of human eukaryotic initiation factors 2B (eIF2Bγ) and 2 (eIF2γ), respectively. The involvement of eIF2Bγ and eIF2γ in HCV IRES-mediated translation was further validated by ribozymes directed against additional sites within the mRNAs of these genes. In addition to leading to the identification of cellular IRES cofactors, ribozymes obtained from this cellular selection system could be directly used to specifically inhibit HCV viral translation, thereby facilitating the development of new antiviral strategies for HCV infection.

Pentameric assembly of a neuronal glutamate transporter

Freeze-fracture electron microscopy was used to study the structure of a human neuronal glutamate transporter (EAAT3). EAAT3 was expressed in Xenopus laevis oocytes, and its function was correlated with the total number of transporters in the plasma membrane of the same cells. Function was assayed as the maximum charge moved in response to a series of transmembrane voltage pulses. The number of transporters in the plasma membrane was determined from the density of a distinct 10-nm freeze-fracture particle, which appeared in the protoplasmic face only after EAAT3 expression. The linear correlation between EAAT3 maximum carrier-mediated charge and the total number of the 10-nm particles suggested that this particle represented functional EAAT3 in the plasma membrane. The cross-sectional area of EAAT3 in the plasma membrane (48 ± 5 nm2) predicted 35 ± 3 transmembrane α-helices in the transporter complex. This information along with secondary structure models (6–10 transmembrane α-helices) suggested an oligomeric state for EAAT3. EAAT3 particles were pentagonal in shape in which five domains could be identified. They exhibited fivefold symmetry because they appeared as equilateral pentagons and the angle at the vertices was 110°. Each domain appeared to contribute to an extracellular mass that projects ≈3 nm into the extracellular space. Projections from all five domains taper toward an axis passing through the center of the pentagon, giving the transporter complex the appearance of a penton-based pyramid. The pentameric structure of EAAT3 offers new insights into its function as both a glutamate transporter and a glutamate-gated chloride channel.

DNA annealing by Rad52 Protein is stimulated by specific interaction with the complex of replication protein A and single-stranded DNA

Homologous recombination in Saccharomyces cerevisiae depends critically on RAD52 function. In vitro, Rad52 protein preferentially binds single-stranded DNA (ssDNA), mediates annealing of complementary ssDNA, and stimulates Rad51 protein-mediated DNA strand exchange. Replication protein A (RPA) is a ssDNA-binding protein that is also crucial to the recombination process. Herein we report that Rad52 protein effects the annealing of RPA–ssDNA complexes, complexes that are otherwise unable to anneal. The ability of Rad52 protein to promote annealing depends on both the type of ssDNA substrate and ssDNA binding protein. RPA allows, but slows, Rad52 protein-mediated annealing of oligonucleotides. In contrast, RPA is almost essential for annealing of longer plasmid-sized DNA but has little effect on the annealing of poly(dT) and poly(dA), which are relatively long DNA molecules free of secondary structure. These results suggest that one role of RPA in Rad52 protein-mediated annealing is the elimination of DNA secondary structure. However, neither Escherichia coli ssDNA binding protein nor human RPA can substitute in this reaction, indicating that RPA has a second role in this process, a role that requires specific RPA–Rad52 protein interactions. This idea is confirmed by the finding that RPA, which is complexed with nonhomologous ssDNA, inhibits annealing but the human RPA–ssDNA complex does not. Finally, we present a model for the early steps of the repair of double-strand DNA breaks in yeast.

Overexpression of a glutamate receptor (GluR2) ligand binding domain in Escherichia coli: Application of a novel protein folding screen

Expression of the S1S2 ligand binding domain [Kuusinen, A., Arvola, M. & Keinänen, K. (1995) EMBO J. 14, 6327–6332] of the rat α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid-selective glutamate receptor GluR2 in Escherichia coli under control of a T7 promoter leads to production of >100 mg/liter of histidine-tagged S1S2 protein (HS1S2) in the form of inclusion bodies. Using a novel fractional factorial folding screen and a rational, step-by-step approach, multiple conditions were determined for the folding of the HS1S2 α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid binding domain. Characterization of the HS1S2 ligand binding domain showed that it is water-soluble, monomeric, has significant secondary structure, and is sensitive to trypsinolysis at sites close to the beginning of the putative transmembrane regions. Application of a fractional factorial folding screen to other proteins may provide a useful means to evaluate E. coli as an economical and convenient expression host.

Escherichia coli RNA polymerase terminates transcription efficiently at rho-independent terminators on single-stranded DNA templates

Several models have been proposed for the mechanism of transcript termination by Escherichia coli RNA polymerase at rho-independent terminators. Yager and von Hippel (Yager, T. D. & von Hippel, P. H. (1991) Biochemistry 30, 1097–118) postulated that the transcription complex is stabilized by enzyme–nucleic acid interactions and the favorable free energy of a 12-bp RNA–DNA hybrid but is destabilized by the free energy required to maintain an extended transcription bubble. Termination, by their model, is viewed simply as displacement of the RNA transcript from the hybrid helix by reformation of the DNA helix. We have proposed an alternative model where the RNA transcript is stably bound to RNA polymerase primarily through interactions with two single-strand specific RNA-binding sites; termination is triggered by formation of an RNA hairpin that reduces binding of the RNA to one RNA-binding site and, ultimately, leads to its ejection from the complex. To distinguish between these models, we have tested whether E. coli RNA polymerase can terminate transcription at rho-independent terminators on single-stranded DNA. RNA polymerase cannot form a transcription bubble on these templates; thus, the Yager–von Hippel model predicts that intrinsic termination will not occur. We find that transcript elongation on single-stranded DNA templates is hindered somewhat by DNA secondary structure. However, E. coli RNA polymerase efficiently terminates and releases transcripts at several rho-independent terminators on such templates at the same positions as termination occurs on duplex DNAs. Therefore, neither the nontranscribed DNA strand nor the transcription bubble is essential for rho-independent termination by E. coli RNA polymerase.

Did homeodomain proteins duplicate before the origin of angiosperms, fungi, and metazoa?

Homeodomain proteins are transcription factors that play a critical role in early development in eukaryotes. These proteins previously have been classified into numerous subgroups whose phylogenetic relationships are unclear. Our phylogenetic analysis of representative eukaryotic sequences suggests that there are two major groups of homeodomain proteins, each containing sequences from angiosperms, metazoa, and fungi. This result, based on parsimony and neighbor-joining analyses of primary amino acid sequences, was supported by two additional features of the proteins. The two protein groups are distinguished by an insertion/deletion in the homeodomain, between helices I and II. In addition, an amphipathic alpha-helical secondary structure in the region N terminal of the homeodomain is shared by angiosperm and metazoan sequences in one group. These results support the hypothesis that there was at least one duplication of homeobox genes before the origin of angiosperms, fungi, and metazoa. This duplication, in turn, suggests that these proteins had diverse functions early in the evolution of eukaryotes. The shared secondary structure in angiosperm and metazoan sequences points to an ancient conserved functional domain.

Thermostable and site-specific DNA binding of the gene product ORF56 from the Sulfolobus islandicus plasmid pRN1, a putative archael plasmid copy control protein

There is still a lack of information on the specific characteristics of DNA-binding proteins from hyperthermophiles. Here we report on the product of the gene orf56 from plasmid pRN1 of the acidophilic and thermophilic archaeon Sulfolobus islandicus. orf56 has not been characterised yet but low sequence similarily to several eubacterial plasmid-encoded genes suggests that this 6.5 kDa protein is a sequence-specific DNA-binding protein. The DNA-binding properties of ORF56, expressed in Escherichia coli, have been investigated by EMSA experiments and by fluorescence anisotropy measurements. Recombinant ORF56 binds to double-stranded DNA, specifically to an inverted repeat located within the promoter of orf56. Binding to this site could down-regulate transcription of the orf56 gene and also of the overlapping orf904 gene, encoding the putative initiator protein of plasmid replication. By gel filtration and chemical crosslinking we have shown that ORF56 is a dimeric protein. Stoichiometric fluorescence anisotropy titrations further indicate that ORF56 binds as a tetramer to the inverted repeat of its target binding site. CD spectroscopy points to a significant increase in ordered secondary structure of ORF56 upon binding DNA. ORF56 binds without apparent cooperativity to its target DNA with a dissociation constant in the nanomolar range. Quantitative analysis of binding isotherms performed at various salt concentrations and at different temperatures indicates that approximately seven ions are released upon complex formation and that complex formation is accompanied by a change in heat capacity of –6.2 kJ/mol. Furthermore, recombinant ORF56 proved to be highly thermostable and is able to bind DNA up to 85°C.

Protein–RNA interactions: a structural analysis

A detailed computational analysis of 32 protein–RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein–double-stranded DNA and protein–single-stranded DNA complexes. The interface properties of the protein–RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein–RNA and protein–DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein–RNA complexes, backbone contacts were more dominant in the protein–DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.

Conserved stem–loop structures in the HIV-1 RNA region containing the A3 3′ splice site and its cis-regulatory element: possible involvement in RNA splicing

The HIV-1 transcript is alternatively spliced to over 30 different mRNAs. Whether RNA secondary structure can influence HIV-1 RNA alternative splicing has not previously been examined. Here we have determined the secondary structure of the HIV-1/BRU RNA segment, containing the alternative A3, A4a, A4b, A4c and A5 3′ splice sites. Site A3, required for tat mRNA production, is contained in the terminal loop of a stem–loop structure (SLS2), which is highly conserved in HIV-1 and related SIVcpz strains. The exon splicing silencer (ESS2) acting on site A3 is located in a long irregular stem–loop structure (SLS3). Two SLS3 domains were protected by nuclear components under splicing condition assays. One contains the A4c branch points and a putative SR protein binding site. The other one is adjacent to ESS2. Unexpectedly, only the 3′ A residue of ESS2 was protected. The suboptimal A3 polypyrimidine tract (PPT) is base paired. Using site-directed mutagenesis and transfection of a mini-HIV-1 cDNA into HeLa cells, we found that, in a wild-type PPT context, a mutation of the A3 downstream sequence that reinforced SLS2 stability decreased site A3 utilization. This was not the case with an optimized PPT. Hence, sequence and secondary structure of the PPT may cooperate in limiting site A3 utilization.

SRPDB (Signal Recognition Particle Database)

Signal recognition particle (SRP) is a stable cytoplasmic ribonucleoprotein complex that serves to translocate secretory proteins across membranes during translation. The SRP Database (SRPDB) provides compilations of SRP components, ordered alphabetically and phylogenetically. Alignments emphasize phylogenetically-supported base pairs in SRP RNA and conserved residues in the proteins. Data are provided in various formats including a column arrangement for improved access and simplified computational usability. Included are motifs for identification of new sequences, SRP RNA secondary structure diagrams, 3-D models and links to high-resolution structures. This release includes 11 new SRP RNA sequences (total of 129), two protein SRP9 sequences (total of seven), two protein SRP14 sequences (total of 10), two protein SRP19 sequences (total of 16), 10 new SRP54 (ffh) sequences (total of 66), two protein SRP68 sequences (total of seven) and two protein SRP72 sequences (total of nine). Seven sequences of the SRP receptor α-subunit and its FtsY homolog (total of 51) are new. Also considered are β-subunit of SRP receptor, Flhf, Hbsu, CaM kinase II and cpSRP43. Access to SRPDB is at http://psyche.uthct.edu/dbs/SRPDB/SRPDB.html and the European mirror http://www.medkem.gu.se/dbs/SRPDB/SRPDB.html

PDBsum: summaries and analyses of PDB structures

PDBsum is a web-based database providing a largely pictorial summary of the key information on each macromolecular structure deposited at the Protein Data Bank (PDB). It includes images of the structure, annotated plots of each protein chain’s secondary structure, detailed structural analyses generated by the PROMOTIF program, summary PROCHECK results and schematic diagrams of protein–ligand and protein–DNA interactions. RasMol scripts highlight key aspects of the structure, such as the protein’s domains, PROSITE patterns and protein–ligand interactions, for interactive viewing in 3D. Numerous links take the user to related sites. PDBsum is updated whenever any new structures are released by the PDB and is freely accessible via http://www.biochem.ucl.ac.uk/bsm/pdbsum.

The European Large Subunit Ribosomal RNA Database

The European Large Subunit Ribosomal RNA Database compiles all complete or nearly complete large subunit ribosomal RNA sequences available from public sequence databases. These are provided in aligned format and the secondary structure, as derived by comparative sequence analysis, is included. Additional information about the sequences such as literature references and taxonomic information is also included. The database is available from our WWW server at http://rrna.uia.ac.be/lsu/.

Transverse relaxation-optimized NMR spectroscopy with the outer membrane protein OmpX in dihexanoyl phosphatidylcholine micelles

The 2H,13C,15N-labeled, 148-residue integral membrane protein OmpX from Escherichia coli was reconstituted with dihexanoyl phosphatidylcholine (DHPC) in mixed micelles of molecular mass of about 60 kDa. Transverse relaxation-optimized spectroscopy (TROSY)-type triple resonance NMR experiments and TROSY-type nuclear Overhauser enhancement spectra were recorded in 2 mM aqueous solutions of these mixed micelles at pH 6.8 and 30°C. Complete sequence-specific NMR assignments for the polypeptide backbone thus have been obtained. The 13C chemical shifts and the nuclear Overhauser effect data then resulted in the identification of the regular secondary structure elements of OmpX/DHPC in solution and in the collection of an input of conformational constraints for the computation of the global fold of the protein. The same type of polypeptide backbone fold is observed in the presently determined solution structure and the previously reported crystal structure of OmpX determined in the presence of the detergent n-octyltetraoxyethylene. Further structure refinement will have to rely on the additional resonance assignment of partially or fully protonated amino acid side chains, but the present data already demonstrate that relaxation-optimized NMR techniques open novel avenues for studies of structure and function of integral membrane proteins.

RNA editing in Trypanosoma brucei: characterization of gRNA U-tail interactions with partially edited mRNA substrates

Guide RNAs (gRNAs), key components of the RNA editing reaction in Trypanosoma brucei, direct the insertion and deletion of uridylate (U) residues. Analyses of gRNAs reveal three functional elements. The 5′-end of the gRNA contains the anchor, which is responsible for selection and binding to the pre-edited mRNA. The second element (the guiding region) provides the information required for editing. At the 3′-end of the gRNA is a non-encoded U-tail, whose function remains unclear. However, the cleavage–ligation model for editing proposes that the U-tail binds to purine-rich regions upstream of editing sites, thereby strengthening the interaction and holding onto the 5′ cleavage product. Our previous studies demonstrated that the U-tail interacts with upstream sequences and may play roles in both stabilization and tethering. These studies also indicated that the U-tail interactions involved mRNA regions that were to be subsequently edited. This raised the question of what happens to the mRNA–U-tail interaction as editing proceeds in the 3′→5′ direction. We examined gCYb-558 and its U-tail interaction with 5′CYbUT and two partially edited 5′CYb substrates. Our results indicate that the 3′-end of the U-tail interacts with the same sequence in all three mRNAs. Predicted secondary structures using crosslinking data suggest that a similar structure is maintained as editing proceeds. These results indicate that the role of the U-tail may also involve maintenance of important secondary structure motifs.