A structural census of the current population of protein sequences

We examine the occurrence of the ≈300 known protein folds in different groups of organisms. To do this, we characterize a large fraction of the currently known protein sequences (≈140,000) in structural terms, by matching them to known structures via sequence comparison (or by secondary-structure class prediction for those without structural homologues). Overall, we find that an appreciable fraction of the known folds are present in each of the major groups of organisms (e.g., bacteria and eukaryotes share 156 of 275 folds), and most of the common folds are associated with many families of nonhomologous sequences (i.e., >10 sequence families for each common fold). However, different groups of organisms have characteristically distinct distributions of folds. So, for instance, some of the most common folds in vertebrates, such as globins or zinc fingers, are rare or absent in bacteria. Many of these differences in fold usage are biologically reasonable, such as the folds of metabolic enzymes being common in bacteria and those associated with extracellular transport and communication being common in animals. They also have important implications for database-based methods for fold recognition, suggesting that an unknown sequence from a plant is more likely to have a certain fold (e.g., a TIM barrel) than an unknown sequence from an animal.

Three novel families of miniature inverted-repeat transposable elements are associated with genes of the yellow fever mosquito, Aedes aegypti

Three novel families of transposable elements, Wukong, Wujin, and Wuneng, are described in the yellow fever mosquito, Aedes aegypti. Their copy numbers range from 2,100 to 3,000 per haploid genome. There are high degrees of sequence similarity within each family, and many structural but not sequence similarities between families. The common structural characteristics include small size, no coding potential, terminal inverted repeats, potential to form a stable secondary structure, A+T richness, and putative 2- to 4-bp A+T-biased specific target sites. Evidence of previous mobility is presented for the Wukong elements. Elements of these three families are associated with 7 of 16 fully or partially sequenced Ae. aegypti genes. Characteristics of these mosquito elements indicate strong similarities to the miniature inverted-repeat transposable elements (MITEs) recently found to be associated with plant genes. MITE-like elements have also been reported in two species of Xenopus and in Homo sapiens. This characterization of multiple families of highly repetitive MITE-like elements in an invertebrate extends the range of these elements in eukaryotic genomes. A hypothesis is presented relating genome size and organization to the presence of highly reiterated MITE families. The association of MITE-like elements with Ae. aegypti genes shows the same bias toward noncoding regions as in plants. This association has potentially important implications for the evolution of gene regulation.

Two-dimensional IR correlation spectroscopy: Sequential events in the unfolding process of the λ Cro-V55C repressor protein

A question often posed in protein folding/unfolding studies is whether the process is fully cooperative or whether it contains sequential elements. To address this question, one needs tools capable of resolving different events. It seems that, at least in certain cases, two-dimensional (2D) IR correlation spectroscopy can provide answers to this question. To illustrate this point, we have turned to the Cro-V55C dimer of the λ Cro repressor, a protein known to undergo thermal unfolding in two discrete steps through a stable equilibrium intermediate. The secondary structure of this intermediate is compatible with that of a partially unfolded protein and involves a reorganization of the N terminus, whereas the antiparallel β-ribbon formed by the C-terminal part of each subunit remains largely intact. To establish whether the unfolding process involves sequential events, we have performed a 2D correlation analysis of IR spectra recorded over the temperature range of 20–95°C. The 2D IR correlation analysis indeed provides evidence for a sequential formation of the stable intermediate, which is created in three (closely related) steps. A first step entails the unfolding of the short N-terminal β-strand, followed by the unfolding of the α-helices in a second step, and the third step comprises the reorganization of the remaining β-sheet and of some unordered segments in the protein. The complete unfolding of the stable intermediate at higher temperatures also undergoes sequential events that ultimately end with the breaking of the H bonds between the two β-strands at the dimer interface.

Conformational influences of glycosylation of a peptide: A possible model for the effect of glycosylation on the rate of protein folding

Improved strategies for synthesis make it possible to expand the range of glycopeptides available for detailed conformational studies. The glycopeptide 1 was synthesized using a new solid phase synthesis of carbohydrates and a convergent coupling to peptide followed by deprotection. Its conformational properties were subjected to NMR analysis and compared with a control peptide 2 prepared by conventional solid phase methods. Whereas peptide 2 fails to manifest any appreciable secondary structure, the glycopeptide 1 does show considerable conformational bias suggestive of an equilibrium between an ordered and a random state. The implications of this ordering effect for the larger issue of protein folding are considered.

Antibody C219 recognizes an α-helical epitope on P-glycoprotein

The ABC transporter, P-glycoprotein, is an integral membrane protein that mediates the ATP-driven efflux of drugs from multidrug-resistant cancer and HIV-infected cells. Anti-P-glycoprotein antibody C219 binds to both of the ATP-binding regions of P-glycoprotein and has been shown to inhibit its ATPase activity and drug binding capacity. C219 has been widely used in a clinical setting as a tumor marker, but recent observations of cross-reactivity with other proteins, including the c-erbB2 protein in breast cancer cells, impose potential limitations in detecting P-glycoprotein. We have determined the crystal structure at a resolution of 2.4 Å of the variable fragment of C219 in complex with an epitope peptide derived from the nucleotide binding domain of P-glycoprotein. The 14-residue peptide adopts an amphipathic α-helical conformation, a secondary structure not previously observed in structures of antibody–peptide complexes. Together with available biochemical data, the crystal structure of the C219-peptide complex indicates the molecular basis of the cross-reactivity of C219 with non-multidrug resistance-associated proteins. Alignment of the C219 epitope with the recent crystal structure of the ATP-binding subunit of histidine permease suggests a structural basis for the inhibition of the ATP and drug binding capacity of P-glycoprotein by C219. The results provide a rationale for the development of C219 mutants with improved specificity and affinity that could be useful in antibody-based P-glycoprotein detection and therapy in multidrug resistant cancers.

Localization of the N terminus of hepatitis B virus capsid protein by peptide-based difference mapping from cryoelectron microscopy

Recently, cryoelectron microscopy of isolated macromolecular complexes has advanced to resolutions below 10 Å, enabling direct visualization of α-helical secondary structure. To help correlate such density maps with the amino acid sequences of the component proteins, we advocate peptide-based difference mapping, i.e., insertion of peptides, ≈10 residues long, at targeted points in the sequence and visualization of these peptides as bulk labels in cryoelectron microscopy-derived difference maps. As proof of principle, we have appended an extraneous octapeptide at the N terminus of hepatitis B virus capsid protein and determined its location on the capsid surface by difference imaging at 11 Å resolution. Hepatitis B virus capsids are icosahedral particles, ≈300 Å in diameter, made up of T-shaped dimers (subunit Mr, 16–21 kDa, depending on construct). The stems of the Ts protrude outward as spikes, whereas the crosspieces pack to form the contiguous shell. The two N termini per dimer reside on either side of the spike-stem, at the level at which it enters the shell. This location is consistent with formation of the known intramolecular disulfide bond between the cysteines at positions 61 and −7 (in the residual propeptide) in the “e-antigen” form of the capsid protein and has implications for why this clinically important antigen remains unassembled in vivo.

The increased level of β1,4-galactosyltransferase required for lactose biosynthesis is achieved in part by translational control

β1,4-Galactosyltransferase (β4GalT-I) participates in both glycoconjugate biosynthesis (ubiquitous activity) and lactose biosynthesis (mammary gland-specific activity). In somatic tissues, transcription of the mammalian β4GalT-I gene results in a 4.1-kb mRNA and a 3.9-kb mRNA as a consequence of initiation at two start sites separated by ≈200 bp. In the mammary gland, coincident with the increased β4GalT-I enzyme level (≈50-fold) required for lactose biosynthesis, there is a switch from the 4.1-kb start site to the preferential use of the 3.9-kb start site, which is governed by a stronger tissue-restricted promoter. The use of the 3.9-kb start site results in a β4GalT-I transcript in which the 5′- untranslated region (UTR) has been truncated from ≈175 nt to ≈28 nt. The 5′-UTR of the 4.1-kb transcript [UTR(4.1)] is predicted to contain extensive secondary structure, a feature previously shown to reduce translational efficiency of an mRNA. In contrast, the 5′-UTR of the 3.9-kb mRNA [UTR(3.9)] lacks extensive secondary structure; thus, this transcript is predicted to be more efficiently translated relative to the 4.1-kb mRNA. To test this prediction, constructs were assembled in which the respective 5′-UTRs were fused to the luciferase-coding sequence and enzyme levels were determined after translation in vitro and in vivo. The luciferase mRNA containing the truncated UTR(3.9) was translated more efficiently both in vitro (≈14-fold) and in vivo (3- to 5-fold) relative to the luciferase mRNA containing the UTR(4.1). Consequently, in addition to control at the transcriptional level, β4GalT-I enzyme levels are further augmented in the lactating mammary gland as a result of translational control.

A transcriptional activating region with two contrasting modes of protein interaction

A C-terminal segment of the yeast activator Gal4 manifests two functions: When tethered to DNA, it elicits gene activation, and it binds the inhibitor Gal80. Here we examine the effects on these two functions of cysteine and proline substitutions. We find that, although certain cysteine substitutions diminish interaction with Gal80, those substitutions have little effect on the activating function in vivo and interaction with TATA box-binding protein (TBP) in vitro. Proline substitutions introduced near residues critical for Gal80 binding abolish that interaction but once again have no effect on the activating function. Crosslinking experiments show that a defined position in the activating peptide is in close proximity to TBP and Gal80 in the two separate reactions and show that binding of the inhibitor blocks binding to TBP. Thus, the same stretch of amino acids are involved in two quite different protein–protein interactions: binding to Gal80, which depends on a precise sequence and the formation of a defined secondary structure, or interactions with the transcriptional machinery in vivo, which are not impaired by perturbations of either sequence or structure.

Intron self-complementarity enforces exon inclusion in a yeast pre-mRNA

Skipping of internal exons during removal of introns from pre-mRNA must be avoided for proper expression of most eukaryotic genes. Despite significant understanding of the mechanics of intron removal, mechanisms that ensure inclusion of internal exons in multi-intron pre-mRNAs remain mysterious. Using a natural two-intron yeast gene, we have identified distinct RNA–RNA complementarities within each intron that prevent exon skipping and ensure inclusion of internal exons. We show that these complementarities are positioned to act as intron identity elements, bringing together only the appropriate 5′ splice sites and branchpoints. Destroying either intron self-complementarity allows exon skipping to occur, and restoring the complementarity using compensatory mutations rescues exon inclusion, indicating that the elements act through formation of RNA secondary structure. Introducing new pairing potential between regions near the 5′ splice site of intron 1 and the branchpoint of intron 2 dramatically enhances exon skipping. Similar elements identified in single intron yeast genes contribute to splicing efficiency. Our results illustrate how intron secondary structure serves to coordinate splice site pairing and enforce exon inclusion. We suggest that similar elements in vertebrate genes could assist in the splicing of very large introns and in the evolution of alternative splicing.

Characterization of the KIF3C Neural Kinesin-like Motor from Mouse

Proteins of the kinesin superfamily define a class of microtubule-dependent motors that play crucial roles in cell division and intracellular transport. To study the molecular mechanism of axonal transport, a cDNA encoding a new kinesin-like protein called KIF3C was cloned from a mouse brain cDNA library. Sequence and secondary structure analysis revealed that KIF3C is a member of the KIF3 family. In contrast to KIF3A and KIF3B, Northern and Western analysis indicated that KIF3C expression is highly enriched in neural tissues such as brain, spinal cord, and retina. When anti-KIF3C antibodies were used to stain the cerebellum, the strongest signal came from the cell bodies and dendrites of Purkinje cells. In retina, anti-KIF3C mainly stains the ganglion cells. Immunolocalization showed that the KIF3C motor in spinal cord and sciatic nerve is mainly localized in cytoplasm. In spinal cord, the KIF3C staining was punctate; double labeling with anti-giantin and anti-KIF3C showed a clear concentration of the motor protein in the Golgi complex. Staining of ligated sciatic nerves demonstrated that the KIF3C motor accumulated at the proximal side of the ligated nerve, which suggests that KIF3C is an anterograde motor. Immunoprecipitation experiments revealed that KIF3C and KIF3A, but not KIF3B, were coprecipitated. These data, combined with previous data from other labs, indicate that KIF3C and KIF3B are “variable” subunits that associate with a common KIF3A subunit, but not with each other. Together these results suggest that KIF3 family members combinatorially associate to power anterograde axonal transport.

Recombinant scrapie-like prion protein of 106 amino acids is soluble

The N terminus of the scrapie isoform of prion protein (PrPSc) can be truncated without loss of scrapie infectivity and, correspondingly, the truncation of the N terminus of the cellular isoform, PrPC, still permits conversion into PrPSc. To assess whether additional segments of the PrP molecule can be deleted, we previously removed regions of putative secondary structure in PrPC; in the present study we found that deletion of each of the four predicted helices prevented PrPSc formation, as did deletion of the stop transfer effector region and the C178A mutation. Removal of a 36-residue loop between helices 2 and 3 did not prevent formation of protease-resistant PrP; the resulting scrapie-like protein, designated PrPSc106, contained 106 residues after cleavage of an N-terminal signal peptide and a C-terminal sequence for glycolipid anchor addition. Addition of the detergent Sarkosyl to cell lysates solubilized PrPSc106, which retained resistance to digestion by proteinase K. These results suggest that all the regions of proposed secondary structure in PrP are required for PrPSc formation, as is the disulfide bond stabilizing helices 3 and 4. The discovery of PrPSc106 should facilitate structural studies of PrPSc, investigations of the mechanism of PrPSc formation, and the production of PrPSc-specific antibodies.

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.