Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Structure and function of a human transcription factor TFIIIB subunit that is evolutionarily conserved and contains both TFIIB- and high-mobility-group protein 2-related domains.

Transcription factor TFIIIB plays a central role in transcription initiation by RNA polymerase III on genes encoding tRNA, 5S rRNA, and other small structural RNAs. We report the purification of a human TFIIIB-derived complex containing only the TATA-binding polypeptide (TBP) and a 90-kDa subunit (TFIIIB90) and the isolation of a cDNA clone encoding the 90-kDa subunit. The N-terminal half of TFIIIB90 exhibits sequence similarity to the yeast TFIIIB70 (BRF) and the class II transcription factor TFIIB and interacts weakly with TBP. The C-terminal half of TFIIIB90 contains a high-mobility-group protein 2 (HMG2)-related domain and interacts strongly with TBP. Recombinant TFIIIB90 plus recombinant human TBP substitute for human TFIIIB in a complementation assay for transcription of 5S, tRNA, and VA1 RNA genes, and both the TFIIB-related domain and the HMG2-related domain are required for this activity. TFIIIB90 is also required for transcription of human 7SK and U6 RNA genes by RNA polymerase III, but apparently within a complex distinct from the TBP/TFIIIB90 complex.

Xeroderma pigmentosum and trichothiodystrophy are associated with different mutations in the XPD (ERCC2) repair/transcription gene

The xeroderma pigmentosum group D (XPD) protein has a dual function, both in nucleotide excision repair of DNA damage and in basal transcription. Mutations in the XPD gene can result in three distinct clinical phenotypes, XP, trichothiodystrophy (TTD), and XP with Cockayne syndrome. To determine if the clinical phenotypes of XP and TTD can be attributed to the sites of the mutations, we have identified the mutations in a large group of TTD and XP-D patients. Most sites of mutations differed between XP and TTD, but there are three sites at which the same mutation is found in XP and TTD patients. Since the corresponding patients were all compound heterozygotes with different mutations in the two alleles, the alleles were tested separately in a yeast complementation assay. The mutations which are found in both XP and TTD patients behaved as null alleles, suggesting that the disease phenotype was determined by the other allele. If we eliminate the null mutations, the remaining mutagenic pattern is consistent with the site of the mutation determining the phenotype.

Substantial narrowing of the Niemann–Pick C candidate interval by yeast artificial chromosome complementation

Niemann–Pick disease type C (NP-C) is an autosomal recessive lipidosis linked to chromosome 18q11–12, characterized by lysosomal accumulation of unesterified cholesterol and delayed induction of cholesterol-mediated homeostatic responses. This cellular phenotype is identifiable cytologically by filipin staining and biochemically by measurement of low-density lipoprotein-derived cholesterol esterification. The mutant Chinese hamster ovary cell line (CT60), which displays the NP-C cellular phenotype, was used as the recipient for a complementation assay after somatic cell fusions with normal and NP-C murine cells suggested that this Chinese hamster ovary cell line carries an alteration(s) in the hamster homolog(s) of NP-C. To narrow rapidly the candidate interval for NP-C, three overlapping yeast artificial chromosomes (YACs) spanning the 1 centimorgan human NP-C interval were introduced stably into CT60 cells and analyzed for correction of the cellular phenotype. Only YAC 911D5 complemented the NP-C phenotype, as evidenced by cytological and biochemical analyses, whereas no complementation was obtained from the other two YACs within the interval or from a YAC derived from chromosome 7. Fluorescent in situ hybridization indicated that YAC 911D5 was integrated at a single site per CT60 genome. These data substantially narrow the NP-C critical interval and should greatly simplify the identification of the gene responsible in mouse and man. This is the first demonstration of YAC complementation as a valuable adjunct strategy for positional cloning of a human gene.

Structure and function of the C-terminal PABC domain of human poly(A)-binding protein

We have determined the solution structure of the C-terminal quarter of human poly(A)-binding protein (hPABP). The protein fragment contains a protein domain, PABC [for poly(A)-binding protein C-terminal domain], which is also found associated with the HECT family of ubiquitin ligases. By using peptides derived from PABP interacting protein (Paip) 1, Paip2, and eRF3, we show that PABC functions as a peptide binding domain. We use chemical shift perturbation analysis to identify the peptide binding site in PABC and the major elements involved in peptide recognition. From comparative sequence analysis of PABC-binding peptides, we formulate a preliminary PABC consensus sequence and identify human ataxin-2, the protein responsible for type 2 spinocerebellar ataxia (SCA2), as a potential PABC ligand.

Ccm1, a regulatory gene controlling the induction of a carbon-concentrating mechanism in Chlamydomonas reinhardtii by sensing CO2 availability

Aquatic photosynthetic organisms, including the green alga Chlamydomonas reinhardtii, induce a set of genes for a carbon-concentrating mechanism (CCM) to acclimate to CO2-limiting conditions. This acclimation is modulated by some mechanisms in the cell to sense CO2 availability. Previously, a high-CO2-requiring mutant C16 defective in an induction of the CCM was isolated from C. reinhardtii by gene tagging. By using this pleiotropic mutant, we isolated a nuclear regulatory gene, Ccm1, encoding a 699-aa hydrophilic protein with a putative zinc-finger motif in its N-terminal region and a Gln repeat characteristic of transcriptional activators. Introduction of Ccm1 into this mutant restored an active carbon transport through the CCM, development of a pyrenoid structure in the chloroplast, and induction of a set of CCM-related genes. That a 5,128-base Ccm1 transcript and also the translation product of 76 kDa were detected in both high- and low-CO2 conditions suggests that CCM1 might be modified posttranslationally. These data indicate that Ccm1 is essential to control the induction of CCM by sensing CO2 availability in Chlamydomonas cells. In addition, complementation assay and identification of the mutation site of another pleiotropic mutant, cia5, revealed that His-54 within the putative zinc-finger motif of the CCM1 is crucial to its regulatory function.

Production of infectious recombinant Moloney murine leukemia virus particles in BHK cells using Semliki Forest virus-derived RNA expression vectors.

We describe a heterologous, Semliki Forest virus (SFV)-driven packaging system for the production of infectious recombinant Moloney murine leukemia virus particles. The gag-pol and env genes, as well as a recombinant retrovirus genome (LTR-psi (+)-neoR-LTR), were inserted into individual SFV1 expression plasmids. Replication-competent RNAs were transcribed in vitro and introduced into the cytoplasm of BHK-21 cells using electroporation. The expressed Moloney murine leukemia virus structural proteins produced extracellular virus-like particles. In these particles the gag precursor was processed into mature products, indicating that the particles contained an active protease. The protease of the gag-pol fusion protein was also shown to be active in a trans-complementation assay using a large excess of Pr65gag. Moreover, the particles possessed reverse transcriptase (RT) activity as measured in an in vitro assay. Cotransfection of BHK-21 cells by all three SFV1 constructs resulted in the production of transduction-competent particles at 4 x 10(6) colony-forming units (cfu)/ml during a 5-hr incubation period. Altogether, 2.9 x 10(7) transduction-competent particles were obtained from about 4 x 10(6) transfected cells. Thus, this system represents the first RNA-based packaging system for the production of infectious retroviral particles. The facts that no helper virus could be detected in the virus stocks and that particles carrying the amphotropic envelope could be produced with similar efficiency as those that carry the ecotropic envelope make the system very interesting for gene therapy.

GRIP1, a novel mouse protein that serves as a transcriptional coactivator in yeast for the hormone binding domains of steroid receptors.

The yeast two-hybrid system was used to isolate a clone from a 17-day-old mouse embryo cDNA library that codes for a novel 812-aa long protein fragment, glucocorticoid receptor-interacting protein 1 (GRIP1), that can interact with the hormone binding domain (HBD) of the glucocorticoid receptor. In the yeast two-hybrid system and in vitro, GRIP1 interacted with the HBDs of the glucocorticoid, estrogen, and androgen receptors in a hormone-regulated manner. When fused to the DNA binding domain of a heterologous protein, the GRIP1 fragment activated a reporter gene containing a suitable enhancer site in yeast cells and in mammalian cells, indicating that GRIP1 contains a transcriptional activation domain. Overexpression of the GRIP1 fragment in mammalian cells interfered with hormone-regulated expression of mouse mammary tumor virus-chloramphenicol acetyltransferase gene and constitutive expression of cytomegalovirus-beta-galactosidase reporter gene, but not constitutive expression from a tRNA gene promoter. This selective squelching activity suggests that GRIM can interact with an essential component of the RNA polymerase II transcription machinery. Finally, while a steroid receptor HBD fused with a GAL4 DNA binding domain did not, by itself, activate transcription of a reporter gene in yeast, coexpression of this fusion protein with GRIP1 strongly activated the reporter gene. Thus, in yeast, GRIP1 can serve as a coactivator, potentiating the transactivation functions in steroid receptor HBDs, possibly by acting as a bridge between HBDs of the receptors and the basal transcription machinery.

Transdominant genetic analysis of a growth control pathway

Genetic selections that use proteinaceous transdominant inhibitors encoded by DNA libraries to cause mutant phenocopies may facilitate genetic analysis in traditionally nongenetic organisms. We performed a selection for random short peptides and larger protein fragments (collectively termed “perturbagens”) that inhibit the yeast pheromone response pathway. Peptide and protein fragment perturbagens that permit cell division in the presence of pheromone were recovered. Two perturbagens were derived from proteins required for pheromone response, and an additional two were derived from proteins that may negatively influence the pheromone response pathway. Furthermore, three known components of the pathway were identified as probable perturbagen targets based on physical interaction assays. Thus, by selection for transdominant inhibitors of pheromone response, multiple pathway components were identified either directly as gene fragments or indirectly as the likely targets of specific perturbagens. These results, combined with the results of previous work [Holzmayer, T. A., Pestov, D. G. & Roninson, I. B. (1992) Nucl. Acids. Res. 20, 711–717; Whiteway, M., Dignard, D. & Thomas, D. Y. (1992) Proc. Natl. Acad. Sci. USA 89, 9410–9414; and Gudkov, A. V., Kazarov, A. R., Thimmapaya, R., Axenovich, S. A., Mazo, I. A. & Roninson, I. B. (1994) Proc. Natl. Acad. Sci. USA 91, 3744–3748], suggest that transdominant genetic analysis of the type described here will be broadly applicable.

The wing in yeast heat shock transcription factor (HSF) DNA-binding domain is required for full activity

The yeast heat shock transcription factor (HSF) belongs to the winged helix family of proteins. HSF binds DNA as a trimer, and additional trimers can bind DNA co-operatively. Unlike other winged helix–turn–helix proteins, HSF’s wing does not appear to contact DNA, as based on a previously solved crystal structure. Instead, the structure implies that the wing is involved in protein–protein interactions, possibly within a trimer or between adjacent trimers. To understand the function of the wing in the HSF DNA-binding domain, a Saccharomyces cerevisiae strain was created that expresses a wingless HSF protein. This strain grows normally at 30°C, but shows a decrease in reporter gene expression during constitutive and heat-shocked conditions. Removal of the wing does not affect the stability or trimeric nature of a protein fragment containing the DNA-binding and trimerization domains. Removal of the wing does result in a decrease in DNA-binding affinity. This defect was mainly observed in the ability to form the first trimer-bound complex, as the formation of larger complexes is unaffected by the deletion. Our results suggest that the wing is not involved in the highly co-operative nature of HSF binding, but may be important in stabilizing the first trimer bound to DNA.

Oligomerization domain-directed reassembly of active dihydrofolate reductase from rationally designed fragments

Reassembly of enzymes from peptide fragments has been used as a strategy for understanding the evolution, folding, and role of individual subdomains in catalysis and regulation of activity. We demonstrate an oligomerization-assisted enzyme reassembly strategy whereby fragments are covalently linked to independently folding and interacting domains whose interactions serve to promote efficient refolding and complementation of fragments, forming active enzyme. We show that active murine dihydrofolate reductase (E.C. 1.5.1.3) can be reassembled from complementary N- and C-terminal fragments when fused to homodimerizing GCN4 leucine zipper-forming sequences as well as heterodimerizing protein partners. Reassembly is detected by an in vivo selection assay in Escherichia coli and in vitro. The effects of mutations that disrupt fragment affinity or enzyme activity were assessed. The steady–state kinetic parameters for the reassembled mutant (Phe-31 → Ser) were determined; they are not significantly different from the full-length mutant. The strategy described here provides a general approach for protein dissection and domain swapping studies, with the capacity both for rapid in vivo screening as well as in vitro characterization. Further, the strategy suggests a simple in vivo enzyme-based detection system for protein–protein interactions, which we illustrate with two examples: ras–GTPase and raf–ras-binding domain and FK506-binding protein-rapamycin complexed with the target of rapamycin TOR2.

Solution structure of a 142-residue recombinant prion protein corresponding to the infectious fragment of the scrapie isoform

The scrapie prion protein (PrPSc) is the major, and possibly the only, component of the infectious prion; it is generated from the cellular isoform (PrPC) by a conformational change. N-terminal truncation of PrPSc by limited proteolysis produces a protein of ≈142 residues designated PrP 27–30, which retains infectivity. A recombinant protein (rPrP) corresponding to Syrian hamster PrP 27–30 was expressed in Escherichia coli and purified. After refolding rPrP into an α-helical form resembling PrPC, the structure was solved by multidimensional heteronuclear NMR, revealing many structural features of rPrP that were not found in two shorter PrP fragments studied previously. Extensive side-chain interactions for residues 113–125 characterize a hydrophobic cluster, which packs against an irregular β-sheet, whereas residues 90–112 exhibit little defined structure. Although identifiable secondary structure is largely lacking in the N terminus of rPrP, paradoxically this N terminus increases the amount of secondary structure in the remainder of rPrP. The surface of a long helix (residues 200–227) and a structured loop (residues 165–171) form a discontinuous epitope for binding of a protein that facilitates PrPSc formation. Polymorphic residues within this epitope seem to modulate susceptibility of sheep and humans to prion disease. Conformational heterogeneity of rPrP at the N terminus may be key to the transformation of PrPC into PrPSc, whereas the discontinuous epitope near the C terminus controls this transition.

Crystal structure of a Staphylococcus aureus protein A domain complexed with the Fab fragment of a human IgM antibody: Structural basis for recognition of B-cell receptors and superantigen activity

Staphylococcus aureus produces a virulence factor, protein A (SpA), that contains five homologous Ig-binding domains. The interactions of SpA with the Fab region of membrane-anchored Igs can stimulate a large fraction of B cells, contributing to lymphocyte clonal selection. To understand the molecular basis for this activity, we have solved the crystal structure of the complex between domain D of SpA and the Fab fragment of a human IgM antibody to 2.7-Å resolution. In the complex, helices II and III of domain D interact with the variable region of the Fab heavy chain (VH) through framework residues, without the involvement of the hypervariable regions implicated in antigen recognition. The contact residues are highly conserved in human VH3 antibodies but not in other families. The contact residues from domain D also are conserved among all SpA Ig-binding domains, suggesting that each could bind in a similar manner. Features of this interaction parallel those reported for staphylococcal enterotoxins that are superantigens for many T cells. The structural homology between Ig VH regions and the T-cell receptor Vβ regions facilitates their comparison, and both types of interactions involve lymphocyte receptor surface remote from the antigen binding site. However, T-cell superantigens reportedly interact through hydrogen bonds with T-cell receptor Vβ backbone atoms in a primary sequence-independent manner, whereas SpA relies on a sequence-restricted conformational binding with residue side chains, suggesting that this common bacterial pathogen has adopted distinct molecular recognition strategies for affecting large sets of B and T lymphocytes.

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.

The crystal structure of the nitrogen regulation fragment of the yeast prion protein Ure2p

The yeast nonchromosomal gene [URE3] is due to a prion form of the nitrogen regulatory protein Ure2p. It is a negative regulator of nitrogen catabolism and acts by inhibiting the transcription factor Gln3p. Ure2p residues 1–80 are necessary for prion generation and propagation. The C-terminal fragment retains nitrogen regulatory activity, albeit somewhat less efficiently than the full-length protein, and it also lowers the frequency of prion generation. The crystal structure of this C-terminal fragment, Ure2p(97–354), at 2.3 Å resolution is described here. It adopts the same fold as the glutathione S-transferase superfamily, consistent with their sequence similarity. However, Ure2p(97–354) lacks a properly positioned catalytic residue that is required for S-transferase activity. Residues within this regulatory fragment that have been indicated by mutational studies to influence prion generation have been mapped onto the three-dimensional structure, and possible implications for prion activity are discussed.